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~The Mechanics of Non-Linear Materials with Memory — 


Pant 


A.E.GREEN & R.S,. RIVLIN 


1. Introduction 


In recent years a number of writers have formulated rheological equations 
of state for non-linear visco-elastic materials. A review of the subject up to 
1952 has been given by TRUESDELL (1952, 1953). RrivLin & ERICKSEN (1955) 

and RIvLINn (1955) have proposed a theory in which the components of stress 
at time ¢# in an element of material depend only on the gradients of displacement, 
velocity, acceleration, second acceleration, ..., (#—1)™ acceleration in that 

element at time ¢#. Some special problems in this theory have been solved by 

_ RIVvLIN (1956). 

Other aspects of the mechanics of non-linear materials have been considered 
by No iy (1955), particularly with reference to the invariance requirements for 
“hygrosteric’’ materials and various special cases of these materials have been 
examined, including the theory of “hypo-elastic’”’ materials proposed by TRUEs- 
DELL (1955). The relation between the theories of hypo-elasticity and plasticity 

_ has been studied by GREEN (1956a, b). 

OLDROYD (1950) suggested that the right invariance properties which must 

~ be satisfied by a rheological equation of state can most readily be recognised 
if the frame of reference is a coordinate system convected with the material. 
Using this procedure he formulated rheological equations of state of various 
types for visco-elastic materials. 

In the present paper we adopt a procedure which is akin to that employed 
by Rrviin & ERICKSEN (1955) in their work on the stress-deformation relations 
for isotropic materials. In § 2 we assume that the stress in a rectangular Cartesian 
coordinate system at time ¢ depends on the displacement gradients at time ¢ 
and at N previous instants of time in the interval 0 to ¢. We then determine 
the limitations imposed on the stress-deformation relation by the restriction that 
if the body is subjected to a rigid rotation and the reference system is rotated 
with the body, the stress is unaltered by the rigid-body rotation. In § 3, the further 
restrictions imposed on the stress-deformation relation by isotropy of the material 
are derived. 

In § 4, we return to the case discussed in § 2 and allow the number N to 
tend to infinity so that certain functions in the stress-deformation relations are 
replaced by functionals. It is then shown that these functionals may be expressed 
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2 A. E. GREEN & R. S. RIVLIN: 
to any desired approximation by the sum of a number of multiple integrals. | 
In §5 the limitations imposed on the form of these integrals by isotropy of the 
material are discussed. 

In §6, we derive sufficient conditions on these integrals for the stress- 
deformation to be of the hereditary type and in §§ 7 and 8 we discuss the relation 
of the stress-deformation relations so derived to those of the theory of Rrvi1n & 
ERICKSEN. Finally, in § 9 some of the equations derived in the present paper are 
transformed to a general curvilinear coordinate system convected with the 
material. 

2. Stress Dependent on Displacement Gradients 


Let every particle of a three-dimensional body be at rest, prior to time T=0, 
relative to a fixed rectangular Cartesian coordinate system x;. Let us suppose 
that after zero time the body is deformed continuously so that a typical particle 
P, of the body moves to P at time ¢ and that at any time t (0S 1S/) the co- 
ordinates of the particle in the system x; are x;(t). We shall employ the notation 


X= %,(0) ~ and 0) (2.4) 


We consider that the functions x,;(t) are specified as single-valued functions of 
X,and 7, so that the deformation of the body is completely defined. We then have 


x A(t) =%,(X,,t) and %*,=%x,(X;,2). (2.2) 


If this deformation is to be possible in a real material then |éx;(t)/¢X, >0. 


We dencte the stress components in the coordinate system x; at the particle 
F, and time t by o;,;(t) and employ the notation 


ig = 0;;(2). (2.3) 


We consider that the components of stress o,; are polynomial functions of 
the 9(N +4) displacement gradients 0x,(t,)/OX, («=0,1,2,...,N) at N+4 


oO 


distinct instants of time t,, T:,..., Ty, T)(=#) between t=O and t=¢#. Thus, 
aes (ae) 

OF Or i oe 4 2.4 

j= 045 | Aa (2.4) 

where «=0,1,2,...,N and o;;=o;;. Now, let us suppose that we superpose 


on the deformation described by (2.2) an arbitrary rigid body rotation. Then, 
the coordinates %;(t) of the particle Py at time t in the coordinate system x, 
will be given by 


Hares (a) (0) (2.5) 


where A,,(t) is any continuous function of t defining the rigid body rotation 
and is therefore subject to the conditions 


A y(t) A), («) =A AU) Ay (tne [Ar (e) tt (2.6) 


The stress components g,; at time ¢ in the coordinate system x; due to the 
resultant of the deformations described by (2.2) and (2.5) are given, from (2.4), by 


Oe 


ax, 


vy vq 


(2.7) 
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Let x;*(z) be a moving rectangular Cartesian coordinate system which coincides 
with the fixed system x; when t=0 and rotates with the rigid body motion 
described by (2.5). The stress components at time ¢ in this system corresponding 
to the resultant of the deformations described by (2.2) and (2.5) are, of course, 


o;;. At time ¢, the coordinate systems x;() and x; are related by 
%, =A, x; (), Aiy=Ai,(t). (2.8) 
Hence, 
Fi5= AipAjg Opa. (2.9) 


Equations (2.4), (2.7) and (2.9) yield 


0%, (Tq) 


Oxga(a 
os HL (2) 


ax, 


(2.10) 


ni Ap Aja pa| 


The arguments in this relation may be considered to be the components of 
3(N-+1) vectors g,(t,), the components of which are 0x,(t,)/0X, in the co- 
ordinate system x%;(t,) and 0%,(t,)/@X, in a coordinate system %;(t,) related to 
%;(T) by (2.5) with t=7,. Taking A;,;(t)=6,; in (2.10) for all t except in a 
small interval about t,, equation (2.10) becomes 


er ae Ga we OF AlN ey 0): 
UN OS DX ee OX 
— o,,| 2% la) 8%, (%) |, 8% (tw) 8%, (0) (2.41) 
GATT OTE Dow Ee? Epp tae De alee ae 


where %,(t,) and x,(t,) are related by (2.5) with t=7,, for all A;;(t,) satisfying 
(2.6) with t=1,. Hence o;; must be form-invariant under all proper orthogonal 
transformations of rectangular Cartesian coordinate axes %;(t,), the axes %,(T,) 
(«=0, 2,..., N) remaining unaltered. Similarly, it can be shown that o;; must 
be form-invariant under all proper orthogonal transformations of rectangular 
axes x;(t,) for each value of t, (x=1,2,..., .N). The transformations for each 
value of « are independent of those for the remaining values of «. A complete 
‘table of typical basic invariants for the orthogonal group is known from classical 
invariant theory [see WEYL (1939) p. 53]. Every even* orthogonal polynomial 
invariant depending on three vectors g,(t,) (for a given «) in three-dimensional 
vector space is expressible, in polynomial form, in terms of six scalar products 
G,(Tx) Js (Tt). Every odd invariant is the sum of terms [g,(Tq) Gs (Tx) Gi (Ta) If, 
where / is an even invariant and [g,(t.) g; (Tx) Gi(T.) | is the scalar triple product. 
The basic invariants are therefore 


Y, (a ‘9s (T..) are er (t) = nine ae . (2.12) 


and 


[9, (Ta) gs (T.) 9: (Ta) | = rst Ve (ta) = rst | Ox 


foten 41.2, &-, N}) where 
(2.14) 


& (Ta) = [Sm n (Ta) 


, 


* An even invariant is also invariant under improper orthogonal transformations 
but an odd invariant is not. The theorem applies to invariants which are polynomials 
in the components of the vectors. 


4% 
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and e,,, is the alternating Cartesian tensor. We thus have 
ATER VOMIT ee 
04; = F4j|8rs(%), 8,5 (Te); =") Les (EN)> Vga), g(t), — Vg (tw), yal , (2.15) 


the polynomial dependence of o;; on the arguments in (2.15) being, in general, 
different from that in equations (2.4) to (2.11). From (2.15) and (2.10), we have 


O57 lers(to) Ve (eye Sa fo) 


=A,,A4; Apa {Bro ta ), Vg (ta), oe a]. (2.16) 


Let a,b be two arbitrary vectors with components a;,}; in the coordinate 
system x;. We form the quantity 


F= a;0;0;;. (2417) 
Then 
oer 
64 = oa (2.18) 


It follows from (2.16) that F is a polynomial scalar invariant * of the three vectors 
g;(t)) and a and b for proper orthogonal transformations. F can therefore be 
expressed as a polynomial in the quantities defined by (2.12) and (2.13) with 
a0 and the further invariantsa-a,b-bandJ,, J,, P.,, Q,,, K and R, defined 


i L490) as eg 
: * TON ‘ taxp? 
P, = [ag,(0) 9.()] = eimai gy ox 
Ors = [B Gr) G.(6)] = enn dj ae Ge» 
K =a-b=a;);, R,=(abg,()] = ¢,;,4;0, on (2.19) 


Since F is homogeneous and bilinear in a and b, it must be expressible in 
the form 


Ate 3 L,Js+@K+ C pars i Pq Coan G, R, a Forty On oe <i fo rs? (2.20) 
where the @’s, ¥’s and Y’s are polynomials in the invariants defined by equations 
(2.12) and (2.13) with «=0,1,2,..., N. 

It follows from re and (2.20) that 

Ox; OtpALOX, OX, OF, 
=e x + 8544 aes ikl aX, EX, timn axe ax, + 


CX FOR OG 
+, 65:4 = ae + Fp ps jmn 2x, ax, ax. t 


oO 


ay 


0%; OX%m, OX» 


+ Gprslimn OX, OX, OX, ° (2.24) 
Now, employing the relation 
bi, Oe Ont 
Cnt Can = | On jr Onms Onn | s (2:22) 
Ors Cbyns Ot | 


* 
i a F * a scalar polynomial which is form-invariant regarded as a function of 
s\To), &, 
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we obtain, with the notation i; (4) =8;;; 


Chas (aeign Chee, aes 
Cinleimn OX, OX, Ox AIOx S107 (Oy Patan hos Core 


Ox; Ox; OX; 0x; 
T OX, aX, 8s OX, OX, bas | 
0%, OX; OX; OX; 
SHC» ade WD Maes) YD 
=f OX, OX » Ear Ox, 0X4 Spr: (2.23) 


Again, using the relation (2.22), we obtain 


OM, VON Oe == 
Cimn aX, ax, ax. Ve @ 


OO On, 


tO Mia OKs CON 


SG. é Sa ie ~ 
jmnCuvw @X, OX, OX, OX, OX, OX, 
OX; 0%; 
ee OX, aX, CrimS&ir8ms ° (2.24) 
Whence, OX; O%m O%, 0%; Ox; 1 3 
Cimn OX» OX, OK, OX Ox Vs) CnimS&ir8ms- (2.25) 


Introducing the relations (2.23) and (2.25) into (2.21) and bearing in mind 
that o;; 1s symmetric, we obtain 
1 Ox, Ox; 
i Veo) Lax, axes be + Ou]: ia) 
where /,, and 7 are polynomials in g,,(t,,) and Vg (2) (head Oke Meee Oth) 

With the notation 


0; 


1 
Wed eae Opals (2.27) 
equation (2.26) may be re-written as 
1 £ Ox; Ox; | 
055 = ag | 8a (ax aX, oo (2.28) 


If the material considered is incompressible, Vg (t,) =1 and o;, is undeter- 
mined to the extent of an arbitrary hydrostatic pressure p, so that equation (2.26) 
becomes 


Ox; OX; 
Oi; = ax, aX, Irs — P93;5 (2.29) 


where /,, is a polynomial in g,,(t.) («=0,1,2,..., N). 


3. Isotropic Materials 


If the material considered is isotropic when ¢=0, equation (2.28) must be 
form-invariant for proper orthogonal transformations of the rectangular Cartesian 
coordinate system x,;. Let %; be a fixed rectangular Cartesian coordinate system 

elated to x; b = 
are constants satisfying the orthogonality conditions 


A,,4j,= Ay, Ay, = 0; 


where 4;; 


fat, (3.2) 


2? 


6 A. E. GREEN & R. S. RIVLIN: 


If %;(t) denotes the coordinates of x;(r) in the system x; and X;=%;(0), we 


have, from (3.1), 


X ; (7) = A;,;%;(t), X;=A;;X;, 

Ae x(t) O#,(T) 

= Pasa) Gee) . 
&;; (7) == ax, aX; A NEG (ay. (3 3) 


Denoting the stress in the system ¥; by o;;, the form-invariance of (2.28) requires 


that re jean yee APA 
O74 ne Vz \/ Leaulcraye V2 (T.) | i] ig 
OX; Ox; 135 3] = = 
=r i ae B Sys 5i,] hes [8pq(ta)» 8 (ta) Ip (3.4) 


From (3.4) and (2.28) and the relation 


61; = AjpAjn Onn (3.5) 


jn 


it follows that 


| a = oe 5:;) Irs leaa (Ta), Ve(e.)] 


OX, ax, 
j 2) m 0d n , tod heaps 
pai Ain Asn Be oe i = vs Smn) is [oq (Ta), | g(T.) | ’ 
and 
P*18p a (Ta) VE (Ta) | =F" Baltes Veal, (3.6) 
where g,,(t) and g,,(t) are related by (3.3) so that 


g(t) = g(t), (3.7) 
and the equations are valid for all A;, satisfying (3.2). 
Employing the notation 


9(t) =\ls,() |, g() =|lg.;(x)|| and A=||4,, 


; (3.8) 


the last of the relations (3.3) can be re-written as 


g(t) =Ag(t) A’, (3.9) 
where A’ denotes the transpose of A. 


With (3.7), the second of equations (3.6) states that /* is a polynomial*® in 
Vg (t.,.) and a polynomial scalar invariant under proper orthogonal transformations 
of the matrices g(t,). It may therefore be expressed as a polynomial in Veg (T) 
(7=0,1,..., N) and traces of products formed from the matrices O(te) (a=0, 
eae se 

Let b;; and 6;, be the components of an arbitrary (not necessarily symmetric) 


tensor in the coordinate systems x; and x; respectively. It is readily seen that 
B,, and B,,, defined by 


LVAD OF 1 
Bas Cal bee OX, &s 5:4) b; 5 


* We note that Ve (t,) is a scalar invariant of g(t) under proper orthogonal 
transformations but is not a polynomial scalar invariant. 
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and 


= Ox; OX; ics zt 
Dey Ps : 
= (See gee wy brs Oui) Bi (3-10) 
are the components of a Cartesian tensor in the coordinate systems x; and x; 
respectively. We obtain from the first of equations (3.6) and equation (3.7), the 
relation 


Berrle te ), Ve.) T)| = a Byles epukts Ve (3.11) 


Denoting the expression on the right-hand side of (3.11) by F and with the 
notation B=||B,,||, equation (3.11) states that F is a polynomial in Vg (t,,) and 
a polynomial scalar invariant of the matrices B and g(t,). It can therefore be 


expressed as a polynomial in Vg (T.) («=0,1,...,.N) and traces of products 
formed from the matrices B and g(t,), (c=0,1,...,.N). Since F is linear in 
the elements of B, /,, is given by 

fs = OF /0B,,. (3.12) 
We see that f/,, are the components of a symmetric matrix polynomial in the 
matrices g (t,) in which the coefficients are polynomials in \/g (a) (= OAs cea) 
and traces of products formed from the matrices g(t,) («=0,1,..., N). It is 


immediately apparent that g,,/,, 1s a scalar polynomial in Vg (t,,) and traces of 
products formed from the matrices g(t,). Hence, from (2.27), f is a scalar poly- 
nomial in \/g (t,,.) and traces of products formed from the matrices g(t,). 


We thus obtain the result that for the isotropic material 


5 [fur + ax aa hs (3.13) 


where /,, are the elements of a Seaeeey matrix polynomial in the matrices 
g(t) (x«=0,1,..., N) and / and the coefficients in /,, are scalar polynomials 
in Vg (t,) and traces of products formed from the matrices g (T,). 


Employing the notation 


O=|lon||}-— © = ||0x,/0X,|| ands f= ||7,,\|, (3.44) 
we may express the result (3.13) in matrix notation by 
1 
6 ——— |i efe], (3.15) 
Ve (2) 


where J denotes the unit matrix and ec’ denotes the transpose of €. 
If the material considered is incompressible, it is readily seen that equation 
(3.13) becomes 
oe Ox; 
=P OG ae ere Phe (3.16) 


where # is an arbitrary hydrostatic pressure and /,, are the elements of a sym- 
metric matrix polynomial in the matrices g (t,) (« =0, 1, 2, ..., N), the coefficients 
in which are scalar polynomials in traces of products formed from the matrices 


g (Tx) 
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4. Tensor Functionals 


We begin by taking an expression for the stress o;; at time ¢ in the form 


(2.28), vez. ) 
Cie Ve® i [pq(Ta) » ), Ve (eo) ia 
t (ae : , : a [ep 4 (ta), Ve (ea), (4.1) 
where 1, denotes N-+4 instants of time t, T:,.-., Tv, To(=#) and 7 and 7,, are 


polynomials in the arguments. Following ideas expounded by VOLTERRA (1930), . 
we pass from tensor functions of the form (4.1) to tensor functionals and write 


eel vaay | {pa ) eC 


ee ax —s 87s i fone [8pa(7), Ve} o (4.2) 


This equation means that /*, /,, and hence o; ; depend on all the values taken 
by gpq(t) and Vea) when t varies in the invest Os ec 2 

Bearing in ie that Vg@ (rz) is a single-valued function of the g,,(t), we see 
that there is no loss of generality in taking /* and /,, in (4.2) as functionals of 
the g,,(t) only and omitting the factor 14 Vg 

We thus obtain 


: : t 
54 = [Boa ()] 8; + (Se le — Bro Oss) brs [Bp oO). (4.3) 
0 x > 0 


7* and /,, denoting new functionals which are different, in general, from those 
which these symbols denote in (4.2). 

In the following discussion we shall assume that p(t) belongs to the aggre- 
gate (C) of functions which are continuous in 0<1t<¢ and that the aggregate 
is compact. Also, we assume that /* and /,, in (4.3) are continuous functionals 


on (C). We denote the Fourier half-range cosine coefficients of Spq(t) by Ge 
Then, 


t 
2 
IMO 
0 


and ; 
~ 7 J eels) de (4.4) 


Following an argument mits to that employed by Hee (1910), we Bee 
approximate the continuous functional /* by a polynomial P (GAG) eee) 


which tends uniformly to f* as «— oo. Now, let us consider s by nical seen 
Gla) G2) Ge) 


bit Fba**> Shpqp I the polynomial P™, where some or all of Oy 5 gs Hay Ole 
may be repeated values of the same integer. From (4.4), we obtain 

(1) (% 
CG ie: ora. ff Ls we = es re 008 ak 


x Sp.a, és Ep.4¢ cz) -++ Sprar (TR) A%4%,...dTR, (4.5) 
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where Cp=2"-?t-*, p denoting the number of the indices «,, %,...,% Which 
are zero. 
t 


It follows that /*[g,,(t)] may be expressed with any desired degree of 


0 
approximation in the form 
t t 
P [8pq(t)] = K() + JEpalt, T) &pq(t) dt + 
0 


t 
oh ID Se (é, T%)> To) 8p, n (7) Bhs q (72) at, 4%, =i oe =e 


00 
t 


ae 
FSS SB pape. baanll Ty, Te, oat OR) x 


0 


X 8,9, (%) Spode (T2) -++ Sppox (tr) 414% g-0.dTR. (4.6) 


t 
Similarly, it is seen that /,,{g,,(t)] may be expressed with any desired degree 
0 
of approximation in the form 
t t 
eset) = L,;(t) ae J Lispall T) &pq(T) dt ae 
0 


tt 
at LI Lr spataal qe (¢, Ty, Tp) Pena ea) S292 (T.) AT, + -- aos 


0 
ane t 
Se ae epee erate beaalts Lan hate one). x 
00 0 
Kaood (Pnerig lle) c- eae ag oR) 24 tec. A TR. (4.7) 


In (4.6) and (4.7) the K’s and L’s are continuous functions of their arguments. 


5. Isotropic Tensor Functionals 


If the material is isotropic, the expression (4.3) for the stress must be form- 


invariant under proper orthogonal transformations of the reference system. We 
t 


t 
have already remarked that the functionals /*[g,,(z)] and /,,[g,,(t)] may be 


9 0 
expressed with any desired degree of accuracy as polynomials in the quantities 
Ge GY, Sioa Ge). Let us write 


q 
t 
f* [8p q(t)] = P (G3?) 
0 


and 


t 
frs(8pq(7)] = PY (GY), (5.1) 
where the index yw takes the values 0,1, 2,...,« and Pe) is symmetric. Then, 
from (4.3), we have 
PONG) js |e Oe 1 (a) 
Gp POUCH) bys + [ae Ge — Ere us) BO (GH) (5.2) 


We now consider the conditions which must be satisfied by P@ and P® for 
(5.2) to be form-invariant under a proper orthogonal transformation of the 
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ployed is similar to that used in § 3, in which 


the case of stress dependent on the displacement gradients at a finite number 
of instants of time is considered. ¥; is a fixed rectangular Cartesian coordinate 


system related to x; by equation (3.1) and Z;;(t) is defined by equation (3.3). 
Analogously with (4.4), we write 


reference system. The method em 


t 
Ge) a +f &»_(T) cos a= at Ate) (5.3) 
0 
and t 
= 4 a. 
Go) = +f Boal) dt, 
0 
and note that Ge = AAG (5.4) 


The form-invariance of equation (5.2) requires that 


= a) (Glu OX; 0%; 15 x) ((H) 
Bij = PO CP) bu; + (Sot Gxt — Bre Bus) Bo CH)- (5.5) 


From equations (5.2), (5.5), (3.3) and (3.5), we obtain 


Nn Wral Ox, Ox. 1 (a) (4) 
rs 6:3] Tee (GY) =a Acalun SS Ox. a = Srs Om APs (Gyo) 


aut Po EY) = PoE). (5.6) 


The second of these equations implies that P® is a polynomial scalar invariant 
under proper orthogonal transformations of the matrices G™ (= ||G¥) ||), w=0, 
1,2,...,a. It may therefore be expressed as a polynomial in traces of products 
formed from these matrices. Again, we can employ the first of equations (5.6) 
in a manner similar to that used in § 3 to show that P®™ (= ||P ||) is a sym- 
metric isotropic matrix polynomial in the matrices GM. 

It is apparent that if, from the terms on the left-hand side of equation (4.5), 
we form a product of traces of products of the matrices GM, G@), ..., G2), 


then, in the integrand on the right-hand side of (4.5), we must replace 


8, q, (71) 8p, , (Te) ++» Sonar (TR) 


by a product of traces of products of the matrices g(t), g (ts), -.., g (Tr). For 
example, we have from (4.5) 


ty Gia) Gia) ty Gas) Gia) space Gia) 


see (41) (a2) (3) (4) (ar) 
a Gr, Pe Gy. Pi Gp, Ps Gp ee Cpr bs 


t ¢t t 
AX, UT Xo TT 
=f [++ [Cxeos i legos = 2... cos “REIR 
ONtor Meo 


X 8p. p, (7) Bp, p, (Ta) 8p, p, (Ts) 8p.p, (Ta) °° * Sprp, (Tr) 


bat t 
a, UT Ky UT. 
= [ ff Cxcos ze ‘cos - Bon COGE 
0 0 0 


t 


X g(t) G(T.) vg (ts) g (Ta) +++ G (TR) dt, dt,-+- dtp. (5.7) 
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In a similar way, we see that if, from the terms on the left-hand side of 
equation (4.5), we form a symmetric isotropic matrix polynomial in @™, G), 
..., GI”), then in the integrand on the right-hand side of (4.5) we must replace 
8p, 4, (T1) 8,9, (Ta) --+ pan (Tr) by a symmetric isotropic matrix polynomial in 
9 (71), G(T), ---, g(t) Which is multilinear in these matrices. For example, 
from a 


tr (Gi) Ge) [Glo Gow). Glen) + Gee), Gl) Gl] 


Oy Ot») 
— Eo) oe ‘ 


(%) (> (%s) 
On, Geos: Gos, ha Coe Ge 1PR Gy | 


eee 
et Ay Ty As UT, “Kare 
=| f- face el ae, -++ COS Boi, (Ty) Boy p, (Fa) X 
00 


x ie, (Ts) &p.95 (Ta) --- Song (Tr) + &p., (Tr) rire pple en Ua) Spey (Fs) 
GAT IA Te a oR 


= = ff [Cx Pecos wet e.. “++ COS oR SE x [vg (t%) g (T2)] X 


x ees + G (Tr) + G(T) ---G (Ta) G (ta) tA T2...dte. (5.8) 


It follows from the above considerations that P! (G )) may be expressed in 
the form 


t 
PO (GE) = 9 7! GY {t,t g(t) dt + 


tt 
se Ab T,T 3, G(T), 9 (T.)$ dt, dT. +--+ 


igo Oe 
Shee gee Tq, T3,+--, tr, J (T)s Ge) oJ (Te)} AH AT... dT, (5.9) 


where py is a homogeneous polynomial scalar invariant of the matrices g (t,), 
g(t.), .--, g(t) linear in each of them, and a function of ¢, 1, T),...,T,. Ina 
similar manner it follows that 


is 
PO (GH) = ©) + [OVE tg (pdt + 
0 
if fé 
=e [OO (2, T,T2, Q(T), G (T2)} dT AT. + + 
00 


oe foe plea GTO (ts) Gea) Ce, Ovo. dite) (5.10) 


where ®") is a symmetric isotropic matrix polynomial in g(t), g(T2), ---, J (T;); 
of degree k in the elements of the matrices and linear in each of them, and a 
BICLIOUOL 7. Tau Con, neck bee 


6. Hereditary Tensor Functionals 


Returning to equation (4.3) in which the functionals /* and /,, are given by 
(4.6) and (4.7), we see that if we employ the notation 


€pq(T) = 8p q(t} = Op7 (6.1) 
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the functionals /* and /,, may now be expressed by 


t fee a, t 
ie [@54(t)] = ye ae " Mp, a. sda ntrge UD ATS ee 7) x 
0 


R= 


X p,q, (Tr) pags (Ta) +++ Sovae (t,)d 7%, 0te0-- 0 T; 
and 


ie eo4( Soaps LN 4S P11 2292+ PKU (t, T, To; seey Tp) x 


es (74) Cpean CAR (%) €%,4%2...4T,, (6.2) 


where the M’s and N’s are continuous functions of their arguments. 

Let us now consider two deformations A and B of the body in which a point 
at X,;in the undeformed state moves to ¥;(t) and x,(t) respectively at time T, 
where 


x, (t + h) = %;(t), (6.3) 


h being a positive constant. Let g,,(t) be the value of g,,(t) and e,,(z) the 
value of e,,(t) for the deformation B. Then, from (6.3) 


Spg(t +h) =g,,(t) and e,,(r +h) =e,,(t). (6.4) 


The stress o;,(f) at time ¢ corresponding to the deformation A is given by 
(4.3) in which /* and f,, are given by (6.2). The stress o;;(¢+h) at time ¢+h 
corresponding to the deformation B is given by 


t+h 


ols (t +h) =P epq(0)] 84; + 


Ox; (t+h) dat +h ake ith 
+ [Ae Sa — Set + M8) frolsoa(el], (65) 


with 
th tth t+h t+h 


R 
Pleo =2 J JS Mpa pete... Pet (¢ + h, Ty, Te, -.-, Te) X 


, 
X ba, (t,) Con a, (T.) ... Cb, og (t,) dt,dT,... dT, 
and 


pes R tth t+h t+h 
/rs Bop (2)] ae i He Nrst.gits da tem ee Celestina Ee 


pea. (4) PACA ind Ge iT Ud Tagen te (6.6) 


Now, since in the deformation B, the material is undeformed prior to time h, 
we have ¢é, ,,(t)=0 for t<h, so that 


j* tah R t+h t+h t+h 
el = nie 
Bpal )] > / oh I Ue ees AU os h, Ty, To; Senu) x 


Crue) oye Ate oe (A dtd Toccoa te (6.7) 
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Employing the second of the relations (6.4) and the notation 7 =t —h, we obtain 


t+h 


(3 [gpa (7) 2 


t ¢ t 


neice ST Moat. peaeb th, +h, Th, ...,t+h) X 
Cnn (T) &p, gs (T2) sont Cp an (T,) aT, AT. eee aT, 


xe 
R tg t 
= sal JF Movatrta.traell+ hy AM, t+ hy ot +2) x 


k 


x ba, (7) Cb» qe (T2) ties ohn ak (T,) at, dT, 2509 At, x (6.8) 


Similarly, we can obtain 


t+h 


frs[8pal(™) 


t 
i Sie BE ie h, Ty h, T9 Wy. aby ory 
000 0 


7 
if TM» 


Saas Pe OS eae Coy ay, (Up) 4 Uy A Ty A Tp. (6.9) 


Introducing the relations (6.3), (6.4), (6.8) and (6.9) into (6.5), we see that if 
the M’s and N’s have the forms 
M p,9,..-Prav (tr Tas +> Te) = Mp.0,...tuge(t — T1, +++, — Tp) 
and 
Ny spy4,... uae (6 Ta +++) Te) = Nesp... teall— TM ---1t— Tt), — (6.10) 
then 


If the relations (6.10) are satisfied then the functionals / and /,, are said to be 
hereditary functionals. 

The stress o;; at time ¢ corresponding to the deformation A is then given 
by equation (4.3) in which 


t t ¢ t 


R 
fae ISpa(7) 2 ats SL Movnb.4---dame ll T,¢ Playas Lea Tp) x 


Kong ATs) ob a (Te) +> + Cpa, (Te) 4 Te... AT, 
and 


fle = 3, Sf al Neonat Jo-- wats 7,4 Patan 1 ep) 
Kaeo gal Ba pags (Talae® Cop agh tae tad Toi nid Ty (6.12) 


It is seen from equations (5.1), (5.2), (5.9) and (5.10) that if the material 


considered is isotropic in its undeformed state, the stress o;; is given by (4.3) 
in which P : 


R t 
File(s Niel ONT ee, Hog tine Rays 
0 k=0 0 0 0 


g (7); J (72), Se OAC) PLT AT, oi 0 Ep 


14 A. E. GREEN & R. S. RIVLIN: 


and t ee Mia ee 
Wts [8 pq(7)1 || ama > i frS BD” fF, Ty, To, +++) TR» 
0 k=0 0 O 0 

G(T). G (Te), G(T} AAT... aT, (6.13) 
where gp) is a homogeneous polynomial scalar invariant of the matrices g(t,), 
g(t), ---, g(t), linear in each of them, and a function 08 2; tH, tans Th and 
@") is a symmetric isotropic matrix polynomial in g (1%), g (T.), ---, g (t,), linear 
in each of these. The coefficients of the matrix products in ®” are polynomials 
in the invariants of the matrices g(z,), g(t2), ---, g(t.) and scalar functions of 
ty hin Cor ony ise 


With the notation (6.1) it is seen that (6.13) may be re-written as 


0 k=0 

Cg) eee) e(t,)}dt,dt,...dT; 
and t Beets ot 

Its [8p 9 (2) = ‘| ot wr ie, qT), Te; ey Tr 

@(c,)} @(6,) a (ee ee aa (6.14) 
where e(t) = ||e,,(z)||, y® is a homogeneous polynomial scalar invariant of the 
matrices e€(t,), €(t,),..., €(t,), linear in each of them, and a function of #, 7, 
T:,---,T, and W") is a symmetric matrix polynomial in e(t,), e(t2), ..., e(T,), 
linear in each of these. The coefficients of the matrix products in Y™ are poly- 
nomials in the invariants of the matrices e(t,), e(t5), ..., €(t,) and scalar functions 
GR arpa cant 


It is easily shown, in a manner similar to that employed in the more general 
case, when the material is not necessarily isotropic, that if the dependence of 


we and U''on fT, , te,..+, T, 1S Of the form 
p™ fi, Uy; Te picey Tay ClUy), © Ea) ec el uayy 
= yp ft iy Vi ncesiay bo Pag Te ClTs), 2 vO (tat 
and (k) § 
IE, Ty Caan 9 Uy CLT Ze eal ieeeen 


= Wt — t,t — ty, = ty, C(t), C(t), -.., €(t,)} (6-45) 


then for two deformations related by (6.3), the stresses are related by (6.11). 


7. Alternative Forms of Hereditary Functionals 


We consider the stress o;; at time ¢ to be given by equation (4.3) in which 
f* and /,, are functionals of the hereditary type given by (6.12). We assume that 
the kernels M, 4. 5,4,... puq, and Ny sid:bede-..Peqe 1 (6.12) become zero if any of 
the arguments ¢—t,,i—t,,...,t—t, is greater than some fixed time T. 
Physically, this means that the memory of the material considered is of limited 
duration 7; 7.e. the deformation at times earlier than T before the instant at 


which the stress is measured has no effect on the stress. We can then replace 
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the expressions (6.12) by 


t 


Ree) 
Pe al Os J lie fM, MIRON pene ek AAA? ceed 9 
0 p= = ge 


x “oa | (a) Chea ATs) oh Cras (Ui) 11d Coren OM, 


R if 
= i ites a Mp, 016, q2- aa (eae To; OOM) Tp) x 
RC. alle Tole Ce 10s). Cd Tost AT, 
and 
t R t t t 
alexaGall — if ate OS av Gio aaaee = Ue pe Cees nh Tp) x 
0 k=1¢—Tt—T t—T 


€p, 9, (T1) ep acs) PAO AT MATRA Up AT, 
Pre af 
dS a ’SPinP2%- Derar one To; OOO) T.) x 


x Cpag lt 1) Cpe, a To) ot Repay Pia th) OT ao. E,. (7-4) 


We now assume that e,,(t) has continuous derivatives up to order y+ 1 in 
the range ¢—T<1Xt, so that e,,({—t) has the Taylor expansion 


== Weave Sg 
epalt = 2) = DATE Set + Soul), 72) 
where z 
eT et peri GE 
Spqlt,2) ==" fw — ry Sel at, (7.3) 


Substituting from (7.2) in the expression (7.1) for /*, we obtain 


fUee(el] = J ff. iM Pi b2 Ge -: eee Cl Saree) x 


(—1) CT DB ep 
x | >, A a wie + Spyall tf dt dt. dt (7.4) 
jj ap=0 


If e,,(t) is a polynomial of degree vy in the range — T<r<t, then Sp, ap (t, T,) 
vanishes identically and (7.4) becomes 


t iQ 50 46 T 
ie [ep q(t)] = yy ig te EV GaG ree tras angels) 
0 jan OO 0 


+P. OCG eee XK eos ox D™%e D“e D*% 
Om ( rGa)) ( %) OO ! Tp) x An ete “et "bn a hay “ioe ATs, (7.5) 
Ay: A+... Hy: ID KE Di@ Di%* 


where, in the summation 2, under the integral sign, the «’s take the values 0 
to y. In a similar manner it can be shown that 


holEoa )] — ff. JN SPiN b29- Peon ae Wie Caress 0 LB) OX 


Oly Ca 
Jaa a: T 1)" (— a “ee a 14, )o% D "bd D eb. ds es D** "bn ak dt, dt,.. db, (7.6) 
G Oy log! ... a! Di Dir Di* 
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where, again, in the summation 2, under the integral sign, the «’s take the values 
0 to y. In (7.5) and (7.6) we may carry out the integrations with respect to 
T1, Ta) +++, T} and we then obtain 


ry R % CH an 
* = oe M3 (ame: DNs D "bs ds ae "beak 
/ [2p q(t)] ae a 7119 Po d2+++ Pk IK Di Dit Di% 
0 k=1 
and 
} R % Ap Ak 
= Dee NV (96 2 = x) D © beds D "bs de Pate D "bx q% (7.7) 
fis [8p (7), Nig hemr L**4S Pi G1 P2do++-PkIK Di2 Di® Dit , 
where 
fee SL ih 
%K (01 Oly». - &) (— ‘Wisk ee tele . ee ae 
iat — - eee 1 2 C 
Me ta ty! Oy! ... og! ff fa 2 --+ Up 
0 0 0 
X Mp, a Page -Peax (T Tp, +++, Tp) AT AT... dT, 
and 


d ira i Ti 
% (01 Oy.» - Of) (Sea eae ee BANG ta Oh 
182191 P2Jo+++PkIK Celtors lemon! 4 ee k 
00 0 


XN, sp:4:toqs.4. been eae tasenar phe ee ee ee oe (7.8) 


We note that /* and f,,, given by (7.7), and hence o;,; given by (4.3), are inde- 
pendent of ¢ and involve only the M*’s and N*’s, which are parameters defining 
the relevant physical properties of the material, and kinematically defined quan- 
tities measured at time t. 


In the case when e,,(t) is not a polynomial function of t, we assume that 
it has continuous derivatives up to order y+ 1 in the range ¢— T < t<t¢ and that 


D%* e;; (t) (a 
Pel <p (a =0,4,...9) 
and 
D+, ;(z) | 
sai |S00; ¢-T<t<p, (7.9) 


where P;# and Q;; may depend on X; but not on ¢ or t and where ¢ is a suf- 
ficiently small positive constant. It follows from (7.3) and (7.9) that 


y+. 
Sp, (be eee (7.10) 


If we expand the integrand in (7.4), we see that /* is given by the expression 
on the right-hand side of the first of equations (7.7), plus terms of the form 
e*B™ (% > 1), where |B@|<C™, a quantity independent of ¢. Similarly, it can 
be shown that /,, is given by the expression on the right-hand side of the second 
of equations (7.7), plus terms of the form «*B!® («= 1), where Bela Cea 
quantity independent of ¢. It follows that if €pq(t) satisfies the conditions (7.9), 
o;; may be expressed in the form (4.3), while /* and /,; are expressible approxi- 


mately in the forms given by equations (7.7) with the neglect only of terms of 
order ¢. 
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8. Alternative Forms for Isotropic Hereditary Functionals 


It has been seen in §6 that, if the material considered is isotropic in its 

_ undeformed state, /* and /,, are given by equations (6.14), in which y™ and 
Ww have the forms (6.15). We assume, as in the more general case discussed 
in § 7, that y and Y become zero if any of the arguments t—71,, t—T,..., 
t—vt, is greater than some fixed time T. It is easily seen, as in the previous 


/ case, that ; pe oe 


f* [pq (7) ] Sy, LE LOO a te te 


e(t—7,),e(¢—T,),...,e(¢ —%)} dt, dt... dT; 
and 


Ieslepa(t) | = 2 
e(, 7, ),€(¢ —7,), -0., C(t —T,)} 41,0 7.-.a%,. (8.1) 


Employing (7.2) and assuming that e,,(t) is a polynomial of degree y» in 
the range ‘—T<rt<t, so that S,,(t, 7) =0, it is readily seen that equations 
(8.1) may be re-written, with the notation e=e(t), De/Di=De(t)/Dt, etc., 


ole 

iol Epa aS fone MVD setae ete ec; Det ts Drei dud ty.08 Tp 

R110" 20, 0 
and 

t i Mage ciaieie 3 

Tye kpg®) Se alsa ye (*) 
0 =0 0 0 0 
KT Tope te, DO Dies. SDE DPEA Ts A Va.2H0 ts 21 (8.2) 

where y*™ is a homogeneous polynomial scalar invariant of the matrices e, 


De/Dt,...,D’e/D?’ and a function of 7, t,..., 7, and W*") is a symmetric 
isotropic matrix polynomial in these matrices, homogeneous in their elements. 
The coefficients of the matrix products in W*™) are polynomial scalar invariants 
in the matrices e, De/Dt,..., D’e/D?” with coefficients which are functions in 
T,, To, -.., t- Carrying out the integrations with respect to 7,, T;, -.., % in (8.2), 
we see that /* may be expressed as a polynomial scalar invariant of the matrices e, 
De/Dt,..., D’e/Dt’ and ||f,,|| may be expressed as a symmetric isotropic 
matrix polynomial in these matrices. We note that neither /* nor ||7,,|| depend 
on t, except through the matrices e, De/Dt, ..., D’e/D?’. 

It has been shown by RIVLIN & ERICKSEN (1955) that if we assume that the 
stress components o;; at time ¢ in an isotropic material are expressible as 
polynomials in the displacement gradients 0x,/0X,, velocity gradients du Ox,, 
acceleration gradients dv)/Gx,,...,(y—1)™ acceleration gradients dv}/0x,, 
where vb ) noe ..., v”) denote the velocity, acceleration, ...,(v—1)™ acceleration 
at time 1, then o (= ||o;;||) must be expressible as a symmetric isotropic matrix 
polynomial in the matrices C, A,, A,,...,A,, defined by 


Ox, Ox, || 
C= [Grell = [axe a3 

20g du, || 
A,= [40 =|52 + 5] 
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ang 4‘) (1) Y4) (1) || 
\| Dz ou Vig | 
A a fl =) Pate Ato) Me ANS Ee : 
a+1 APe || I De cs om O%q qm OX» i (8 3) 


Now, we have 


Dépy _ 009 O%m , OUD Om 
ID OX» OX, 

_ OU) OX O%m ag OUp O%n O%m 

Gee, OX» ox OR, OX, OX» 
OXm OXy AW) (8 4) 
~ _- OUD = mn’ . 
OX» OX, 


q 


It follows by induction that 
D% epg _ O%m O%n g(a) (8.5) 


Di ~ tke 


where A), is defined by the recurrence formula (8.3). We assume that the 
relation (8.5) is valid for «=a (say). Differentiating once with respect to time, 
we readily see that it is valid for «=a +1. Since equation (8.4) states that 
(8.5) is valid for a =1, it is valid generally. 

It is easily seen from (8.5) that the trace of any product formed from the ma- 
trices e, De/Dt, ..., D’e/Dt’ may be expressed as the trace of a product formed 
from the matrices C, A,, A,,...,A,. Consequently, since /* is a polynomial 
scalar invariant of e, De/Di,..., D’e/Dt’, and can therefore be expressed as a 
polynomial in traces of products of these matrices, it can be expressed as a 
polynomial in traces of products formed from C, A,, A,,..., A,. 

Also, from (8.5), it is seen that any symmetric matrix polynomial in e, 
De/Dt,..., D’e/D?’ may be expressed in the form 


where Q,; is a symmetric matrix polynomial in the matrices C, A,, Ay, ..., A,. 
Since ||f,,|| is a symmetric matrix polynomial in e, De/Dt, ..., D’e/D?’, it may 
be expressed in the form (8.6). 

Thus, from (4.3), we see that if the material is isotropic and e,,(t) is a poly- 
nomial function of 7, of degree » the stress matrix may be expressed as a sym- 
metric isotropic matrix polynomial in C, A,, A,,..., A,, in agreement with the 
theory formulated by Riviin & ERICKSEN (1955). 

In the case when e,,(t) is not a polynomial function of t, we assume, as 
in §7, that it satisfies the conditions (7.9). We can then show, in a manner 
similar to that adopted in § 7, that the stress matrix may be approximated by 
a symmetric isotropic matrix polynomial in C, A,, Ay, ..., A,, with the neglect 
only of terms of order e. 


CE Gay BOE 


(8.6) 


| ? 


9. Transformation to Convected Coordinates 
We now consider that the position of a typical particle Py of the body in 
its undeformed state is defined in a general curvilinear coordinate system #. 
We further suppose that as the body is continuously deformed the curvilinear 
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system moves with the body so that the coordinates of the particle P, in the 
system # are unaltered. The metric tensors of the system % at time t are denoted 
by y;;(t) and y*’(r) and these are defined by* 


_ Ox" (t) Ox” (t) Ariat a lu Ooi 
Vii (7) Gt ey ve Y 1(z) On” (2) eox™ (z) (9.1) 
We employ also the notation 
Vi O=Vi, PO=", 
750) =Lj, 770) =P". (9.2) 


We may define at each instant of time t the covariant base vectors y;(z) at the 
particle Py along the coordinate curves of the # system. Thus, y;(t) has com- 
ponents 0@x/(r)/0%* in the rectangular Cartesian coordinate system x*. The 
contravariant base vectors y‘(t) at the particle P, and time 7 are then defined by 


WT (a): (9.3) 


From (9.3) and (9.1), we see that y’(r) has components 0#‘/0 x/(r) in the system 
x'. We shall employ the notation 


TAU Sees OUP ie 
Y¥0)=N, yi(o)=I". (9.4) 


The contravariant stress tensor at the point # at time ¢ in the coordinate 
system # is denoted by z‘’. The corresponding covariant and mixed stress ten- 
sors are m,; and a. a'’,,; and a} may, of course, be related by the metric 
tensors y;; and y"’. 

We shall now determine the expression for the contravariant stress tensor 
zi in terms of the metric tensors and base vectors of the coordinate system #. 
We have ; 

wn 28 09 


Introducing the-expression (4.3) for o?% into (9.5), we obtain 


wii = 1*[soq(0)| 7! +(28 Bor — ten) ileal]. 0.8) 
Now, : ; 
q(t) = Yun) Soe axe (0.7 
Hence, 
81s =Vmn oe a (9.8) 


* We may, of course, raise at will the indices on any symbol defined in the 
rectangular Cartesian coordinate system, without altering the significance of the sym- 
bol. In this section we shall do so in order to obtain the conventional balance of 
indices in tensor equations. 


oF 
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From (9.6) and (9.8), we obtain 
ati = PUBpe(o) + 
+ (65th — 44 yan) 95 2 flee: (0.9 


If the functionals /* and /,, are taken in the forms (4.6) and (4.7) respectively, 
we obtain immediately with (9.7) that 


R 
P* [epg tl = Def Son fete ie Te Tapa 
0 00 0 


x Vim, Ny (7) Vm Ne (T2) rea Ymg nt (T) aT aT. ee dt, 


and 

GU OE ; . 7 STF mm my Maha «EM t ) x 

axe axe trs[Boe =X JIS ‘ (EE a Eppes Ee 

0 k= 
X Moin, (Et) Vomit; (Vo) <= Veo Ep) O04 @ Gon 
where 
KM Male MEM — apm oom dom OM Wake OOM OO Kei b2 42° ++ Pe 
aXh: AXn EXP: OX2 OXPr OX 

and 


iC? MN IN, Mz Ny ++ * MEN 


— CO oo" oom ee Om: eee Com OOM L801 2 92°** PEGE (9.11) 
OoxXt OXS OXh: OXU OX: 0X2 OXbr OX 


We note that the quantities 00"/@ X’ occurring in the above expressions are the 
components of the contravariant base vectors I” of the curvilinear system 9% 
in the undeformed state of the body. 

If the material is isotropic so that the functionals /* and /,, in (9.9) take the 
forms (5.9) and (5.10), then, with (9.7) and the notation 


V*(7) = ||P" ye,(0)| and P= ||P" ], (9.12) 


we see that 
es ees PAE Uy Mgt es ¥*(t%), ¥*(te), +++, ¥*(t)} dt, dty--- dt, 


R t t 


t 

— i i k 

ae sdiece DEO Maa tse 
OO: 0 


gm ag 
ax? axe trs 


k=0 
OP GCA coe ailcs WAL Meee 2c). (9.13) 


where g*") is a polynomial scalar invariant of Yaya), V7 einear 
in each of them, and a function of t, t,, t),..., tT and O@* isa symmetric isotropic 


matrix polynomial in y*(t,), y*(t2),..., y*(t,), linear in each of them and a 
PUCCIO Ol weer rain T, 


ks 
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The isothermal Expansion of a Gas Cloud 
into a Non-uniform Atmosphere 


M. H. ROGERS 
Communicated by G. C. McVITTIE 


A study is made of the isothermal motion of a spherically symmetric shell of gas 
surrounding a uniformly expanding sphere. The undisturbed gas ahead of the 
expanding region is assumed to have a density distribution given by the law (0 On timce 
and the flow is examined for different values of «. It is found that in all cases the 
flow is necessarily headed by a shock wave, and that there is a minimum value of 
the shock strength, which depends on «, in order that this type of flow can exist. 
The other effect of the density variation ahead of the shock front is on the thickness 
of the expanding shell of gas; as « increases up to its limiting value of 3 the thickness 
of the expanding region increases indefinitely. Approximate solutions are also ob- 
tained for the flow in certain cases and these are in good agreement with the exact 
solutions, obtained by numerical integration of the equations governing the motion. 


Introduction 


An idealized problem describing the hydrodynamic effects due to an early 
Q-star in an interstellar gas cloud is considered in this paper. As an O-star is 
forming its ultra-violet radiation will ionize the hydrogen in its immediate 
neighborhood, and this is referred to as the HII region. This then expands 
into the surrounding unionized gas and various estimates [see for example, 
J. H. Oort (1955) ] of the velocity of expansion are of the order of 10 kms./sec. 
As this region expands it compresses some of the surrounding gas into a shell, 
while at sufficiently great distances the gas is undisturbed; a representation of 
the motion is shown in Fig. 1. The isothermal expansion of a spherically sym- 
metric gas cloud has been considered by SAVEDOFF & GREENE (1955); the 
model employed consisted of an expanding HII region surrounded by a dense 
HI region, and this was taken to be expanding into an undisturbed H I region 
and separated from it by a shock front. The velocity of expansion of the dense 
shell was about 14 kms./sec. and its density was of the order of two hundred 
times that of the undisturbed region. However, these authors considered the 
case when conditions in the undisturbed H I region were uniform; in the present 
paper a study is made of the effects on the motion of the dense HI region due 
to non-uniform conditions ahead of the shock front. In this investigation the 
H II region is replaced by a uniformly expanding sphere, across which there is 
no mass flow; when the HIT region is present, there may well be some mass 
flow across the ionization front separating the two regions, but this will not be 
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considered. When the density and pressure of a gas vary in some fashion ahead 
of an advancing shock front, the flow behind the shock front can be affected 
in several ways. In certain cases it has been found (ROGERS 1957) that cavitation 
can occur; while in a study of the propagation of weak shocks, WHITHAM (1953) 
has shown that the strength of a shock wave can increase for particular density 
distributions. In the present paper it will be seen that the effect of a decreasing 
density in the undisturbed gas is to impose minimum values on both the shock 
strength and the velocity of expansion of the shock front. The density is taken 
to vary as a power of the radial distance, and this makes it possible to obtain 
similarity solutions for the flow. This particular choice of density distribution 
is also convenient for considering the case when a shock 
Wave moves into a cloud of increasing density; such a 
situation could possibly arise following the collision of 
two interstellar gas clouds, although it should be borne wp 
in mind that the restriction of a similarity solution has 
been imposed. Ce ae 
One feature of the motion of interstellar gas clouds is 
that the flow, in general, is not adiabatic. Now a complete 
discussion of the dynamics of the HI region would have By aean se Oey ee os 
to include energy losses due to radiation, collisions and center is an expanding H Il 
: 4 region of radius of R;; this 
other effects so that determination of the flow becomes js surrounded by a shell of 
extremely difficult. In the present work it will be assumed BElaee Rees 
that the temperature of the neutral hydrogen region re- to the undisturbed gas cloud 
mains constant, so that the flow is isothermal. This is 
equivalent to assuming that the time taken by the gas to reach equilibrium with 
the radiation field is short compared with the length of time for which the motion 
persists, and a discussion of this point is given by SAVEDOFF & GREENE (loc. cit.). 


The Equations of Motion and the Boundary Conditions 


In view of the results obtained by Sir GEOFFREY TAYLOR (1946) it seems 
likely that the expanding HI region will always be headed by a shock wave, 
and indeed this will be seen to be the case. Accordingly, exact solutions of the 
equations of gas dynamics are sought in which the velocity and density of the 
gas are of the form 


(1) WC I(2) 
{2) @ = e9(R) A(z), 


where c is the velocity of sound in the gas, z=7/ct is a nondimensional similarity 
variable and 0,(R) is the density of the undisturbed gas at the boundary of the 
expanding HI region. At this boundary there are two possibilities: in the first 
place a shock wave may be present, which necessitates abrupt changes in the 
velocity and density of the gas, or, alternatively, the transition between the two 
regions may be gradual, so that the velocity vanishes and the density of the 
gas is continuous across the boundary. Since 7 is the radial distance and ¢ is 
the time, a particle moving in such a way that its value of z remains constant 
is necessarily moving with a constant velocity, and this fact will be used in 
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considering the motion of the uniformly expanding sphere. In view of the 
assumption of spherical symmetry, the equation of motion is 


au Ou. 1 ap. 
8) i a 
while the equation of continuity is 
Comm 00 du , 20% 
(4) ie) ay ! or “i y ee 


where p is the pressure of the gas. To complete the equations governing the 
motion of the gas an energy equation is necessary, and as explained in the 
introduction the flow is taken to be isothermal so that the third equation is 
simply 

(5) p= ott: 


It should be pointed out at this stage that this equation is inaccurate near the 
shock front and a complete treatment of the problem would have to take account 
of radiation losses. The motion of a fluid in a field of radiation has been con- 
sidered by L. H. THomas (1930) and RossELAND (1926), while the effect of 
radiation losses within the shock front itself has been investigated by R. G. SACHS 
(1946). However, an important difference between these three papers on the 
one hand and the assumption expressed in equation (5) on the other should be 
noted; when the radiation terms are omitted from the results obtained by each 
of the three authors mentioned above, the equations of classical aerodynamics 
are recovered. In the present case this is not so, since the flow examined here 
corresponds to the singular case y= 1. 

On substituting equations (1), (3) and (5) into equations (3) and (4) the 
following two ordinary differential equations are obtained: 


(6) (-af=—*, 
(7) (an +A +2240) =0, 
where 
(8) — _Vt_ deo 
Qo(R) aR? 
a prime denotes differentiation with respect to z, and 
; _ aR 
(9) Karnes 


is the velocity of the outer boundary of the expanding HI region. If this 
boundary is a shock wave, then certain conditions must be satisfied here and 
these will now be discussed. In the first place mass must be conserved, and this 
can be expressed by the equation 


(10) 0o¥ = 0,(u— m), 


where the subscripts 0 and 4 denote conditions immediately ahead of and behind 


the shock front respectively. Secondly, momentum is conserved across the shock 
front, and this can be written in the form 


(11) 01%, (u, — v) = py — py. 
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Finally, the usual conservation of energy condition will not be used but, in 
accordance with the assumption of this paper, it will be replaced by the condition 
that the flow remains isothermal so that 


01 05 


Equations (10), (11) and (12) can be solved for the velocity and density of the 
gas immediately behind the shock front so that 


(12) Pi Po 


ce 
(13) m =0(1 =< 
and 
(14) oa 


It is then possible to fit these boundary values to the flow, if the motion of 
the shock front is such that it is defined by a fixed value of z, which will be 
denoted by z,. It follows that 


(15) go RiCt 
and so from this relation and equation (9), 


(16) Vie 


s¢ 


Equation (16) shows that the velocity of expansion of the shock front is uniform 
and, further, that the value z, is simply the ratio of the velocity of the shock 
front to the velocity of sound in the gas, or in other words the Mach number 
of the flow. The function ® which occurs in equation (8) can now be found, 
since in the present model the density of the undisturbed gas is taken to vary 
as a power of the distance from the center so that 


(17) 0o(7) = 0,7 * 


where 9, and « are constants. Of course, such a density distribution is unstable 
unless there is some force present which is capable of maintaining this con- 
figuration. However the point of this paper is to investigate the effect on the 
flow of a variation in density ahead of the shock front. Even if this particular 
density law is valid only over a limited region of space and for a limited period 
of time, the results obtained can still be expected to give some indication of 
the possible effects on the flow, and contribute to a better understanding of 
the hydrodynamics involved. It should be noted at this stage that if the effect 
of gravitation is included, either assuming the gas to be self gravitating or assuming 
the gravitational field to be caused by a central mass, then similarity solutions 
do not exist and the resulting flow would have to be determined from the inte- 
gration of partial, rather than ordinary differential equations. 


By using equations (15), (16) and (17) it can be seen that 
(18) D=—«%. 


Different values of « can now be substituted into equations (6) and (7) and 
the resulting velocity and density distributions found by numerical integration. 
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If the outer boundary of the expanding shell is not a shock front, then the 
boundary conditions would be, by equations (1) and (2) 


(19) F(Z) = 0 
and 
(20) h(z,) = 1. 


At the inner boundary, it is assumed that the expanding H II region is 
replaced by a uniformly expanding sphere, and consequently there will be no 
mass flow across this surface. If the radius of this sphere at any time ¢ is given 
by the equation 


(21) R, = 0,4 
then 

R; =. Ui 
(22) saa a 


and is therefore a constant. Furthermore, the velocity of the gas at z=z; is, 
from equation (1), cf(z;) and since this is the same as the velocity of the ex- 
panding sphere 


(23) cf (%;) = %; 
or, by using equation (22) this is equivalent to 
(24) I (2;) = 25. 


The value of the density at the inner boundary is not known, and will have 
to be determined from integration of the equations governing the motion. 


Variation of Velocity Distribution With a 


When equations (6) and (7) are combined, it is possible to obtain a first 
order differential equation for the function f and this can be written as 


af 2(f/z) — 
25 — es = 
(25) ars eae 
In general it does not appear possible to integrate this equation in a closed 
form: however it is possible to discuss the general behavior of solution curves, 
and then particular solutions can easily be found by numerical processes. 
From equation (25) it can be seen that the velocity gradient df/dz is infinite 


along the lines f=z-+1 in the (z, f) plane; physically this means that a shock 


wave must occur in the gas. In addition is zero along the line f=ye. 
Lg 


This line will intersect the lines / =z -+ 1 at two points which are therefore singular 
points of the differential equation (25). Since we are only interested in positive 
values of z, the relevant singular point is found to be 

(26) = 


Be | 
A wars 


’ 


(27) {j= tS 


= oe 
where the positive sign is to be taken when « is less than 2 and the negative 
sign when « is greater than 2. The nature of the singularity can be determined 
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from the classical theory of non-linear differential equations (STOKER 1950); 
when « is greater than 8/3 the singularity is a spiral point, otherwise it is a node. 


4 
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Fig. 2a—d. Velocity fields for different density distributions in the undisturbed gas. In each diagram the line 1 is the line 
1 
{=z and represents the surface of the uniformly expanding sphere. The curve 2 is given by f=z— oF and represents the 


shock front; the region between the two represents the expanding shell and the solution curves of physical interest are 
i shown as unbroken lines. a) «=0 


ule 


a 7 zs 3 ¥ r/ct 5 
Fig.2b. O0<x%<2 


It is now possible to determine the velocity field for different values of « and the 
results are displayed in Fig. 2. In each case, a particular solution starts on 
the line /=z, which represents the surface of the uniformly expanding sphere, 
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and can be followed until either it crosses the z-axis, which means that the 
velocity of the gas vanishes and we have reached the edge of the disturbed region, 
or alternatively the velocity gradient becomes infinite. This latter condition 


u/c 


go 7 | 
a J ¥ 7r/cr J 
Fig.2d. 2<a<3 


means that a shock wave must exist in the flow, and from equation (13) it is 
a simple matter to determine where the shock front must occur. By using equa- 
tions (1) and (16) the shock wave boundary condition can be written as 

(28) [Ge oe 


&s 
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This curve has been plotted on each of the four cases in Fig. 2 and two results 
are immediately apparent. Firstly, no matter how small the velocity of ex- 
pansion, a shock wave must always occur: no solution curve starting from the 
line f= can cross the f-axis without first crossing the line f/=z—1. This result 
is in agreement with that obtained by Sir GEorrREY TAYLOR in the case of a 
uniform atmosphere. Secondly, in all cases except «=0 (which is the case of 
a uniform atmosphere ahead of the ¢— 
expanding shock front) there is a lower 
limit to the shock strength in order 
that this type of flow may exist. The 
solution curve corresponding to this 
critical value passes through the origin 
and meets the shock curve at the point A. 
If the shock strength is weaker than 
this critical value, then it can be seen 
from Fig. 2 that a solution in which 
the flow is continuous all the way from 
the shock front to the surface of the ex- 
panding sphere does not exist: the fact 
that for very weak shocks the solution 
curve crosses the z-axis and gives nega- 
tive values of / would seem to indicate | 
the possibility of a back flow in this 
region. An approximate form for the re | 
critical solution can be obtained using ? 7 2 7 3 


— 


ae 


! 
| 
1 
| 
| 
| 
| 
4 
! 
| 
! 
| 
! 
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a power series development for /, and 
by calculating the coefficients by direct 
substitution into equation (25); the 
leading terms are found to be 


Poe = 24 =($ 1) 24 


Fig. 3. Curve showing critical value of Mach number plotted 
against «; the case ~=3 corresponds to an infinite value 
for the mass of gas in the initial configuration. Note that 
when the density decrease is less than the inverse square 
law, the critical Mach number is quite small; large values 
are attained only as % approaches the limiting value of 3 


ae ee \ 5% ig gale elt. 
tale VS 3] ies 


Now at the shock front z=z,, and from equations (28) and (29) it follows that 
z, and « are related by the equation 
is 
3 


Ay 28 oe (< 
: LOSES 
The terms in the series on the right hand side decrease quite rapidly provided 
that z is not too large, and in the case «=1 a value of z,=1.271 was obtained 
by retaining the first two terms only, as against a value of 1.2844, obtained by 
numerical integration. However for larger values of z, numerical methods were 
used and the results are given in Fig. 3. This figure shows the variation of « 
with z,, the critical Mach number for the flow. It can be seen that as « ap- 
proaches the value 3 the Mach number increases rapidly, and in fact for the 
limiting case «= 3 (which corresponds to an infinite mass of gas in the original 
configuration) the solution reduces to the blast wave solution, discussed in a 
previous paper (ROGERS 1957). 


(30) Bs 5 = 7% = (5 oes “e 


“9 
2) 
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Density Distribution 


The density of the gas is, from equations (2) and (17), 


(31) e =e. R*h(z). 


The function h(z) can be found from equation (6) by numerical integration 
once f(z) has been computed from equation (25). Numerical integrations were 
performed, in the case «=1, for different values of the shock strength and the 
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Fig. 4. Density Profiles in motion surrounding uniformly expanding sphere. The single curve meeting all the profiles 
represents the shock front, hs=2z3. The value «1 has been used here 


results are shown in Fig. 4. These calculations and all the other numerical 
integrations in this paper were done on the University of Illinois digital computer. 
It is interesting to note that although the value of the density at the surface 
of the expanding sphere is not known it is still possible to start the integration 
at this inner boundary with an arbitrary value of # and then correct it at the 
shock front by using this boundary condition, equation (14). The reason for this 


Table 1 

M eee | aR ee. | mile w(2)/23 
1.2844 1.319 [o.e) .799 .O -500 
1.30 1.589 1.482 940 | 524 500 
se S@) 2.399 413 1.066 1.061 500 
2.0 4.218 -162 1.055 Wefan .500 
5.0 25.247 021 1.010 2 899 500 
10.0 100.249 .005 1.002 ).950 -500 
15.0 225250) .002 jtcOOd 4.967 500 


is, of course, that the equations governing the r otion are homogeneous in 0, 
and hence in h, due to the fact that the flow is isechermal. To obtain a physical 
picture of the expanding region it is convenien. to compute the ratio of the 
density at the inner boundary, that is at the surface of the expanding sphere, 
to its value immediately ahead of the shock front and these results are displayed 
in Table 1, Column 2. It can be seen that as the Mach number of the flow 
increases, the density ratio increases roughly as the square of the Mach number. 
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Column 4 in this table gives the ratio of the density at the inner boundary to its 
value immediately behind the shock front; it is seen that this ratio is close to 
unity in the cases except near the critical solution, shown in the first row. 

The total mass contained in the expanding region can be calculated quite 
simply: because of the similarity assumption it is evident that this must be 
equal to the total mass of gas originally contained in a sphere of the same 
radius R. Hence from equations (2) and (17) if follows that 


R R 
(32) 4a fo,r*rdr=4n fordr 
0 Ri 


and by using equations (15) and (16) this is equivalent to 


as 
3 
(33) sop = [PAW az =o). 
2G 

The integral w(z,) was evaluated numerically and the results are shown in 
column 6 of Table 1; this serves as a useful check on the accuracy of the calcu- 
lations. 

The density distribution in the case of the critical solution has been found 
by numerical methods, but it is worth nothing that a series expansion is also 


Fig. 5. Density Profiles at successive times, in the case ~=1 


possible. When the expression for /, equation (29), is substituted into the equation 
of motion (equation 6), the following series expansion for / is obtained, in the 
case «= 1 


(34) h=h[1tse?+ eet O()]. 


In this expression , is a constant of integration which can be determined from 
the boundary condition on the density at the shock front, given by equation (14). 
Succeeding terms in this expansion decrease quite rapidly, so that equation (34) 
gives a good approximation to the density distribution in the critical case. Some 
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idea of the accuracy of the expansion may be gained by comparing the approxi- 
mate value 1.337 of f,, obtained by including the first three terms only of 
equation (34), with the value 1.319 obtained by numerical integration; an error 
of less than 2%. The actual value of the density at any time is of course given 
by equations (2) and (17) 

(635) 0 = 0,R-*h(2) 


and in order to get a clearer picture of the nature of the flow, it is of interest 
to evaluate this at different times during the motion. Accordingly, in Fig. 5, 
the density profiles at successive intervals of time are shown for the critical 
solution, and this demonstrates the effectiveness of the shock wave in producing 
a more uniform density distribution within the gas cloud. 


Approximate Solution for Strong Shocks 


As can be seen from Fig. 2, the thickness d of the expanding region is small 
compared with the radius R; of the expanding sphere for strong shocks, provided 
that « does not assume values close to the limiting case « = 3. Column 3 of Table 1 
gives computed values of the ratio of d to R; for the case « =1, and it is seen that 
even for a Mach number of 2 this ratio is quite small. Since the chief interest 
lies in flows headed by shock waves having a Mach number of the order of ten, 
it is worth while to investigate an approximate solution of equations (6) and (7), 
which is applicable to such cases. 


A new coordinate # is therefore defined by the relation 
(36) 2=2,(1 +8) 
so that at the shock front, where # assumes its maximum value @,, it follows that 


a! eG 
(37) a, pee nad Te 


If series expansions for f and h are now substituted into equations (6) and (7) 
and the boundary conditions (24) employed, the following expressions for the 
velocity and density are obtained: 


(38) he + (x —2)%+ 3 —a) +... 
ads. (a — 3) (a —2 a — 3) (8 —3a 
(39) h=h1—4 Aa — 2) 5292 Bas poe 4 | 


Oo 


where h, is a constant of integration which can be determined from the boundary 
conditions at the shock front. From equations (13), (14) and (16) we have 


(40) ie == Bs — aie 
and - 
(41) hows 


When these results are substituted into equations (38) and (39) it is found that 


; 1 
a) Gia so Oe 
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and 


(43) h=#|14+ °F* 9,4 18+ 009). 
2) 


Equation (42) shows that this approximate solution will be valid as long as 
(3 —a)~1z;* is small: this method is applicable to the motion behind strong 
shocks in all cases except when («—3) is small. Equation (39) shows that the 
density distribution is parabolic to a first approximation, except for the case 
%—=2 when a cubic profile results. Table 2 gives computed values in the case 
% = 1 for varying shocks strengths: on comparing these with the results in Table 1, 
which of course were obtained by numerical integration of the exact equations, 
it is seen that the accuracy of this method is good down to shocks of Mach number 
two. Column 5 in both Tables1 and 2 gives the ratio of the velocity of the 


Table 2 
M 0/20 d/Ri=9s5 0i/01 u;/c 
—__—— SS 
2D 4.19 rie ALO ihe Fh | 
3 92500" ox. 068 1.30 2.82 
5 25.25 -021 1.01 4.90 
10 100.25 005 1.003 9.95 
15 0a 20525 Bl O02 Bi 0Gt_ | 14.07 


surface of the expanding sphere to the velocity of sound in the gas, and for shocks 
with high Mach numbers this is very nearly the same as the Mach number of 
the shock itself. 


Discussion of Results 


One of the chief effects of the variation in density of the undisturbed gas 
ahead of the shock front is on the thickness of the compressed expanding region. 
From equation (37) it is seen that the ratio of this thickness to the ratio of the 
uniformly expanding sphere is simply given by #,, and the approximate value 
of this is given by equation (42) in cases when this is small. However, when « 
assumes values close to the limiting value of 3, this approximation is no longer 
valid and it can be seen from Fig. 2 that the thickness of the shell increases 
considerably in such cases. This indicates that when an expanding wave of the 
type considered moves into even a limited region in which there is a sudden 
decrease in density, considerable expansion of the compressed region will occur. 

The existence of the minimum shock strength is also interesting since the 
boundary of the flow, the uniformly expanding sphere, does not exist in this 
critical solution. The corresponding density profile shows the effectiveness of 
even weak shock waves in smearing out non-uniform density distributions in a 
gas cloud, which could for example be caused by turbulence or by collisions of 
interstellar gas clouds. Another general feature of the solutions presented here 
is the accuracy of approximate solutions under certain restricted conditions; 
the complete equations governing the motion of a compressible gas present a 
serious obstacle in the study of cosmical gas dynamics, but the present results 
demonstrate that in certain cases a relatively simple approximate solution can 
give excellent agreement with the exact solution. 
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Finally, it may be remarked that by fitting non-homogeneous boundary 
conditions to progressive wave type solutions of the equations of gas dynamics 
it is possible to give a physical interpretation to a large class of solutions, and 
this method is capable of extension to other flows of compressible fluids. 
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The Pi Theorem of Dimensional Analysis 


LovuIs BRAND 


Communicated by C. TRUESDELL 


1. Introduction 


The fundamental principle in dimensional analysis is known as the Pi Theorem. 
While the ideas involved were used by earlier authors, BUCKINGHAM [/] stated 
the theorem essentially as follows: 

If an equation in n arguments is dimensionally homogeneous with respect to 
m fundamental units, it can be expressed as a relation between n —m independent 
dimensionless arguments. 

It is well known that this useful practical rule is not strictly correct [2]. We 
shall give a simple, constructive proof of the theorem in its rigorous form. This 
makes no use of partial differentiation [3] or abstract spaces [4, 5] and is based 
upon a simple idea already used in dealing with homogeneous functions. Although 
dimensional quantities may be regarded as elements of a vector space defined 
by a set of postulates [6], the Pi Theorem needs no elaborate logical setting. 
The proof given in § 4 depends only on matrix algebra. However we shall sketch 
two other recent treatments, due to BiIrRKHOFF [7] and Drosor [8] respectively, 
at the end of this paper. 

In applying the Pi Theorem to actual problems, considerations beyond the 
domain of pure mathematics may enter. Thus certain “paradoxes” have arisen 
which have occasioned much discussion [9]. These paradoxes are not due to 
any failure of the Pi Theorem; for this is a straightforward mathematical pro- 
position of universal validity. The successful application of dimensional analysis 
nearly always depends on a real understanding of the essential variables involved 
in the problem. But for every physical assumption, the Pi Theorem gives a 
corresponding answer; and experiment alone can decide which of several answers 
most nearly matches the facts. We proceed therefore, without physics or meta- 
physics [9], to lay down a few definitions essential for a precise statement of the 
Pi Theorem. 


2. Dimensional Matrix 


We consider problems in which the physical quantities X; involved have 
positive measures x; which depend upon a system of m fundamental units 
U,, U;,..., U,,. When these units are changed to 


m* 


(1) G=Ul; (>), 
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the positive variables x; also change. If the new value x; is related to the old 


by the equation 


/ 


(2) of = I ge. Him x, 
we say that x; (strictly X;) has the dimensions 


(4:1, Ajay ss Lim) 


in the units U,, U;,..., U,-. If all a;;=0, %; is said to be dimensionless. In any 


case the dimensions a,,; are real numbers. 
The dimensions of 2 quantities X; may be arranged in a rectangular 1 xm 


matrix: 


| U YU, U, 
% | 41 Ne Am 
| 
xe, | de, A229 Gm 
ele aoa ie 
Xn ay, ab ane shee ay m 


We denote this dimensional matrix (a;;) by A. 

In mechanics the fundamental units are of length L, time 7, and mass M, 
and m=3. For example in dealing with the speed of sound v in a gas, the addi- 
tional variables involved are the pressure #, density @, and viscosity of the gas; 
and the dimensional matrix has the form: 


| yah) Raga A 
ae ieee a 
bh aAipsaitont 
0|\—3 @. .4 
wi\—-1—-1 1 
Lemma. I/f %,, %3,..., %, have the dimensional matrix A and the products 
(3) Cae AL CREE BaD Gatto 26D) 
have the exponential matrix (b;;)=B, then y,, Y9,..., Vp, have the dimensional 
matnx BA. 
Proof. From (2) we have 
A — ba 5°42 he x7 Pin = Ue Le aur on Vi 


where C;;=6;,4;+---+6,;,4,;. Thus the dimensional matrix (c;;) of the y’s 
is C ==.BAy 

A matrix is said to be of rank r if it contains at least one non-zero determinant 
of order 7, while all determinants of higher order which the matrix may contain 
are zero. If the matrix contains no determinants of order r+, it will be of 
rank 7 if it contains one non-zero y-rowed determinant. Thus if m<n, thenxm 
dimensional matrix A is of rank m when it contains a single non-zero m-rowed 
determinant. This is the usual situation in dimensional analysis. 


The Pi Theorem 47 


3. Isobaric Functions 


Let the measure x of the physical quantity be expressed as a function of 
SCRE eet ete 


aU eee ay, |: 


If this equation holds for all changes of units, it becomes 


x! = f (1, Xa, --+5 Xp) 
when U; = U,/t;. Hence if x has the dimensions (a,, a, ..., G,), We have from (2) 
tr tts... ttm x — f (Mute... tmx, ...). 
Thus under the given unit change the function satisfies the identity in ¢,, fg, ..-, by! 
(4) TORE UCU AROS aR GAIN bie 2M Pere 5 rare aid 2, oe 
Such a function is said to be zsobaric with the dimensions (a,, a), ..., @,,). We 


use the term “‘isobaric’’ (Greek zsos, equal; barus, heavy) in preference to the 
more cumbersome phrase “‘dimensionally homogeneous’. 


We note that isobaric functions relative to m fundamental units U,,..., U,, 
are also isobaric with respect to any subset of these units, say U, U,,..., U, 
(k<m); for in the identity (4) we need only put 4,.,=4..=-:-=4,=1. Thus 
the isobaric functions relative to U,, U,,...,U, form a larger class than the 
functions isobaric relative to U, U;,...,U, (m>khk). In brief, the fewer the 
fundamental units, the larger the class of corresponding isobaric functions. We 
shall later have occasion to use this fact in explaining an alleged paradox. 


If the sum 


(5) fH=O% CoX- Rate Cy Xn 


with dimensionless coefficients 1s isobaric, each of its terms must have the same 
dimensions. For if we apply (4) and put ¢;=7¢ and all other ?’s equal to 1, we get 


O17 Cy xy 5 WIC» Noni ein 18 Cy Lh a i% f 
or 
£011-— D6, 4 OI Co Hg ts HIGE 6, x, =f. 


Since the right-hand side of this identity in ¢ is independent of ¢, the same must 
be true of the left-hand side; that is 


yj = Oyj = 87 = 1,2,...,m). 


All rows of the dimensional matrix are now identical. 
The product 
k. 


(6) eee ae oe 


is an isobaric function whose dimension a; in any unit U; 1s obtained by multiplying 
the dimensions of its factors by their exponents and adding the results: 


(7) hy Ay; + Re Gaj + 0+ + Ry Ong = 4;- 
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For if we change only U,, say U =U//t, we have 
(3 5)" (7 Xo)" Pie” see ph14ij + ea; free Vip 


Since this is true for all £>0, we see that f is isobaric with the dimension a, 
given by (7) with respect to the unit in question. 


4. The Pi Theorem 


Let the function f in an equation 


(8) A ae Oe =O 

with n arguments be isobaric with respect to m fundamental units U,, Up, ..., U,- 
Then tf the nxm dimensional matrix of %1,,..., %, 18 of rank r=n—k, the given 
equation 1s equivalent to 

(9) CO mene eh te chen) 


im which the first y arguments are 1, and the x's are n —r independent and dimen- 
sionless products formed from x1, ..., X%p- 

Proof. Let us first consider a special case. In ordinary analysis a function 
is said to be homogeneous of degree a, if for any {>0, 


(10) bX Kn) 05 tq) =F) tere aa ee 
is an identity in ¢. From our present point of view / is an isobaric function of 
dimension a with respect to a single unit U,; and its arguments i eae eee 


have the dimensional matrix consisting of a column of 2 ones. Since (10) is an 
identity in ¢ it will hold when t=1/x,; then we have 


z a = 
i sa ise | ae PE (ei te 453, ele 


Ma Ap 
Hence the equation f/(x,, %,..., x,)=0 is equivalent to the equation 
f(t, et 2) = 0 
ie Hy 
in the 2—1 dimensionless products 2, = %q[ 1, ..., Uy, =%,_/%,. Since the 


dimensional matrix is of rank 1, we have proved the Pi Theorem for this case. 
The proof in the general case consists in a natural extension of this reasoning. 

By suitably numbering the arguments x; and the units U; we can bring a 
non-singular 7x7 matrix P into the upper left corner of the dimensional matrix 
A =(a;,;). Then A may be written as a matrix of matrices 


1 i 
(11) A=(é or where det P+ 0. 


If n=r-+hk, the sub-matrices POURS S rate respectively 7x7, kxr, rxk, kxk. 
We first show that the sub-matrix 
(12) Die ORR 


Since A is of rank 7, the last k rows of A, namely (Q, S), are linear combinations 
of the first y rows (P, R). Thus there exists a Xr matrix C such that Q=CP, 
S=CR; hence C=QP4 and $= OP AR: 

After this preliminary, the proof in the general case is based on the same 
idea as that used in the preceding special case. In the defining equation (4) 
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of an isobaric function we set the first 7 arguments equal to 1; if we put 


4,4,=+::=1t,,=1 in these equations we can then determine ¢,,..., t, uniquely 
since det P=+:0. These ¢-values then define a change of units which converts 
Sy Megas 3s Sp, X, 475-2, X)\intoa multiple OF LUsAS4.. $4; 245. .2, 2;). .Lloweyer 


we present the argument in a modified form which is easier to follow. 
If the yxr matrix P4=(b,,), consider the r quantities 


(13) ie Aa Ane ail 44,2). 7. 


Since %,, %),..., x, have the yxm dimensional matrix (P, R), the lemma shows 
that the quantities y; have the yx m dimensional matrix 


PAP Rea PER), 
where J, is the unit 7x7 matrix. Now take 
= — = — ewe —- t —— es) t = 
Vy ) ty Vo ? , t, Vy ed don ib m 1 

in the defining equation (4) of an isobaric function. Equation (8) is then equi- 
valent to 

Xx. x x. Xy4+ x. 
14 . (=, EES Sil sid oe er ety 
(14) i Bay? is Fees Orta) | 
where the » quantities 
(15) ee VET se NaCe, ee rar Wi: 

nee be 3 , 

Since their x7 exponential matrix 3( } the lemma shows that the z’s have 
the dimensional matrix 


P Athy if Pate oa ABUR 
(g)%? Bale foals 5 


in view of (12). Thus %,..., 2, have the same dimensional matrix as x,, ... 
consequently all the arguments x;|z; in (14) are dimensionless. 


? Xn> 


Finally, if we substitute y,; from (13) into (145) we have 


(16) Fe ARO I ge Veh Ds parece 


Av fie 
To find the xv exponential matrix we have from the lemma 


k 


(cJen=(5 - ) 
C; Q QP‘R 


poses Cpe Pe Ges ue Oe OC = OP 
The z’s in (16) therefore have the exponential matrix 
l, 


that is, z;= x; (¢+=1, 2,...,7) whereas 2,,,,..., %, have the exponential matrix 
Q P+, Equation (14) now becomes 


(18) #14, 24, 2, ..., 22) =0, 


PT 
eyAy En 
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an equation with n —r dimensionless arguments, which are obviously independent 
since each contains an x, (in the numerator) that none of the others contain. 


On writing 
x 


4 Att x, Nai PA ecg dd 

( 9) Byta 

in (18) we obtain the equation (9) required by the theorem. These dimensionless 
arguments are readily found since z,,,,...,2,1, are given by (16) with the 


exponential matrix QP. 
The above proof now enables us to state the Pi Theorem in a still more specific 
form. 


Let x1, %,,..., %, have the nxm dimensional matrix of rank r=n—k: 
et gc 
(20) ans est 
(CRO YB § 
where P is a non-singular rxr matrix. Then if f(x,, X,..., X,) is an isobaric 


function with respect to m fundamental units, the equation (8) is equivalent to equation 
(9) in which 
Th he Ne at, toe eyes 


are k=n—r independent and dimensionless quantities with the kxn exponential 
matrix 
(21) LOL ae 
where I, is the kX k unit matrix. 
In brief, the Pi Theorem states that 
(22) EA=0. 


As an example, consider the 6x4 dimensional matrix of rank 3: 


On ae Una et), 

Gh teunO At tOa el 

Xo | O-—4 yl 2 
eee 1 | » 
an 1 pes 2) 

Xs 1 2a 

Veal eoeel 1 2 4 


To verify that A is rank 3 we note that the vectors [4], [%5], [%,] in the last 


three rows are linear combinations of [#1], [%el, [4¥9)3 for example 


[%4] = —$ [4] —2 [x2] — 4 [oy). 


We may take the 33 matrix in the upper left corner as P since det P=6. 
Then Q is the matrix in the lower left corner and 


= sae ara Ye We ads 2 ot, pees 

OPT = 162 ao? | eee me gales Ope ee 
SNL pe 17) ah ae 

gall Renna S: 
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The dimensionless products have the exponential matrix 


PAO d 10% 
OP yet Ol Wg i Ole Oe] 
3-$-$ 001 
hence ay Gre Mike 4 8 yb yo} 
Uieeaed Ue 2 BAe AOC 8 ee oe MENG Eh gy Ver ge %G 


If we wish to avoid fractional exponents, z%, as i, may be taken as the di- 
mensionless products. 

We may readily verify that EA =0, thus checking the calculation. The fact 
that S=QP1R, also shows that A is of rank 3. 


5. Fundamental Units 


In a certain area of knowledge let the quantities X involved belong to a 
certain class @ which includes the real numbers. The quantities U,, U,,..., U,, 
is said to form a set of fundamental units in the class @ if every X €@ may be 
expressed uniquely in the form 


(23) X = xU%Us...U%m, x > 0; 


here x is the measure of X and aj, a,,...,a,, are real numbers which give the 
dimensions of X in the units U,, U,,..., U,, respectively. The dimensions of X 
form a vector which in MAXWELL’s notation is denoted by 


(24) eC he AGP baer UL 


It is customary, however, to write x for X in (24); then [] denotes the dimensions 
of the physical variable X whose positive measure is x. We shall adopt this 
convention. 

When the units are fundamental X = x, a pure number, when and only when 
X is dimensionless, that is, when [X]=0, the zero vector. Thus when the units 
are fundamental, 


(25) [xa 0. implies .X6== x; <a. positivernuinber. 
If U,, U,,..., U,, form a system of fundamental units, the system 
Ve Up Ue 2 U beer, (=A On 


is also fundamental provided det 0;;=- 0. For the system of linear, homogeneous 
equations which state that 


X = x VnVoe,.. Vom 


is dimensionless in the units U; has the solution «, =a,=-::=«,,=0 when and 
only when det 6;;=-0; then X =~, a positive number. 

In mechanics units of length, time and mass (L, T, M) are usually chosen 
as fundamental. But engineers often prefer to regard units of length, time and 
force (L, T, F) as fundamental. This is permissible since #=LT~?M and the 
determinant 


Pie Or SO 
foe ti 2OieeO. 
Pee 
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To illustrate the importance of the choice of fundamental units let us consider 
the following problem [7] which leads to an alleged paradox in dimensional 
analysis. 

A ball of diameter d is fixed in a stream of liquid and kept at a temperature 
# above that of the liquid at a great distance from the ball. If the velocity of 
the stream is v and the liquid has the heat capacity ¢ (per unit volume) and 
thermal conductivity k, find the rate y at which heat is transferred from the ball 
to the liquid. 


Solution 1. Units of length, time, mass and temperature (L, T, M, H) taken 
as fundamental. 

The dimensions of calories are those of energy, L?7~?M. Since c is given in 
calories per degree per cubic centimeter, 


icv (2 772M) AL = bey i, 


The thermal conductivity k is given as calories per second per square centimeter 
per unit of temperature gradient; hence 


VEST EOE 


aa 


The rate v of heat transfer is given in calories per second; hence 
(7) Sears i ae 


Thus we have the dimensional matrix: 


| £ i Viaawel 
alt On0r ea 
a) a (Ee S= 40) 

ote 

k pri.arg £24 
Gir) OU" o-Or Von 1 
Re See re a) 2 
r soe Yh mene 


This matrix has the rank 4; for the 4x4 matrix P has det P=—1. Hence 
there are 6 — 4=2 dimensionless products. Since 


{ ORORO 

Pi ee ae! x OP = olbrear tan tea. 
Di = 3 ed 4, One. Waele 
NOD OR OR 


B= (OP 1) =( ih kes aaron an 
Aly 0, —1, —1, O, 4| : 


* To find P4 partition P into four 2x2 matrices & a 


have P= = ( A 0 Oe 
=D GAs 


i then since B=0O we 
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they are 
Ty =AVRAC, To =a k19 47, 


Thus a relation between the six variables may be put in the form 2,—f(z,), or 
(i) r= f(a,)ake: 


Solution 2. Units of length, time and mass (L, T, M) taken as fundamental. 


If the temperature of a body is regarded as a measure of the kinetic energy 
of its molecules, we have H—L?T-2M. Our dimensional matrix is now 


hg oN 
d er a) 
ese tt FON 
0 ee eR 
eS ae Ore 
v Ae Amal cere 
Y 2 A 


The rank is 3 since the 3x3 matrix P has det P=—1. The Pi Theorem now 
gives 6 — 3 =3 dimensionless products. Since 


| 
pe — 
| 
iS) © 
= j=) 
O 
a 
ny 
| 

no YH 
—_ —- © 
iS) eS 


RE ORM mt! 
(OP eee = 40, 08 taa0 
Ef 1 4 026-4 


they are 
TUNA. Oe ea That Soe Un oe Use he ete 


A relation between the six variables may be put in the form 2;=g(z,, 2) or 
(ii) Pap (wie) a hee. 


Since functions of two variables form a much wider class than functions of 
one variable, the added knowledge that H=L?T~?M given by the molecular 
theory of heat seems, as Lord RayLeicu said, to “put us in a worse position 
than before’’ [12] when the desired relation was given by (i). The explanation 
of this apparent paradox is given in § 3; for solution 2 with three units admits 
a larger class of isobaric functions than solution 1 with four units. 

The paradox may be explained also as follows [13]. The assumption that 
L, T, M, H are fundamental units implies that 


Bice co Bek Bh MG 9 he 


is dimensionless when and only when 


(iii) (ape = a0! 
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In solution 2 the assumption that L, T, M are fundamental units implies that 
X =x L*T’ MM (L?T?M)* 

is dimensionless when and only when 

(iv) == 2d, 20=20, lo 2. 


Now conditions (iv) are included in (iii) when d=0, and are therefore less 
restrictive than conditions (iii). Thus the Pi Theorem gives formula (i) under 
one theory, formula (ii) under another; only experiment can decide which is 
correct. 

6. Recent Developments 


In 1914 BucKincHam [1] proved the Pi Theorem for functions capable of 
being expanded in Maclaurin series. This restriction is not imposed in BRIDGMAN’S 
proof, but the functions are assumed to be differentiable and the proof depends 
on the solution of a linear, partial differential equation of the first order [3]. 
The problem, however, is strictly algebraic and differentiability is not a relevant 
requirement. Purely algebraic proofs have been given by LANGHAAR [4] and 
BrirKHoFF [7]; but neither author gives the theorem in the specific form stated 
in § 4, where the precise function of the dimensionless arguments is given as 
well as their exponential matrix. BrrRKHOFF’s argument is difficult to follow in 
view of its extreme concision and a variety of misprints; moreover his definition 
of “unit-free’’ functions f as those for which “‘the locus defined by f= 0 is in- 
variant under all transformations [of units]’’ is not well adapted to a sharp proof. 

S. Drogot in his recent paper On the Foundations of Dimensional Analysis [ 5} 
aims “‘to construct the Dimensional Analysis by means of quite simple algebraic 
methods, namely using the theory of linear space’. In part II of his paper the 
usual postulates for a linear (or vector) space 2 over the field of real numbers 
are given. The operations involved are the addition of elements of »' and their 
multiplication by numbers. He then proves two theorems on the form of func- 
tions whose arguments, as well as the function itself, are elements of X. In part 
III this entire theory is carried over bodily to a ‘“‘multiplicative form of linear 
space’’ [I whose postulates are precise analogues of those for X. The elements 
of /T are called “‘dimensional quantities’, and the positive numbers are regarded 
as a subclass of //. The operations involved are the multiplication of elements 
in JT and the raising them to real powers. If 4,, Ay,..., A,, are ‘dimensionally 
independent” elements (7.e. a system of units) any plement PB. Se IT can be uniquely 
represented in the form 


(i) B= 1; Afr... APim 


where 7;, fj, ---, Pjm are real numbers and mt; 1S positive. The two theorems of 
part II are now combined and stated in thes following multiplicative form ap- 
propriate to an n-dimensional space J7: 

Dee D (As rn Apr eel bem dimensionally invariant and homogeneous 


function. If vt a mae for <n) are dimensionally independent and the P. have 
the form (1), then 


(ii) D(A,,..., Ay; B,..., PB) = pAb... Alm 


m? 
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where the positive coefficient p does not depend on the A’s, and the real exponents 
}; depend neither on the A’s nor on my, ..., Ty. 

This theorem, which the author regards as basic, gives the form of a function 
element of J/. It states what m does not depend upon; and later the author 
casually remarks that w does depend upon the positive coefficients 7,,..., 2, 
and the exponents #,, of the elements 2,..., P that appear in ®. But since 
the nature of this dependence on #;, is never stated, this theorem is no sub- 
stitute for the Pi Theorem. Indeed in the examples cited, the coefficient is 
a. function of the 7 dimensionless quantities 7,, 7, ...,%, Which appear in the 
P’s. Since there are m-+y7 arguments in ® and m units (the A’s), this suggests 
the Pi Theorem in the form first given by BuckiNGHAM in which the rank of 
the dimensional matrix is not considered. 


Drogor concludes his paper with a number of illuminating examples. One 
of these disposes very clearly of Lord RAyLeEicuH’s paradox; and others deal 


with the theory of quality control by sampling and particular solutions of partial 
differential equations. 
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Uber ein System partieller Differentialgleichungen 


ERWIN KREYSZIG 


Vorgelegt von S. BERGMAN 


1. Einleitung 


In der vorliegenden Arbeit wird eine Klasse von Systemen partieller Dif- 
ferentialgleichungen 


(1.1) Ua Oy yy F Pata So) 0; C= 412s 


Ky Xe 

untersucht. Hierbei sind gy, und g, nicht identisch verschwindende analytische 
Funktionen der jeweiligen (zunachst als reell angenommenen) Variablen. Die 
genannte, nachstehend definierte Klasse besitzt die folgende bemerkenswerte 
Eigenschaft: Fiir jedes System (1.1) aus dieser Klasse existiert eine unendliche 
Menge unabhangiger Partikularlésungen, deren jede gewdhnliche lineare Dif- 
ferentialgleichungen (mit algebraischen Koeffizienten) in jeder der vier in (1.1) 
vorkommenden Variablen befriedigt. Die Ordnung dieser Differentialgleichungen 
hangt nur von g, bzw. gz, ab, ist also fiir jede der genannten Partikularlésungen 
dieselbe. Die Menge dieser Loésungen bildet eine vollstandige Approximations- 
basis (vgl. Abschnitt 5). Auf diese Weise kann die Theorie der gewéhnlichen 
linearen Differentialgleichungen fiir die Untersuchung von Lésungen partieller 
Differentialgleichungen herangezogen werden. 


2. Formale Vorbereitungen 


Indem wir auch komplexe Werte der vier Variablen zulassen, kénnen wir die 


Funktionen q, und , analytisch ins Komplexe fortsetzen. Wir fiihren nun die 
Variablen - 
By = Xe Vy Bq = Xe — 4 V es a= 1,2, 


ein, die fiir komplexe Werte der urspriinglichen Variablen voneinander unab- 
hangig sind. Unter Benutzung von 


a) 1 0 0 0 AP ienc 5 

a = sla? ae = =F 

02 2 \ OX, ae 0%%, A ee By) ghee 
geht (4.1) in 
(2.4) Unw t+ H(%,2f)4=0, t,.2+ Ble, 24) u=0 
tiber. Hierbei ist 


ee 7 ok z Be as Jie ry comes 
CAC Ge Pe oh = u (ata 44, 2g + % By 74 
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und Pa psig 

18 Pica) i pe (5, a"), Ce = 11, Be 
Wir setzen fortan stets voraus, daB die Funktion F,(z,, 2*) in einem Gebiet 
B, x Bx regular ist, das den Nullpunkt (z,, z£) =(0, 0) enthalt. 

Es sei angemerkt, daB (2.1) zwei Systemen von je zwei partiellen Differential- 
gleichungen fiir den Real- und den Imaginarteil der komplexen Lésung u aqui- 
valent ist. Diese Gleichungen haben genau dann paarweise dieselbe Form, wenn 
die Funktionen F, fiir reelle Argumentwerte reell sind. 

Im folgenden kénnen wir die zu betrachtenden Systeme in der Form (2.1) 
zugrunde legen. 


3. Integraldarstellung der Lésungen durch Bergman-Operatoren 

Es existieren beliebig viele Operatoren, die analytische Funktionen einer 
komplexen Veranderlichen in Lésungen partieller Differentialgleichungen von 
zwei Veradnderlichen transformieren. Bergman-Operatoren [1]* besitzen die 
Eigenschaft, daB sich die Zuordnung der genannten Lésungen zu den analyti- 
schen Funktionen, die zunachst lokal gegeben ist, ins GroSe fortsetzen laBt. 
AuBerdem iibertragen sich eine Reihe grundlegender Eigenschaften der analyti- 
schen Funktionen auf die Lésungen. Deshalb kann man aus bekannten Ergeb- 
nissen der Funktionentheorie entsprechende Satze iiber partielle Differential- 
gleichungen mit analytischen Koeffizienten gewinnen. BERGMAN und SCHIFFER 
[2], [3] haben diese Theorie vor kurzem auf Systeme partieller Differential- 
gleichungen iibertragen. Die Verfasser bewiesen, daB sich jede im Nullpunkt 
regulare Lésung von (2.1) in der ate 


(3.1) (2%, 2%, 22,22) = Pl he) = J SAG (1,21, 4) Bs (22,22, te) x 


x [A Gal tt), $ 22(4 t3)} 
+ fotear(t tt), 3 23 (1 t3) } (1— 4) #(1— 4) * dtd, 


darstellen 1aBt. Hierbei sind die ,,Zugeordneten“ f, und f, der Lésung u im Null- 
punkt regulare analytische Funktionen zweier komplexen Verdnderlichen. Die 
Erzeugenden E, und E, des Operators P sind unabhangig von der Wahl der 
Zugeordneten. E, hangt nur von Ff, und E£, nur von Ly ab, vgl. (2.1). Die Dar- 
stellung (3.1) gilt in dem Gebiet B=B,xB,» XB, By, wobei die Faktoren 
dieses Produktes die im vorigen Abschnitt angegebene Bedeutung haben. 


4. Die Klasse © 


Wir untersuchen eine Klasse von Systemen (2.1), die wir als die ,,Klasse ©" 
bezeichnen und folgendermaBen definieren: 

Ein System partieller Differentialgleichungen (2.1) gehért der Klasse € an, 
wenn man in der Darstellung (3.1) der zugehorigen Losungen ,,Erzeugende vom 
Exponentialtyp™ 
(4.1) E, = exp Q,, Oe aes a? Sat) oS a=1,2, 

u=0 
wahlen kann. 


* Siehe Literaturverzeichnis am SchluB der vorliegenden Arbeit. 
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Die Systeme dieser Klasse lassen sich relativ einfach iiberblicken. Es ergeben 
sich naimlich die folgenden notwendigen und hinreichenden Bedingungen fiir die 
Zugehérigkeit eines Systems (2.1) zur Klasse C: 

Satz 1. Ein System partieller Differentialgleichungen (2.1) gehort genau dann 
der Klasse © an, wenn sich dessen Koeffizienten in der Form 


1 Oda 


q Oda ; ; * its 
(4.2) (a) Fey, ee) == ie oe eae oder (b) Fe ee) a Doe zx’ a AWG 
darstellen lassen*® ; hierbet ist 

[3(ma—1)] [4 ma] 


COGN Ron een Brey eee 1) uals 2a) = 2 Daa 2 a= 1,2; 
A=0 A= 


die beiden Koeffizienten a, . und b,, kénnen analytische Funktionen von zz sein, 
wahrend die tibrigen Koeffizienten konstant sind. 


Beweis. Wie aus [2, S. 85] hervorgeht, wurde (3.1) durch Verallgemeinerung 
eines Bergman-Operators gewonnen, mit dessen Hilfe sich die Lésungen einer 
einzelnen partiellen Differentialgleichung 


U, «+ F(z, 2*)u =0 
in der Form 


(4.4) u(e,2) = f EG, 2*,) /2(t —P)} (1 — 2)Mdt 


darstellen lassen, vgl.[Z]. Es gilt also fiir die Lésungen jeder einzelnen der 
beiden Differentialgleichungen (2.1) eine derartige Darstellung, und man kann 
die in diesen Darstellungen vorkommenden Erzeugenden EF als Erzeugende E, 
und £, in (3.1) verwenden, vgl. [2]. Demnach sind £, und £, so zu bestimmen, 
daB (4.4) mit EF, bzw. FE, als Erzeugender eine Loésung der ersten bzw. der zweiten 
Gleichung (2.1) darstellt. Nun sollen E, und E£, nach Definition der Klasse © 
die Form. (4.1) besitzen. Diese Wahl der Erzeugenden zieht Bedingungen iiber 
die Form der Funktionen / und fj in (2.1) nach sich, und zugleich ergeben sich 
Beziehungen zwischen den zunachst noch willkiirlichen Koeffizienten q, ,,(z,, 22) 
in (4.1) und den Funktionen /, und #. Durch Einsetzen der Darstellung (4.4) 
mit £ =E, in die erste und mit E = E, in die zweite der Differentialgleichungen 
(2.1) sieht man, da8 E,, «=1, 2, der partiellen Differentialgleichung 


an . A 
49) 0-005 1 Meese (PE tnn)mo, «nna 
gentigt. Man fiihrt nun den Beweis auf einem kiirzlich angegebenen Wege weiter, 
vgl. Kreyszic [4]: Jede in (4.5) auftretende Ableitung von E, laBt sich als 
Produkt aus FE, und einer Funktion darstellen, die Potenzen von it, und die 
Funktionen q,,, und deren Ableitung enthalt, vgl. (4.1). LaBt man den sich so 
ergebenden gemeinsamen Faktor E, weg, dann bleibt auf der linken Seite von 
(4.5) ein Polynom in ¢, iibrig. Die erhaltene Gleichung gilt identisch in ¢,. Durch 


* Dies ist so zu verstehen: A, kann z.B. die Form (b) und F, z.B. die Form (a) 
[oder auch die Form (b)] haben, so da® es insgesamt 4 MOglichkeiten gibt. 
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Nullsetzen des Koeffizienten einer jeden vorkommenden Potenz von i, erhalt 
man ein System von m, + 3 linearen partiellen Differentialgleichungen 2. Ordnung, 
das man nach der in [4, S. 915—919] angegebenen Methode lésen kann. Auf 
diesem Wege ergibt-sich die Aussage des vorliegenden Satzes und zugleich der 


Zusatz. Die tibrigen, in Satz 1 noch nicht genannten Koeffizienten pod Ciena) 
der Darstellung (4.1) haben die Form 


[amor] 5 
(4.6) J=%— >, a 2) see kousts: a = 12. 
al 
3(mo—1)] 


Bs 
Gu outa—(— 2)" (3°5...(2¢+4))* Do AA—1)...A—p+1)a,,2"4 


A= 


(wenn (4.2 a) gilt), 


Ge,ae1=0 (wenn (4.2b) gilt), w=1,2,...,[8m—1)], «= 1,2, 


[3 ma] 

Goon ee 4a) A 1) (An— 2) oni (A — ob AO, 425, 
A=Ké 

Te ee le eee ed, 


Damit ist die Klasse © vollstandig charakterisiert, und man kann fiir jedes 
System, das dieser Klasse angehért, zugehdrige Erzeugende E, und E£, sofort 
angeben. Das System 


(4.7) rl Ota ye. eq + U) = Uy + Zu =0, a=1,2, 
ist ein Beispiel eines Systems der Klasse ©, und die einfachsten Erzeugenden sind 


(4.8) ar sOD |e 27,\y 9) OAT Os 


wie aus Satz 1 folgt. 


5. Lésungen von Systemen der Klasse ©, die gewéhnliche Differential- 
gleichungen befriedigen 


Wahlen wir in (3.1) die Erzeugenden (4.1) und die Zugeordneten 


61) Atw=#4, fla)=A8", By.0,1=0,1,..., 


so erhalten wir eine unendliche Menge M unabhangiger Partikularlésungen der 
Systeme (2.1) der Klasse ©. Wir wollen zeigen, daB jede dieser Lésungen auBer 
(2.1) vier gewohnliche lineare Differentialgleichungen befriedigt. 


Zuvor vermerken wir die wichtige Tatsache, daB die Menge M in folgendem 
Sinne eine vollstdindige Approximationsbasis bildet: Jede Lésung eines der be- 
trachteten Systeme (2.1), die in einem konvexen (reell-vierdimensionalen) Gebiet 
B,,,,, das den Nullpunkt enthalt, regular ist, 1aBt sich in jedem abgeschlossenen 
Teilgebiet von B,,, durch eine endliche Linearkombination von Elementen aus 
M gleichmaBig approximieren. Dies folgt aus den Untersuchungen von BERGMAN 
und SCHIFFER [2] in Verbindung mit dem Approximationssatz von HAMMER- 
STEIN [4] fiir analytische Funktionen von zwei komplexen Verdnderlichen. 
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Um gewohnliche Differentialgleichungen zu gewinnen, die durch die Lésungen 
aus M befriedigt werden, haben wir diese Lésungen in gewissen (reell-zweidimen- 
sionalen) Ebenen des (reell-achtdimensionalen) %,y,%_V.-Raumes zu betrachten, 
vel. (4.4). Wir wahlen die vier Scharen derartiger Ebenen, die sich ergeben, 
indem man jeweils 3 der 4 genannten komplexen Variablen konstant annimmt. 
Dann haben wir den Vorteil, daB die Lésungen analytische Funktionen der nicht 
konstant angenommenen Variablen sind. Es gilt 

Satz 2. Jede Lésung U(x1, V1, %2, Yo) =U(%_ +11, %— ty» Xo +4 V2, X2—1 Ye) 
eines Systems (2.1) der Klasse © mit einer Zugeordneten von der Form (5.1) gentigt 
vier gewohnlichen Differentialglerchungen 


hj Biz 
G 


(5.2) > BO) sae 0, ee te ee 
2=0 


Us 


Hierbei ist v,=%,, Vg=Vy, V3 =%_, Vg—YVo, Und die jeweils nicht angegebenen 3 
der 4 Variablen sind konstant; die Koeffizienten haingen von der Wahl dieser kon- 
stanten Werte ab und sind rationale Funktionen von q,,, 6 =0, 1, --., My, =1, 2, 


und den Ableitungen dieser Funktionen, vgl. (4.1). Die Ordnungen der Gleichungen 
sind unabhangig von dem Wert der Exponenten B,y,o und t in (5.1); es gilt 


(5:3). Ree 2 (41), 7 Say ore Re Se th ee 


Beweis. Infolge der Form der Erzeugenden (4.1) und der Zugeordneten (5.1) 
hat (3.1) die Gestalt 


U= Uy Ug, UW = MyyMy2, Ug = Ug Uy2 
mit 
1 


OAC Seah ak (exp Q,) tha 4h,» Ai Gia, Ds 
=i 
wobei 
2 \B 2) p—4 2 \? Oy oe 
Aiea) =(B)0-AP, hele 4) =(2) ar 


fal. =(2)0-A4, feta (Fya—a 


ist. Es geniigt, eine der vier Differentialgleichungen (5.2), etwa die dem Wert 
let entsprechende, zu betrachten; fiir die anderen verlauft der Beweis analog. 
Die Funktion My, (%, 24) = %41(%,+79,, % —7y,) =U, (x, y,) befriedigt eine ge- 
wohnliche lineare Differentialgleichung 


~ 14 68 U 
(5.4) Ly (#43) =1,(U,,) = De (* 98) <= 0, I,=t, 9 =92=konst. 
tes 


Ox? 
von der Ordnung 7,< m,+1; die Koeffizienten dieser Gleichung sind rationale 
Funktionen VON %3,, @=0,1,...,m,, und den Ableitungen dieser Funktionen, 
wie aus [6, S. 808], Satz 2 folgt. Da in (5.2), (= 1), Wuibag., undsa. aiveniess 
und z konstant sind, so ist auch u, = W414. eine Lésung von (5.4). Entsprechend 
gentigt m2, (Zz, 21°) = U5, (x, y,) einer gewdhnlichen linearen Differentialgleichung 

Le} 


g ~ \ ce 
(5.5) Le(toy) = L, (U2) = » Iss Cas) iO H,=1, = y= konst. 


reMetee acre 
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von der Ordnung 7, < m,+1, deren Koeffizienten ebenfalls rationale Funktionen 
von q,, und den Ableitungen dieser Funktionen sind. Dasselbe gilt demnach 
auch ftir %,=4,;U25. So haben wir insgesamt drei Gleichungen, 


(5.6) U =U, + Ue, L,(u;) =0, Lg(%) = 0, 


aus denen wir durch Elimination von uw, und uw, und den Ableitungen dieser 
Funktionen eine gewéhnliche lineare Differentialgleichung fiir w von der Form 
(5.2) gewinnen kénnen. Da hierbei jede der Gleichungen (5.4) héchstens 
(2m,-+2)mal zu differenzieren ist, so ergibt sich die obengenannte Schranke fiir 
k,, vgl. (5.3). Entsprechend beweist man die Existenz der iibrigen Gleichungen 
(5.2). Man beachte dabei, daB Q,(z5, 23, t,) ein Polynom in ¢, vom Grade mg ist; 
deshalb sind die Schranken fiir k, und k, von denen fiir k, und k, verschieden. 

Zum. Beispiel ergeben sich fiir das System (4.7) bei Wahl der Erzeugenden 
(4.8) und der Zugeordneten (5.1) die Losungen u=u,-+ u,. mit 


maal B+ a)F Ot ala) (gy 024) 1 (44). 


Jy" (2Y" hae) h(a). 


Dies folgt unmittelbar aus der bekannten Darstellung 


Dp aT (o -- s\F le + 2) (4 


1 


[exp (tv cos ¢) (sin g)?"dg, R(n) > a ; 


0 


mn 


7 
Jy ”) Dy Jar (n + 4) 


der Bessel-Funktionen (vgl. NIELSEN [7, S.51]), indem man cos p=#, bzw. 
cos p=t, setzt. Die Differentialgleichungen (5.4) und (5.5) haben im vorliegenden 
Fall die Form 


L, (Uy) = Ua + On1 Ya + On 9 Mz = 0, a eae 


mit 


bar = (% a + 2i¢, v1) /(24 + 97) 


Gp — 6, 1 = 1,2. 


2 
bao = (x1 a y \/ i+ 9) 
wobei Striche die Ableitungen nach x, kennzeichnen. Hieraus gewinnt man fiir » 
eine gewohnliche lineare Differentialgleichung 4. Ordnung in *,, die sich in der 
Form 


| oe c d 0 0 
Li(u G c+d’ d 0 

| a € emcee dC 26 d= 0 
aie ea 1 0 
0 b30 bao + b2a1 boy 1 


darstellen 1aBt, wobei c=b,,—6,, und d=b,, —),, ist. 
4% 
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6. Folgerungen 


Wir schlieBen mit einigen Bemerkungen tiber Folgerungen aus den gewonnenen 
Ergebnissen. 


(1.) AuBer der bisher betrachteten Menge M von Losungen der Systeme (24) 
der Klasse & existiert eine weitere Menge unabhangiger Lésungen mit ahnlichen 
Eigenschaften. Diese Lésungen entsprechen den meromorphen Zugeordneten 


(6.1) f,@1,%2) = CRD a Pe Se is (21,22) = (4 — ae (= 'GNn 3 


Ag, Va At Date, 


wobei a,, b,,¢, und d, beliebige, von null verschiedene Konstanten bedeuten. 
Jede dieser Losungen befriedigt ebenfalls gew6hnliche Differentialgleichungen 
von der Form (5.2), deren Koeffizienten rationale Funktionen von 4,,, “=0, 
1,..., M,, &=41, 2, und den Ableitungen dieser Funktionen sind. Die Ordnung 
der (5.2), 7=1,2 bzw. 7=3,4 entsprechenden Gleichung betragt hochstens 
2(m,+3) bzw. 2(m,+3). Dies beweist man ahnlich wie den Satz 2 der vor- 
liegenden Arbeit unter Heranziehung des Satzes 1 in [8]. Durch Kombination 
dieses Resultates mit Satz 2 ergibt sich der 


Zusatz zu Satz 2. Lésungen von Systemen (2.1) der Klasse ©, deren Zuge- 
ordnete beliebige im Nullpunkt regulire rationale Funktionen sind, befriedigen ge- 
wohnliche Differentialgleichungen von der Form (5.2). 


(2.) Die in (6) angegebene Charakterisierung der Singularitaten einer Lésung 
einer einzelnen partiellen Differentialgleichung mit Hilfe gewohnlicher Differential- 
gleichungen laBt sich auch auf den Fall von Systemen (2.1) der Klasse € iiber- 
tragen. Ohne hierauf im Augenblick naher einzugehen, erwahnen wir die folgende 
bemerkenswerte Tatsache: Die Koeffizienten der gewdéhnlichen Differential- 
gleichungen (5.2) sind im allgemeinen Fall rationale Funktionen der q,,,, “=0, 
1,..., My, £=1,2, und threr Ableitungen. Sind nun F und £& rationale Funk- 
tionen, so sind die Koeffizienten von (5.2) sogar rationale Funktionen der jeweiligen 
Variablen; dies beweist man unter Benutzung von Satz1 und Satz2. Auf 
Grund des vorstehenden Zusatzes erhalt man demnach den 


Satz 3. Fiir Systeme (2.1) der Klasse © mit rationalen Koeffizienten F, und F, 
gilt: Die Singularititen der Lésungen mit rationalen Zugeordneten f, und fy [ugl. 
(3.1)] legen auf (reell-sechsdimensionalen) algebraischen Mannigfaltigkeiten des 
212 22¢-Raumes. 


(3.) Befriedigt eine Lésung einer partiellen Differentialgleichung eine gewéhn- 
liche Differentialgleichung, so kann man diese Lésung mit Hilfe der letzteren 
auch auferhalb des Giiltigkeitsbereiches ihrer Darstellung durch einen Bergman- 
Operator untersuchen. Dies hat BERGMAN [9] fiir den Fall einer einzelnen 
linearen partiellen Differentialgleichung erstmals gezeigt. Es ist zu erwarten, 
da8 sich diese Bergmansche Theorie mit Hilfe der vorstehenden Ergebnisse auf 
Systeme partieller Differentialgleichungen tibertragen 14Bt. 
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Indroduction 


Ce mémoire est consacré a l’étude des inductions électromagnétiques en rela- 
tivite générale. On y trouve exposée une démonstration du principe de FERMAT 
basée sur les propriétés des caractéristiques des équations de MAXWELL. 


Les champs et inductions électriques, magnétiques sont introduits comme l’ont 


fait GoRDON, WEYL, LICHNEROWICZ, a l’aide de deux champs de tenseurs anti- 
symetriques d’ordre 2: le tenseur champ électrique — induction magnétique 
H,,3 et le tenseur induction électrique — champ magnétique G,,. Des équations 


de liaison expriment les relations linéaires entre inductions et champs. C’est 
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l'ensemble des deux tenseurs H,,, G,g qui constitue le champ électromagnétique. 
Le schéma énergétique ot il intervient est un schéma <fluide chargé conducteur» 
dont l’étude a été faite antérieurement* et dont on rappelle les résultats essentiels 
sur l’intégration des équations du champ. 


En présence d’inductions, les variétés caractéristiques des équations de MAx- 
WELL ne sont pas identiques aux variétés caractéristiques des équations d’ EINSTEIN. 
On sait que ces derniéres sont tangentes en chacun de leurs points aux cénes 
élémentaires ©, de l’espace-temps. Les cones caractéristiques ©, pour les équa- 
tions de MAXWELL sont intérieurs aux cénes élémentaires. Les variétés carac- 
téristiques tangentes aux cénes ©, sont orientées dans le temps et les bicaractéris- 
tiques sont les géodésiques de longueur nulle de la métrique associée** 


roe Ln ees, B 

aSe=— (ep (1 =) is ug) dx dx 
ou u* désigne le vecteur vitesse unitaire d’univers et les scalaires ¢, w le pouvoir 
diélectrique et la perméabilité magnétique en chaque point du milieu considéré. 


Or les variétés caractéristiques des équations de MAXWELL jouent le réle de 
surfaces d’ondes électromagnétiques et les bicaractéristiques, celui de rayons 
électromagnétiques correspondants. On est conduit a introduire naturellement 
la variété riemannienne &, définie par la variété différentiable portant l’espace- 
temps et munie de la métrique associée ds”. 

Les équations de MAXWELL peuvent s’exprimer directement dans cette variété 
ou elles affectent une forme simple symétrique de celles de la théorie électro- 
dynamique de Lorentz. Les rayons électromagnétiques sont géodésiques de 
longueur nulle de &,. L’étude géométrique des rayons électromagnétiques dans 
Vespace fournit l’énoncé du principe de FERMAT. 


Pour faire cette étude, nous avons commencé par généraliser au cas du fluide 
parfait chargé conducteur la notion de mouvements permanents lée a |’existence 
pour l’espace-temps %, d’un groupe connexe a un paramétre d’isométries globales 
a trajectoires orientées dans le temps et ne laissant invariant aucun point de B,***. 
Si le mouvement du fluide considéré est permanent, il est défini dans la variété 
¥, un groupe d’isométries induites par le groupe de l’espace-temps. Nous sommes 
dans un cas ot! la«méthode de descente» de LICHNEROWICZz**** s’applique. En 
projetant les géodésiques de longueur nulle de &, sur la variété quotient de &, 
par la relation d’équivalence définie par son groupe d’isométries, nous obtenons 
un théoréme qui généralise le principe de FERMAT en relativité. Ce théoréme 
est valable pour un milieu en mouvement permanent quelconque. En particulier 
dans le cas d’un espace-temps sans gravitation de MInKowskI, il fournit comme 
conséquence une démonstration de la formule relativiste de la composition des 


vitesses. 


x Etude électromagnétique et thermodynamique d’un fluide relativiste chargé 
[Jour. Rational Mechanics and Analysis 5, No. 3, 473 —583 (1956)]. 
xx GORDON a trouvé cette métrique par voie algébrique. 
xxx TICHNEROWICZ, A.: Théories relativistes de la gravitation et de l’électromagné- 
tisme, chap. IV, III, pp. 83—90. Masson 1955. 
xxkxk Tbid., Livre II, chap. 1%. $ 
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Notations employées 


é Se ee 
0,= Eas Cup = apisaek © dérivées partielles 
We : dérivée covariante 

&,. 6, «»», tout, indice grec.== 0,4; 2,3. 


4,4, ..., tout indice latin = 1, 2,3. 


I. Inductions électromagnétiques. 
Intégration des équations de MAXWELL 


1. Les équations de Maxwell dans l’espace-temps %, 


Soit dans la variété espace-temps %, de la relativité générale, munie de la 
métrique d’univers 


(1.1) Ast == 95nd bs Nh = OG oe) 


un domaine  occupé par une distribution matérielle schématisée sous forme de 
fluide — champ électromagnétique. % désigne le vecteur vitesse unitaire en 
chaque point % de . On appelle repéve propre en x, un repére orthonormé dont 
le premier vecteur V coincide avec % et dont les trois autres vecteurs Jt 
orientés dans l’espace sont normés par la condition 


gap VO*VOG— — 4, 


Les phénoménes électromagnétiques sont caractérisés par les deux champs 
de tenseurs antisymétriques d’ordre 2: le tenseur champ électrique — induction 


magnétique H,, et le tenseur induction électrique — champ magnétique Gup, 
dont les composantes relatives 4 un repére propre au point x considéré ont pour 
valeurs 
0 | Sr danine! Dre aI Sk- 0 De Dee IDE 
cot vec 0) B, —B —D 0 =H, 

(H,,,) zat eu 3 2 (G,.3) = 1 3 2 
Wass” eae 0 B, =D), pails 0) if 
=a.) Bede D. 10) = D3ny) Hot 

et vérifient les relations 
(1.2) Gi=eHy;, H,;—=uG,,  (,7=4,2,3) 


ou les scalaires ¢ et  représentent respectivement le pouvoir diélectrique et la 
perméabilité magnétique du milieu considéré. 


Nous introduisons les tenseurs adjoints 


(1.3) He — oT? eee G2 — 4 72079G , 


oU Hagys est le tenseur complétement antisymétrique attaché a la forme élément 


de volume de &,. Les relations (1.2) peuvent alors s’écrire sous la forme in- 
variante 


(1.4) Gapu* = cH, gu, 
(1.5) [Gap =H, yw 
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valable dans un systéme de coordonnées locales quelconque. Ces relations sont 
appelées les équations de liaison. 


Les deux champs de tenseur H,, et G,, doivent satisfaire aux équations de 
MAXWELL qui s’écrivent 
(1.6) UH 0 
(7) VG sane: 


J® est le vecteur courant électrique. Le premier groupe des équations de MAXWELL 
peut encore s’écrire 


£ bro 7 Hy, = 0. 


Il exprime qu’il existe localement un champ de vecteur gy, dont H,; est le rota- 
tionnel, c’est-a-dire 


Ay = Oy, Op Og Pa 
L’évolution du champ électromagnétique est déterminée, si l’on connait J? 
c’est-a-dire la distribution de ]’électricité. On peut supposer pour un milieu non 


> 
conducteur que / est colinéaire au vecteur vitesse 
PH0us 


6 s’appelle la densité propre de charge électrique, et le courant électrique est 
un courant de convection. Plus généralement, on est conduit a faire ’hypothése 


J* = 6u* + ou, H 


ou 6 est encore la densité propre de charge électrique et o un scalaire carac- 


térisant la conductivité électrique du milieu. ir posséde alors une composante 
colinéaire 4 % et une composante 1% =ou,H®* orthogonale 4 uv. La premiére 
représente le courant de convection et la seconde, le courant de conduction qui 
satisfait 4 lhypothése d’Oum. 

Les équations (1.4), (1.5), (1.6), (1.7) constituent les équations de |’électro- 
magnétisme en présence d’inductions dans la matiére. Dans le vide, on a l’égalité 


(1.8) | 
et les équations de MAXWELL deviennent 

(1.9) V,H=! = 0, 
(1.10) Gt? 2-0. 


tandis que les équations de liaison se réduisent a 


(1.11) CoH eee. 


eels 
Feel . > FS . . 
Nous désignerons dans la suite par & et @ les vecteurs qui figurent au pre- 
mier membre des équations de MAXWELL, soit 


(1.12) 6° =V,H*?, 
(1.13) GP = V,G**. 
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On démontre que leurs divergences sont nulles: 
(1.14) Vola 0), 
(45) ow, 


Ces deux équations sont appelées les conditions de conservation relatives aux 
équations de MAXWELL. Elles expriment la conservation de l’électricité. Ainsi, 


on tire de (1.7) et (1.15) V, J*=0 


2. Expression des G,, en fonction des H, 


Les équations de liaison (1.4) et (1.5) traduisent le caractére linéaire des 
relations entre inductions et champs. Elles montrent que les deux champs de 
tenseurs H,,, et G,g ne sont pas indépendants I’un de l'autre. On peut exprimer 
les G, en fonction des H,,. 


En effet, A partir de (1.4) nous pouvons former l’égalité 
(2-4) (Gap Uy -+ Gp, Uy) ue = &(H, pu, + Hg, Uq) u?. 
D’autre part, (1.5) peut s’écrire sous la forme équivalente 
1 
Gap My + Gey Ue + Gyq Mg = a (Hyg Uy + A, p Uy + Hq Up) 


vraie pour tout groupe de valeurs données a4 «, 8, y. Par multiplication contractée 
de cette relation avec w’, puis en retranchant (2.1) de l’égalité ainsi obtenue, 


nous avons 


Sas i (Hg Uy, a Hp, Un “a oie Us) wee e(H. pu, ae Ht, Uy) uP. 


ya (ed 
Nous en déduisons 
(2.2) Gap =) Hag + ——** (Hyg tty — Hyp tt). 
C’est la relation cherchée. En composantes contravariantes, on a 


(2.3) Gee “ He? eee (H°%u,u® — HP uw, u*). 


3. L’intégration des équations de MAXWELL-EINSTEIN 


Le champ électromagnétique (H,,, G,,) et la métrique ds? =g,,dx%dx° sont 
liés par les 6quations de MAXWELL-EINSTEIN. Si le milieu meublant le domaine 


D4 considéré est schématisé sous forme de fluide partait chargé et conducteur, 
les €quations d’ EINSTEIN sont * 


(3.1) Sap = Rag — FR bas =v Typ; 


Ln3= (o 1p) Uy Up - P 8up (Ux Ip Ug Ia) rat pee (1 — &p) Tao U° Ug, 
Tg = Pee (Geel) 7 (Coie Li 


(3.2) (———— 20, B (ge — ue Uy) 


x Cf. Etude clectromagnétique et thermodynamique d’un fluide relativiste chargé 
[Jour. of Rational Mechanics and Analysis 5, No. 3, 473 —583 (1956) ] 
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ou ? est la pression et # la température en chaque point du fluide. Les équations 
de MAXWELL sont 


(3.3) _ &= kn ?V, Hy, =0, 

(3.4) Dy = 8°° V,, Gog = 6g + ou" H,, 9. 
A ces équations s’ajoutent les équations de conservation 
(3.5) yal. 

(3.6) V9" = cota, 3 ——u*d, 0, 
(3.7) VAdul + ou, H°*) = 0. 


Les scalaives x,c,l, &, 4,0 sont supposés donnés; ils caractérisent le fluide 
envisagé qui admet de plus l’équation d’état 


(3.8) 0 = 9(,%). 


Les variables de champ sont constituées par l’ensemble 


9 (Sx. ELIS dv, uU™, p, )) 
oti le vecteur u* est normé 


(3.9) Syphwre= +1. 


Le probléme qui se pose est celui de |’intégration des équations du champ. 
On peut l’étudier au moyen d’une analyse du probléme de Cauchy. Pour cela, 
on se donne sur une hypersurface © orientée dans l’espace, représentée localement 
par 

C10), 
les valeurs des quantités 
6 (S43 , Go 8a B: VeET es a, 60) 


et on se propose de déterminer les divers champs # (g,,, Hug, 8, u%, p, 6) en 
dehors de © dans son domaine d’existence. I] suffit d’étudier la possibilité de 
calculer sur © les valeurs des différentes quantités introduites et de leurs 
dérivées successives. 

L’espace-temps %, étant une variété différentiable de classe (C?, C? par 
morceaux), on supposera donc g,, de classe (C1, C? par morceaux), H,, de classe 
(C°, C2 par morceaux) et # de classe (C?, C+ par morceaux). 

Pour g°°=-0, les équations d’EINSTEIN sont équivalentes a l’ensemble des 
deux systémes 


(3.10) Rip = — 5 Op bij + y= x (Zi; — 3 T bi) 


(3.11) S2=7[(@ +P) WM — Ba — (Udy + Hyg?) + Oy — (1 — ef) Do? Hy] 


ou les F;; et S° ont des valeurs connues sur ©. Les équations (3.11) jointes au 
caractére unitaire de u* et a l’équation d’état fournissent les quantités , u*. 
(3.10) déterminent alors @9g;; si g°°=-0. 
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Les équations de MAXWELL sont équivalentes a l'ensemble du systéme 


6412) 9,= (0° E(t Sep) (w®)2) i pe ; (2° =(1— en) # u') 0)H;;+ @D; 


=6u,+ou'H,;, 
(3.13) fF Oo 0 
et des deux identités 
(3.14) P = du + ou, H*°, 
(3.15) 6) = 47/10, =0 


ot les @, et Y ne dépendent pas des 0,H,, mais dépendent des ¢)u* tandis que 
la quantité G° ne dépend pas de 4),H,, ni de 6yu*. (3.15) exprime qu'il existe un 
potentiel vecteur local pour H;; sur S. L’équation (3.14) détermine 0 si w° +0. 
Pour avoir @)H,,, il faut chercher a déterminer d’abord les dérivées de u*: soit 
d,u*. Cette détermination se fait simultanément avec celle de df, Gy9% au 
moyen des équations de conservation relatives aux équations d’EINSTEIN (3.5) 
auxquelles on adjoint le caractére unitaire de wu, l’équation de conduction ther- 
mique (3.6) et l’équation d’état. Puis la derivée @)6 se calcule par ]’équation 
de conservation du courant électrique qui peut s’écrire 


Oe 
ou Q dépend de @)u* mais ne dépend pas de @)H,.,. 
Les 0,u* étant calculées, on les porte dans (3.12) et (3.13) qui fournissent 


enfin les 09H, si 
g°° — (1 — e44) (m9)? + 0. 


Si hypersurface S portant les données de Cauchy @ n’est pas exceptionnelle, 
il résulte des équations (3.10), (3.12), (3.13), (3.5), (3-6), (3.7) que les quantités 
Ooo8i;» Colas, Cpo%, AU, AP, G5 sont bien déterminées et nécessairement 
continues a la traversée de l’hypersurface GS. Les mémes conclusions s’étendent 
aux dérivées d’ordre supérieur de ces quantités si on suppose les données déri- 
vables 4 un ordre plus élevé que celui de l’hypothése. 

La détermination des quantités précédentes ne fait pas intervenir les équations 
(3.11), (3.14), (3.15). Or celles-ci ne contiennent aucune dérivée oblique des 
données de Cauchy; celles-ci sont donc astreintes 4 vérifier sur la variété © les 
trois équations (3.11), (3.14), (3.15) ou leurs équivalentes 


(I) Q=S8? yal Ren N(0 
Po = @ — (bu + ou, H*) =0 
C=) 


ou l’on a posé 


Cap a5 =X lags = 9, — (ou, ow H,,). 


Considérons maintenant un ensemble # (Sp. 15,0, u%, p, 6), solution des équa- 


tions (3.10), (3.12), (3.13), (3.5), (3.6), (3.7), correspondant a des données de 
Cauchy @ satisfaisant sur G aux équations ( I). En vertu du caractére conservatif 
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des premiers membres des équations de MAXWELL-EINSTEIN et des équations 
de conservation (3.5), (3.7), on a 


O50) fe Vo bee0) Ai Veo 0. 
Ces identités se réduisent en vertu de (3.114), (3.12), (3.13) aux équations 
°° 8,05 = AP 2,03 + BLO} | 
QP) = — C0, Po — (0,C' + 1%5C*) Pe 
0,6 = — 1) 


ou les 43°, BE, C* sont des fonctions continues. Ces équations sont linéaires et 
homogénes par rapport aux inconnues Q}, P®, °. Comme 0° = P® = &° =0 sur 
©, elles n’admettent pas d’autre solution que la solution identiquement nulle. 
Il en résulte que si les équations (I) sont vérifiées sur S par les données de Cauchy, 
elles sont également vérifiées dans tout le domaine d’espace-temps considéré par 
la solution #(g,,, H,,, 0, u%, p, 6) des équations du champ. 

Le probleme de l’intégration des équations du champ consiste finalement dans 
le choix des données de Cauchy rendant compatibles les équations (3.11), (3.14), 
(3.15) qui permettent de calculer u*, #, 6, puis dans l’intégration du systéme des 
équations (3.10), (3.12), (3.13), (3.5), (3-6), (3.7) qui permettent d’étudier 1’évo- 
lution des champs # (g,,, Hig, 0, u%, p, 0). Si les données du probléme étaient 
analytiques réelles, on pourrait établir a l’aide du théoréme d’existence de CAUCHY- 
KOWALEWSKI pour les équations aux dérivées partielles, qu’a un changement 
de coordonnées prés conservant point par point l’hypersurface © et les données 
de Cauchy sur G, le probleme admet une solution analytique réelle et une seule, 
solution dont nous connaissons le développement suivant les puissances de x?. 
Sous des hypothéses de simple différentiabilité, la méthode de Mme Foures 
permet d’établir l’existence et l’unicité de la solution. 


4. Les variétés caractéristiques B)" 


Sur les équations (3.12), on voit que si hypersurface © portant les données 
de Cauchy est telle que sur S 


gio (1 em) (m?)* = 0, 


les Cérivées 6) H,; du champ électromagnétique peuvent étre discontinues a la 
traversée de ©. Il peut exister une infinité de solutions distinctes des équations 
de MAXWELL correspondant aux mémes données de Cauchy. La variété © est 
une variété caractéristique pour les équations de MAXWELL. Une telle variété 
sera désignée par Bf". 

Dans un systéme de coordonnées locales arbitraire quelconque, les variétés 
caractéristiques 83 définies par f(x*) =0 sont les variétés satisfaisant a l’équation 


(4.1) (9%% — (1 — eu) u* 0") 0, f, f = 0. 


Ces variétés a la traversée desquelles peuvent se produire des discontinuités du 
champ électromagnétique, constituent l’extension relativiste des fronts d’ondes 
électromagnétiques classiques. Nous supposerons que ces fronts d’ondes sont 
orientés dans le temps, ou a la rigueur tangents au cone élémentaire ds? =0 
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de %4; nous verrons que cette hypothése est bien en accord avec les exigences 
physiques relativistes. S’il en est ainsi, 


A,f = gO, f Op f = (1 — ep) (Ww, f)? <0. 
On en déduit 


(4.2) EEN Leas Be 


Ceci posé, la généralisation de l’hypothése d’ HUGONIOT permet d’évaluer ce 
qui constitue ici la vitesse de propagation des ondes CRctrOna gn consi- 
dérécs. Pour cela considérons deux surfaces d’ondes voisines (¥3‘)> et (B3')» 
définies par les équations 


Hx) =0, fe)=0 
et prenons # pour infiniment petit principal. 
La ligne de courant issue du point x de (W347), coupe (B37), en un point défini 


aux infiniment petits d’ordre supérieur prés par x+y, 4 étant donné par la 
relation 


(4.3) HU 0, | =e 
Soit 7% le vecteur normé (72 —1) normal en x a la surface d’onde (¥3),. Ila 
pour composantes covariantes en x 
Gf 
4.4 n SS 
(4.4) “Te 


La trajectoire orthogonale des %}/ issue de x coupe (%4‘), en un point qui, 
a des infiniment petits d’ordre supéricur prés, s’écrit «+ 7,7, y, étant déterminé 
par la relation 


mn d,f=9 
On cn déduit 
ov OY — 628 a, f Opi —# 
4.5 YY, = a SS oe ve B — a 
oy nah Of 710 GR Of PPL gat agp o,f 


Introduisons le vecteur ¢ = nu — Ta He En vertu de (4.3) et (4.4), on a 


n (u -72) =n 
et 


b= (qi — 1H) =H -R) +m =0. 


> 
Le vecteur ¢ est donc tangent a la surface d’onde. II est orienté dans le temps, 
car son carré 


mo = (t)? = 9? — i — 2m -%) = 9? + 73 
est positif. 

Le vecteur 7¥ apparaft ainsi comme la somme de deux vecteurs, l'un ortho- 
gonal a la surface d’onde et orientée dans espace, l’autre tangent a cette surface 
et orienté dans le temps. La vitesse de propagation V de l’onde se trouve définie 
comme la limite du rapport des modules de ces deux vecteurs, soit 


ph, 
No 


WY = laa 


0—0 
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On a ainsi 


2 
V2 = lim ; 
930 ne 


soit, en remplacant 7, et 7) par leurs valeurs 


/ 


La vitesse de propagation des ondes électromagnétiques est done 4/]/eu. 
Cette valeur appelle deux remarques: d’abord elle généralise la valeur obtenue 
en électromagnétisme classique. De plus, dans nos hypothéses (eu = 1), la 
vitesse de propagation V est inférieure 4 une vitesse limite c=1; cette valeur 
limite coincide avec la valeur de la vitesse de propagation des ondes électro- 
magnétiques dans le vide (eu“=1). 


II. Etude des caractéristiques 
5. Définition d’une métrique associée 


L’intégration des équations de MAXWELL fait intervenir le champ de tenseur 
contravariant symétrique 


(5.1) g~F — p%B — (4 — ey) ut uP 
dont la forme quadratique associée représente la forme caractéristique des 


équations de MAXWELL. L’étude des variétés caractéristiques 83! de celles-ci 
devient plus suggestive si l’on introduit la métrique riemannienne 


(5.2) ds" =P aa 


dont la matrice des coefficients (%,,) est la matrice inverse de la matrice (g*?). 
On obtient facilement 


(5.3) 20 = Sap — (1 a) Map 


T2| 


en effectuant les calculs en repére propre. Et, si g et g représentent respectivement 
le déterminant de la matrice (g,,) et celui de la matrice (g,,), on a la relation 


(5.4) g = eug. 


La métrique ds? sera dite métrique associée. Elle joue un réle fondamental 
dans l’étude des variétés caractéristiques des équations de MAXWELL. La métrique 
d’univers ds? =g,,dx%dx* est du type hyperbolique normal. Rapportée a un 
repére propre, elle prend la forme canonique 


dst= (w°)? eie (w1)? oes (w?)? as (w?)? 


ot les (w*) sont un systéme de formes de Pfaff locales linéairement indépendantes. 
La métrique associée prend elle-méme la forme 


as? = (5p = (1 —_ a Uy X%9| w* ow! 
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ot 6, =0 Si a==B, dog= +1 et 6;;= —1 et u=1, 4;—0 dans le repére propre. 
On en déduit 


= wm? \2 
ds? = =) — (wl)? — (@?)? — (8), 
Eft 
ce qui montre que la métrique associée est également de type hyperbolique 
normal. z 
Nous désignerons dans la suite par &, la variété riemannienne définie par la 
hice se : ieee a 
variété différentiable portant &, et munie de la métrique associée ds. Et nous 
distinguerons par une barre supérieure les quantités définies relativement a B,. 
Nous appellerons céne élémentaire associé ©, en un point ~, le cone réel de direc- 
tions tangentes a %, défini par l’équation ds? =0. 


6. Etude des bicaractéristiques 


Dans l’espace riemannien &,, les variétés caractéristiques des équations de 
MAXWELL, définies localement par f(x*)=0, sont solutions de l’équation aux 
dérivées partielles du premier ordre 


(6.1) A,f = 8° 0,f Onf = 0. 


Elles ont tangentes en chaque point au cone élémentaire associé ©,. Les cénes 
élémentaires ©, de &, sont donc cones caractéristiques pour les équations de 
MAXWELL et celles-ci admettent pour variétés caractéristiques les variétés tan- 
gentes a ces cones. Mais dans l’espace-temps %, les cénes caractéristiques des 
équations de MAXWELL sont en général différents des cénes élémentaires €, 
(ds*=0). Ils ne coincident avec ces derniers que dans les régions vides de matiére. 

Une variété caractéristique B37, c’est-a-dire une solution de (6.1), peut étre 
engendrée au moyen des bandes caractéristiques de (6.1). Une telle solution 
peut étre engendrée au moyen des bandes de &, constituées chacune par l’en- 
semble d’une courbe & et d’une famille 4 un paramétre de 3-plans élémentaires 
tangents a ces courbes. Les courbes &9 sont appelées les bicaractéristiques des 
équations de MAXWELL. 


Pour les déterminer, posons 
2H(x", y,) = 8%? y, Vp 
et considérons |’équation aux dérivées partielles 
(6.2) Af S2h a= GC 


ou C est une constante arbitraire. Relativement aux variables x*, f, Yg, les 
bandes caractéristiques des équations de MAXWELL (3.1) et (3.2) sont données 
par les solutions du systéme différentiel 


axe eye ak af dV o Vs __ d 
OH CERRO Rim ig ee 
OV OYs 8x ax8 


qui satisfont a l’intégrale premiére 


2A yh aC 
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pour la valeur C de la constante. Si l’on introduit la variable auxiliaire u, les 
fonctions x*(u), yg(u) sont données par le systéme canonique 


(6.3) dx oH AVy oH 


du Ovne du ox 


relatif a la fonction hamiltonienne H(x*, y,,). Le premier groupe des équations 
(6.3) s’écrit explicitement 


*a a8 OY, Ae ax 
(6.4) %% — BrP yn (z = ia) : 
Inversement 
(6.5) Vp = Bap”. 


Cela posé, les solutions x*(u) de (6.3) sont extrémales de la fonction lagran- 
gienne L définie par 2L = Ek 8 
puisque, par passage des variables (x%, %”) aux variables canoniques (x*, Ye) qui 
leur sont liées par (6.4) et (6.5), on a entre H et L la relation classique 

H=2° 2 T=. 


Ox 


Ces solutions sont les extrémales satisfaisant 4 l’intégrale premiére 
(6.6) ZC: 


pour la valeur C de la constante. Or, d’aprés l’existence de cette intégrale pre- 
miére, les extrémales ainsi définies sont aussi les extrémales de 


V2L = Vane? 


satisfaisant 4 (6.6). Il en résulte que les x*(u) définissent des géodésiques de By. 
Si C=0, le systéme différentiel aux caractéristiques de (6.1) admet l’intégrale 
premiére {=const. et les variétés B/ peuvent étre engendrées par les bandes 
de &, définies par les géodésiques de longueur nulle &, le 3-plan associé étant 


le plan tangent au cone élémentaire ©, le long de la tangente a &y. 


Nous avons démontré le théoréme: 


Théoréme. Les bicaractéristiques des équations de MAXWELL sont les géo- 
désiques de longueur nulle de la variété riemannienne B, munie de la métrique 
associée iA=Z,,dedx. 

Dans le langage de la théorie de la propagation par ondes, les variétés caracté- 
ristiques &3/ jouent le rdle de surfaces d’ondes électromagnétiques. Les bi- 
caractéristiques 2, sont les rayons électromagnétiques associés. Nous pouvons 
donc énoncer le résultat suivant: 

Théoréme. Dans un milieu isotrope de constantes diélectrique et magnétique 
é, « variables, les rayons électromagnétiques peuvent étre considérés comme géo- 
désiques de longueur nulle de l’espace riemannien 8, muni de la métrique 


is = @ “3 1 Om 
as" = 8.4% dn? = (Eup i= i) Matty) 4 dx? 


ou gy3 est le tenseur métrique fondamental et u, le vecieur vitesse wnitarre d’univers 
définis en chaque point du milieu. 
Arch, Rational Mech, Anal., Vol. 1 
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7. Les équations de 1’électromagnétisme dans la métrique associée 
Définissons sur B, le champ de tenseur antisymétrique H,,, tel qu’en chaque 


point (x*) 


Hye = Hea 
En élevant les indices de H,, l'aide du tenseur métrique de ¥B,, nous avons 
HP — ge pF FT, = (g2* — (1 — es) uu") (g7? — (1 — eu) uu") H,, 
d’ot 
He? — o%@ pF 4+ (1 — em) (g°% Hy, uu? — 2°? Hue u*). 


En comparant cette égalité avec la relation (2.3) donnant l’expression de Gee 
en fonction de H,,, nous voyons qu’en chaque point (%*) 


Gxb— 1 Fab. 
bb 
L’étude du probléme de Cauchy relatif aux équations de MAXWELL dans &, 
nous suggére d’écrire les équations suivantes dans %: 
(7A) Vee — 0% 
(7.2) V,,G*? — J8 


ot G,,, est un tenseur proportionnel 4 H,,. Le premier groupe (7.1) s’écrit encore 
7, fee — VE 0080, ,=0. 


I] exprime qu'il existe localement un champ de vecteur g, dont H,, est le rota- 
tionnel. Comme g=eug--0 et que H,,—=H,, au point considéré, on voit que 
les Equations (7.1) sont équivalentes aux équations de MAXWELL du premier 
groupe (16) dans ¥,. Et nous pouvons identifier les deux potentiels vecteurs 
Pa Ot Pa: 

Quant aux équations (7.2), on peut les écrire 


ou 


1 ya Ge ye 
aera Oy (Viel) = JT . 

Viel Yeu) Yeu 

En comparant cette équation avec les équations de MAXWELL du second groupe 
(1.7) dans &, écrites sous la forme 


1 
——— 0, (/|¢| G2") = J*, 
Viel 
on voit que les équations (7.2) sont équivalentes aux équations (1.7) si l’on prend 


(Geld = Veu GP et J® = Veu J°. 
On est conduit a poser 


(7.3) Ch ie ee 
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En particulier, si l’on fait l’hypothése 


J’ = bw 4+ ou, H%*, 
on vérifie que 

J’ = 0 + 0 t,he 
ol u* = d x*/ds. On observera que les scalaires ¢, uw, 0, 6 sont les mémes dans &, 
et By. 


L’étude précédente montre que les équations de MAXWELL peuvent s’exprimer 
directement dans la métrique associée comme 


ou G,, est lié a H,, par la relation simple 


Cre 2 ap 
B Lu B 
et ott le vecteur /® satisfait a l’hypothése 
Ji = d? + ot, H*. 


On retrouve facilement a partir de ces équations, les résultats établis aux para- 
graphes précédents. 


Ainsi, la formulation de la théorie électromagnétique de MAXWELL se fait 
de facon équivalente dans %, et dans ¥,. Il convient d’observer que les équations 
de MAXWELL conduisent en Physique classique a |’étude de l’opérateur linéaire 
hyperbolique du second ordre 

ATG o? oO Oo 
VE Or Ox® Oy2 Az?” 


Cet opérateur reste invariant par le groupe de LORENTZ, c’est-a-dire le groupe 
de transformations qui laisse invariante la forme quadratique de différentielles 


V2d?—dx—dy—d2 
et celle-ci n’est que la traduction en repére propre de la forme métrique associée 
qui s’écrit ici 
2 
d3? = ~ d@—dxt—dy?— dz. 
El 
On voit que 


Ve 


El 


La quantité eu qui se présente dans notre étude peut étre interprétée comme 
Vindice de réfringence du milieu considéré. Nous poserons 


ep 


m est un nombre positif = 1 sans dimensions. 
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III. Mouvement permanent d’un fluide partfait chargé 
8. Espace-temps stationnaire dans un domaine 


Considérons un domaine déterminé D, 4 quatre dimensions de &, et supposons 
que la variété riemannienne définie par D4, munie de la métrique d’univers 


ds Sg Ha 


admette un groupe connexe d’isométries globales 4 un paramétre, ne laissant 
invariant aucun point de D, et dont les trajectoires z soient orientées dans le 
temps. Nous supposons de plus que: 

a) les z sont homéomorphes a la droite réelle H; 

b) l’on peut trouver une variété différentiable a trois dimensions 23, satis- 
faisant aux mémes hypothéses de différentiabilité que %,, telle qu'il existe un 
homéomorphisme différentiable de méme classe de 2, sur le produit topologique 
D3 R dans lequel les z s’appliquent sur les lignes facteurs. 

Dans ces conditions, nous dirons que l’espace-temps riemannien YW, est 
stationnaire dans D,. Les trajectoires z sont appelées lignes de temps. La variété 
Dz quotient de D, par la relation d’équivalence définie par le groupe est appelée 
espace. 


> 


Soit € le générateur infinitésimal du groupe d’isométries. Aucun point de D, 


n’étant invariant, a est +0 en tout point de 2,. On sait que ce vecteur satisfait 
aux équations de KILLING 


(8.1) X £53 = Vers V5 é.=0 


ou X désigne l’opérateur de dérivation de Lie relatif au vecteur = 


On peut définir dans D, des systémes de coordonnées locales (x*), dits adaptés 
au caractére stationnaire, de la maniére suivante: les (x') sont un systéme de 
coordonnées locales arbitraire de 23. La donnée des (x‘) détermine une ligne 
de temps. Pour déterminer un point sur cette ligne, on se donne la variété 
«® = const. a laquelle il appartient, ces variétés étant les variétés homéomorphes 


, Py te ns ) , . ee . 
a Q; définies par l’homéomorphisme b) et telles que les composantes de €& soient 

&=1, &=0. 
Dans ces systémes de coordonnées adaptés, on a 

ea ma Sc0 

et les équations de KILLING se traduisent comme 

ai 

X 843 = Oo Sap = 0. 


Ainsi, les g,, sont indépendants de la variable x°. 


Dans la suite, nous n’introduirons que des systemes de coordonnées locales 


adaptés. En effectuant la décomposition en carrés de la forme quadratique 
fondamentale 


ds? = gi gd x* ax 
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a partir de la variable directrice dx®, nous obtenons 


(8.2) dst=" (,,dx*)2+ de 
‘ S00 
ou 
ine atic] cer (iy eee OFS 07 fh aie 
(8.3) dS? =, dx dx = @ Fe jax dx! 


est indépendante de x et définit sur D, une métrique riemannienne définie 
négative. 


9. Mouvement permanent d’un fluide parfait chargé 


Considérons un fluide parfait chargé et conducteur en mouvement dans un 
domaine 2,. Le mouvement de ce fluide est dit permanent si l’espace-temps 
riemannien associé ¥%, est stationnaire dans le domaine ®, et si le groupe d’iso- 
métries laisse invariants le vecteur vitesse unitaire u*, la pression p, la tem- 
pérature 7, la densité de charge électrique 6, le vecteur courant de chaleur q, et 
le champ électromagnétique induit H,,, G,g. 


Les quantités g,¢, H.5, Gug, u*, 0, 9%, b, 0 restent donc constantes le long des 
trajectoires d’isométries ou lignes de temps. Dans un systéme de coordonnées 
adapté, les égalités précédentes deviennent 


Oo Gap = 9% Hap = % Gap = Oy U = 0 B= % Ix = MP = 0 = 0; 
les g48, Hag, Gug, U*, B, %,, P, 6 ne dépendent pas de la variable +?. 


Considérons un mouvement du fluide envisagé tel que 

a) Vespace-temps riemannien associé %, soit stationnaire dans Q,, 

b) le groupe d’isométries laisse invariants H,,, 0, c, 1, x, €, Ml, 0. 
Nous allons montrer que s’il en est ainsi, le mouvement du fluide est permanent. 
Or les hypothéses a et b se traduisent en coordonnées adaptées par les conditions * 


Oo Gap = 0 Hap = OP = 09% = Oy6 =F 1l=HeE=Ay w= G0 =0. 


I] nous suffit de montrer que 09 G,.g = 09 u* = 0% g, = Oy Pb = % 6 =0. 


Soit x un point arbitraire de D,. Choisissons un systéme de coordonnées 
adapté (x®, x’) tel que le point % appartienne a la variété G d’équation x° =0. 
Supposons que la variété S soit orientée dans l’espace et qu’elle ne soit pas une 
variété exceptionnelle du probléme de Cauchy relatif aux équations du champ 
correspondant au fluide considéré. Nous savons alors (cf. §3) que pour un 


* Dans un mémoire antérieur «Sur une théorie relativiste des fluides thermo- 
dynamiques» [Ann. di Math. pura ed applicata, ser. IV 38 (1955)], nous avons étudié 
les mouvements permanents d’un fluide thermodynamique pur en supposant en plus 
des conditions 0 £48 = 6) = 0 l’hypothese 


l 
009° = CU? And — Vile a 
En fait, cette derniére hypothése est une conséquence de @)=0. Les conditions 


X £48 = X 0 = 0 suffisent pour assurer que le mouvement du fluide thermodynamique 
est permanent. 
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systéme de données de Cauchy © (8,8, % 83 Fag: 0, 6,9) portées par © et 
satisfaisant a 
(I) Or Oye Pe Oe 


le systéme des équations du champ admet une solution H (848, H5, 8, U%, p, 9) 
bien déterminée. De plus les équations (I) qui sont vérifiées sur S sont également 
vérifiées dans tout le domaine d’espace-temps 2, considéré. 


~; 


Considérons maintenant la variété ©’ d’équation 
bow) 


qui correspond point par point a la variété © dans l’homéomorphisme qui applique 
, sur le produit topologique D, x St dans lequel les lignes de temps z s’appliquent 
sur les lignes facteurs. En vertu de l’hypothése a) et b), la solution # (g,,, Hy 3, 
3, u*, p, 0) est telle qu’aux points z de mémes coordonnées locales (x') de S et 
G’, les quantités @(g,5, 0 Sup; Hug; 9, % 08) ont des valeurs égales. Pour la 
solution #, ces quantités vérifient des équations (I’) identiques a (I). Si donc 
on se pose le probleme de Cauchy avec les données @ précédentes portées par S’, 
on doit calculer d’abord les quantités u*, , 6, a partir des équations (I’), puis 
intégrer les équations du champ. Comme toutes les équations sont identiques et 
que les données sont identiques, on obtient une solution #’ identique a la solution 
H. Autrement dit, aux points de mémes coordonnées locales x’ les quantités 
(81, 11g, 0, u*, p, 6) ont mémes valeurs ainsi que leurs dérivées: elles sont donc 
invariantes le long des lignes de temps. 


I] en résulte que les hypothéses a) et b) entrainent 


Oy ut = 8p =%6=0 


et en vertu des équations de définition de G,, et ¢, 


Par conséquent Go Gap = 05.9, = 0= 


X B43 = X Hyg = X Gye X Up, Xo AGe Xp XO OF 

Le mouvement considéré du fluide est donc permanent. On peut énoncer le 
théoréme: 

Théoréme. Etant donné dans un domaine D4, un fluide parfait chargé et 
conducteur, pour que le mouvement de ce fluide soit permanent, il faut et il suffit que 

a) Vespace-temps riemannien associé &, soit stationnaire dans Da; 

b) le groupe d’isométries laisse invariants les champs ® et FH, et les coefficients 
C, 1, u, &, fh, 0. 


Remarque. — S’il existe pour le champ électromagnétique un potentiel 
vecteur global c’est-a-dire un champ de vecteur Q,, tel que 


Hy p= Pp — B Px 


ce qui est en particulier le cas si le domaine 2, est simplement connexe, on 
remplacera I’hypothése concernant le champ H,, par Vhypothése équivalente 


pan le potentiel vecteur: on supposera Y, invariant par le groupe d’iso- 
métries. 
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10. Isométrie induite dans B, 


Nous considérerons dans la suite un domaine simplement connexe D, occupé 
par un milieu chargé et conducteur. Nous supposerons plus généralement que 
le mouvement de ce milieu soit tel que l’espace-temps riemannien associé ¥, est 
stationnaire dans 2, et que le groupe d’isométries laisse invariants le vecteur 
vitesse unitaire % et V'indice n= eu du milieu 


(10.1) X 84g = XU = Xn=0. 


Il en est ainsi en particulier si le milieu considéré est un fluide parfait chargé 
en mouvement permanent. 


D’aprés (10.1), les quantités 
es 1 
Sap = Sup Fea (1 Yea al Uy Up 


sont invariantes par le groupe d’isométries de %,. I] en résulte que le champ 


s 
de vecteur contravariant € générateur du groupe d’isométries de B, détermine 
dans la variété riemannienne %, définie par la variété différentiable portant D, 
et munie de la métrique associée 


as? => Bp bx dx? 


un groupe connexe d’isométries globales ne laissant invariant aucun point du 
domaine , correspondant, et pour lequel le systéme de coordonnées (x®, x*) est 
un systéme de coordonnées locales adapté. On peut prendre pour générateur 


infinitésimal de ce groupe le vecteur © qui a pour composantes contravariantes 
Coes 10: 


Il est manifeste que le carré de ce vecteur a pour valeur 


(10.2) eon (1 va) (0) 


Introduisons maintenant la grandeur d’espace du vecteur % relativement a 
> 

la direction de temps €. Soit 

— w= gu'w’. 

A 5 5 > 

En vertu du caractére unitaire de “, on a 

at “)\2 os ey eee 

Bap Ui UP = —— (Bo _,U%)* + 8,0 ul = 1. 


800 


On en déduit 
(Uo)? = So0(4 + w?). 


En portant cette valeur dans (10.2) et en y remplagant — par V?, il vient 
(e a epee bop © = =) 


On voit que le signe de (¢)? peut changer. Dans &,, le vecteur & peut étre 
orienté dans le temps, dans l’espace ou étre isotrope. I] en est de méme des 
trajectoires d’isométries de ¥,. 
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IV. Etude géométrique des rayons électromagnétiques dans l’espace 
11. Un probléme du calcul des variations 


Nous nous proposons d’interpréter géométriquement les rayons ¢lectro- 
magnétiques dans l’espace a trois dimensions. A cet effet, nous commengons 
par rappeler briévement un probléme du calcul des variations. 

Etant donnée une variété différentiable &B,,.,, soit YW,,,,;3) espace fibré des 
vecteurs tangents aux différents points de ¥,,,. Si l’on adopte sur ¥,., des 
coordonnées locales (x%), chaque élément de Y®,,,,;4) sera constitué par la réunion 
des coordonnées (x%) du point x correspondant de &,,,, et des m +1 composantes 
(“6) du vecteur % dans le repére naturel en x associé aux (x%). Une structure 
de variété finslérienne sur ¥,,,, est définie par la donnée d’une fonction # (x, %) 
A valeurs scalaires dans YW,(,13) telle que pour *% fixe, Y(x, 1%) =AL(x,%). 
En coordonnées locales, une telle fonction est représentée par LY (x%, %°) et est 
homogéne et du premier degré par rapport aux #°. 

Considérons une variété différentiable ¥,,,, munie d’une structure de variété 
finslérienne et supposons qu’elle admette un groupe connexe a un parameétre 


oe 
d’isométries globales, de générateur infinitésimal €, ne laissant invariant aucun 


point de %,., (E+ 0). Supposons de plus que les trajectoires z du groupe sont 
homéomorphes a la droite réelle Ri, et soit B, la variété de B,., par la relation 
d’équivalence définie par le groupe. On sait qu’il existe des systémes de co- 
ordonnées locales (x°, x’) adaptés au groupe d’isométries, tels que dans le repére 


naturel associé, € ait pour composantes contravariantes 
cP ey = Ghee) 


et que les (x') soient un systéme de coordonnées locales arbitraire de &,. La 
donnée des (x*) détermine une trajectoire z. Pour déterminer un point sur cette 
trajectoire, on se donne la variété x° const A laquelle il appartient. 

Dans un systéme de coordonnées adapté (x°, x‘), ’hypothése d’isométrie se 
traduit par le fait que la fonction ¥ est localement indépendante de la variable x°: 


L = L (xi, ti, x9). 


Nous allons montrer qu'il est possible de douer la variété quotient &, de structure 
de variété finslérienne au moyen de fonctions L(z, 4) de facon qu’aux géodésiques 
de %,,,, extrémales de l’intégrale 


x. 


114A fr. , _ ae 
(11.1) i LAK) CU | eye aa 


correspondent par projection sur &, des extrémales de 


Z. 
i 2 dz 
Glo? 9 9 = 
( ) JE@,8) au, ar Prat 
Dans la suite, tout indice grec =-0,4,'2,...,.%; tout indice latin ==17 20 a7et 


nous supposons 


Inductions en relativité générale Ms 


Donnons nous une extrémale de (11.1) par une représentation paramétrique 
x*(u), w désignant un paramétre arbitraire. Le systéme différentiel aux extrémales 
aer(44'.1) 


ax 5 
stat ey 
(11.3) aa 
ou x” satisfait a 
(11.4) d 0% 2 _, 


CE TY 
est caractérisé par le fait d’admettre l’invariant intégral relatif 
(14.5) 0 = dt =a, Ldi+ a Lax. 


ie 


En vertu de l’hypothése 0, =0, on a l’intégrale premiére 
(11.6) Cals 
Comme 0),/==0, on peut résoudre (11.6) par rapport 4 %°; on a l’équation 
équivalente 
(11.7) O(n, 4, h) 
ot est une fonction homogéne et de degré 1 par rapport aux x’ et m dépend 
effectivement de h. 

Considérons la famille des extrémales (€,) correspondant a une valeur déter- 
minée de la constante #. Pour cette famille, le dernier terme de wm a la valeur 
hdx® et définit un invariant intégral relatif. Il en résulte que cette famille 
d’extrémales admet invariant intégral relatif 
(11.8) OS bx" 

Or, d’aprés l’homogénéité de #, on a 

Ki ORL 0, SF = 
Par suite, pour toute solution (11.6) ou (11.7), la quantité %* 0; Y peut s’exprimer 
par une fonction L des variables x kh 


(11.9) EG aah) = LK se) ola 2 h)| —h ola, %', h) 
oe aL =8:£ +0, L 0,9 —hd,9 =f. 
Ainsi, d’aprés (11.8), les projections des (€,) sur %, sont définies par un 
systéme différéntiel qui admet l’invariant intégral relatif 
as Od x", 
Autrement dit, elles sont extrémales de Vintégrale 
(11.10) f L(x*, x!) h) du 


%o 


ou fa la valeur choisie. 
On appelle descente la correspondance qui a la fonction # (x*, #4), 2): fait 
correspondre la fonction L(x*, x’, h). Le probléme inverse est possible*. 


* Voir LicHNEROwIcz, A.: Théories relativistes de la gravitation et de l’électro- 
magnétisme, Livre II, chap. I. 
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12. Projection des géodésiques de longueur nulle de la variété riemannienne B, 


Nous supposons que la variété riemannienne %, satisfasse aux hypothéses du 
paragraphe 10. La fonction @ est définie par la relation 


egy af hy he 


oti le second membre est une forme quadratique non dégénérée puisque g= 
det (g,) 0. Etudions d’abord les extrémales correspondant aux valeurs de %* 
pour lesquelles le second membre est positif. On sait d’ailleurs qu il suffit qu’une 
géodésique le rende positif en un point pour qu'il en soit de méme tout le long 
de la géodésique. 


Premier cas: %)) ne s’annule pas dans le domaine étudié. — Le procédé de 
descente nous conduit a former l’équation 
1 Ss i RY TR 
(42.2) £0, L* = Bug k  Poi% SNS 


et a éliminer x° entre cette Equation et 
(12.3) Eas Fee 


En décomposant ¥? en carrés a partir de la variable directrice %°, il vient 


1 1 Ci CO res 
Lr = — (5% £3) + ,, 444i 
oO 2D) 7 
oe £00 
ou 10n pose 
Tega ae &0i 807 
O19 Sy O47 = 
S00 


, 3 CS arn ‘, E aes Shs Pe 5 
et l'on voit que g;;%' x! est négative si {99 >0 et positive si y)7<0. Dans le premier 
cas on prendra h >max 99. Comme 406,4?=h¥F, on tire ]’équation 


TS te 
12.4 end BRS RE 
(12.4) ee - 


qui fournit Y en fonction des variables x*,%’, h. De (12.2) on tire ensuite 


(12.5) Perse nm Ee 
S00 800 
On en déduit d’aprés (12.3) et en vertu de (12.4) 
(12.6) L=el/(1 — 2 )8,, #4 +h Bit 
Pepuest 8o0 


ou ¢€ est le signe de Z,. 


i oa Eee une fonction de x*, x’, h homogéne et du premier degré par rapport 
aux *%/, FE ‘fini iété B iété fi 
ae ; e définit sur la vaniete quotient &; une structure de variété fins- 

enne. Inversement, ¢tant donnée localement dans &, la fonction L(x*, #’, h) 
precedente, on démontre facilement qu'il existe une fonction L(x", £/, 4°) 


homogéne et de degré 1 par rapport aux %%, qui par descente reconduit a L et 
que cette fonction est 


Les courbes extrémales correspondantes sont donc des géodésiques de &,. 
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Ainsi, les géodésiques de la variété riemannienne &%, qui correspondent a 


Vintégrale premiére 
Ch 


se projettent sur la variété quotient &, selon les extrémales de l’intégrale 


z 


(12.7) [ 
% 
ou a la méme valeur. Ces extrémales coincident avec celles de 


[en on E = 
is \/( — 5 \é,,#9 a pp Boi F, 
S00 ) §00 


(12.8) 


Le long de ces extrémales, on a d’aprés l’expression de %°: 
Bode 


h 4 ig ee 
(i) == SLB hee — 
a ot h? oy S00 
Loo 


Ceci étant, on peut définir les géodésiques de longueur nulle de 8, comme 
les courbes limites vers lesquelles tendent les géodésiques orientées dans le temps 


lorsque #->0. De la relation 
hf = Lox x 
Or 


il résulte que h—> co lorsque #0 et ha le signe de gy, % 


1 =e 2 oon 
aS (Boa %*)? + &:;% 4) = 0. 
00 


oa! 


On en déduit que 2),%* a une valeur non nulle et garde un signe constant. 
R, sur 


D’aprés (12.8), les projections des géodésiques de longueur nulle de 


¥, sont les extrémales de l’intégrale 
du. 


By de: 
s h2 oe on ee g 1 at 
hm Gs ye <5 =} i aS 
[|i h S00 ii f00 
0 
En passant a la limite, on en déduit le résultat suivant: 
Lemme 1. Les géodésiques de longueur nulle de B, se projettent sur B, selon 'es 


extrémales de l’intégrale 
ay 
/ A aoe a. xt 
i [« é' / eee i ie oe) du 
§00 S00 


4 


(12.10) 


ou € est le signe de Bo, et &' le signe de 8), x" 
D’aprés (12.9), le long de ces extrémales on a 
Boi ax4 


dx = ee’ = dk dxi — 
oor Soo 


(12.11) 


On remarquera que dx°=Ldu. 
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Second cas: 2))9=0. — On a alors 
(12.142) G2 = 275,044 8,,0 0. 


Nous supposerons &,;*’ fini ==0. Le procédé de descente nous conduit a éliminer 
Fete entre (12.12)-et 


(12.13) Bit =hL, 
(12.14) La Gh. 
De (12.13) on tire 


On en déduit 


L’équation (12.14) détermine alors L 


meet gg 
Lice B08* py, SU 
(12.15) Tero 


Inversement, a toute fonction L de la forme précédente correspond par 
montée la fonction ¥ définie par (12.12). On note que L se présente par rapport 
aux variables %* comme le quotient d’une forme quadratique par une forme 
linéaire. 

Ainsi, dans le cas 9,0, les projections des géodésiques de &, sur ¥, sont 
les courbes extrémales de 


pour la valeur correspondante de la constante h. Ces extrémales coincident avec 
celles de 


Comme précédemment, les projections des géodésiques de longueur nulle sont 
définies par 
a 
ie me 
i lim ( elias, )au 
hoo 2h? 284%" 


& 


On en déduit en passant a la limite, le lemme 
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Lemme 2. Dans tout domaine de WB, 0% Boo =0, les géodésiques de longueur 
nulle de 8,4 se projettent sur %B, selon les extrémales de V’intégrale 


ay 


TE 
(12.16) i ee ay 
é 280i %" 
% 
Le long de ces extrémales on a 
dx = a 
20504° 


13. Le principe de FERMAT 


Nous avons établi que les rayons électromagnétiques sont géodésiques de 
longueur nulle de la variété riemannienne %,. Nous pouvons les interpréter 
géométriquement dans l’espace si le milieu considéré est en mouvement per- 
manent. En effet les lemmes1 et 2 fournissent une démonstration immédiate 
du théoréme suivant: 


Théoréme. Sz le mouvement du milieu considéré est permanent et tel que 
Boo foo’ a Vw) = 0, 


les rayons électromagnétiques dans l’espace sont des lignes réalisant Vextrémum de 
Vintégrale 
; fe oe eae as ee Zoi x 
13 A ieee Fee ap ant AMEN HY, 
( 3 ) / 80 Si §o0 


00 


0) 


pour des variations a extrémités fixes, ou € est le signe de So et &' le signe de By, x* 
Le temps mis par un rayon pour aller du point %9 au point z, est donné par 


(13.2) fae=f pees g,,46%) — ERE Ey 
é S00 S00 


Ce temps est extrémum. 


Dans le cas ot £9) =0, on obtient un énoncé analogue en remplacant (13.1) 
et (13.2) respectivement par 

? Bish! a 
(13.3) eee du 


Zo 


et 


Ny Bit 
(13.4) fe / 281% ae 


Par le théoréme précédent se trouve démontrée |’équivalence du principe 
géodésique et du principe du moindre temps. 

En particulier, si l’univers est statigue au sens de LEvi-Civita, c’est-a-dire 
si les lignes de courant coincident avec les lignes de temps, l’espace-temps %, 
est orthogonal. Soit 


ds®* = U(dx)? + g,,dxi dx 
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la métrique d’univers de &,. Les ui étant nuls, on en déduit la métrique associée 
Do Ald oe 
ds? = 5 (a0) + ganda. 


On peut mettre (13.2) sous Ja forme 


oe tee 


ou l’on a posé 2 
do® = — g,,dx' dx. 

Cette intégrale rappelle une intégrale qui se présente dans l'étude des fluides 
holonomes. Or U est le potentiel principal de gravitation. Le champ gravi- 
tationnel agirait si l’on ose dire, comme une espéce de pression sur les rayons 
électromagnétiques. 

Si U=1, on démontre que l’espace-temps ¥, est euclidien. L’énoncé du 
théoréme devient 

6fdx®=06fndo=0. 

Nous retrouvons |’énoncé exact du principe de FERMAT en optique. Le théoréme 
que nous avons établi, en constitue donc l’énoncé généralisé en Relativité générale. 


14. Interprétation du signe de gp, &* 
Nous avons interprété %° comme la variable temps. L’interprétation du signe 
é’ de £,,%* est simple. En effet, l’équation 
(C,) Lau? = = (oud)? + 8,,dx'dxi =0 
00 


représente le c6ne caractéristique ©, au point x des équations de MAXWELL. Les 
deux nappes de ce cone sont symétriques par rapport 4 l’hyperplan élémentaire 


(t,) Page S05 


_Désignons par M (x*) le sommet de ce céne ©,. Prenons un couple de points 
voisins de M, ayant pour coordonnées spatiales (x'+dx') et appartenant res- 
pectivement aux deux nappes du céne et symétriques par rapport a z,. Soient 


My (xi +dx, ®4+dx), Mi(ei+ dx, 2° —d'x). 


On peut dire que MM, représente aux infiniment petits d’ordre supérieur prés, 
le deplacement infinitésimal associé 4 un rayon électromagnétique allant du 
point d’espace A (x’) au point d’espace A’(x’-+-d x‘) dans le temps dx°. De méme, 
M My, peut étre considéré comme représentant le déplacement infinitésimal associé 
a un rayon électromagnétique allant du point A’(x!+dx') au point A(x‘) dans 
le temps: d’x°. 

Les deux points M, et Mj sont symétriques par rapport a l’hyperplan élé- 
mentaire z,, on doit avoir 


Cond ee a Cond x. 
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On en déduit é 
dx = dx +2 Guide 

800 
Cette relation montre que, sauf dans le cas statique, le temps mis par un rayon 
pour aller du point d’espace A(x’) au point d’espace A’(x’-+dx') n’est pas le 
méme que le temps mis par un autre rayon pour aller de A’(x‘+dx‘) a A(x’). 


15. Espace-temps de MINKOWSKI et loi relativiste de la composition des vitesses 


Placons-nous dans le cas d’un espace-temps sans gravitation de MINKOWSKI, 
rapporté a un systéme de coordonnées galiléennes réduites. Nous avons la 
métrique d’univers 


(15.1) as? == (dx°)? — (da)? fd. x7)? —-(d x82: 


u représente dans le cas actuel le vecteur vitesse unitaire d’univers dont les 
composantes sont déterminées classiquement a partir de la vitesse d’espace B si 
la vitesse limite c est prise comme unité. 
Un calcul facile donne la métrique associée 

(15.2) as = (@xrt2t = Bd di — > (dzi)!— (B.d.x!)2 
Cette métrique est du type hyperbolique normal comme la métrique d’univers 
(15.1). Il est cependant intéressant de noter un changement de l’ordre dans la 
signature de cette métrique au passage de V?—f?. On le met facilement en 
évidence en choisissant l’axe des x1 paralléle a la vitesse B du milieu. On a ainsi 
la métrique 


(15.3) ia (dx9)24 2 Bax dx 1a aaa 


que l’on peut mettre sous la forme canonique par une décomposition en carrés. 
Si V2 6? on obtient 


a 


ds2 = sae peat yay esihere VE op pale a8) VE (d x1)? — (d.x?)? — (d x8)? 


2 p2 | 4 — p? ! 1— p? V2— p 
et l’on voit que pour V?> f? cette métrique a pour signature «+ — — —» et 
pour V2<f? elle a pour signature «—-++ ——». Pour V?=f?, on obtient 


d5*=2V dx de — (4 + V2) (dx)? — (d x2)? — (dx). 


On vérifie également que cette métrique a pour signature «+ ———» en la 
mettant sous la forme 
ds2 Se (axe)2 eee [4 + V2) dxt+ Vdx?]? — (dx?)? — (dx?) 


A partir de la métrique associée (15.2), cherchons a exprimer le théoreme de 
FEerRmatT en prenant l’arc o du rayon électromagnétique comme parametre. Nous 
avons a remplacer dans (13.2) ** par 


es 
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“ " +/ ayy: ae: 7 : eer ; ey : 
(15.4) hth e | te OS [v2 — p+ (1 — V%) (6; a')?] (1 — ue bdo 


et l’on peut en déduire 


ie 4. LU 2S ee 
oF ae ee pay (+ VBA — 
Si V2—f? +0, cette relation donne 

(15.5) 1—8— (= pyWhe (VW eee 


=> bad i < 
Si on interpréte V comme vitesse absolue et W comme vitesse relative de 
propagation de l’onde électromagnétique considérée, on a manifestement 


(145.6) Ve=——* [w+ p2+ 2W-B + (W-B)?—W?B?]. 

ce LW By 
On veérifie (bed oe calcul direct 4 partir de (13.4) que cette relation reste valable 
dans le cas ob V2 =f?. C’est la formule relativiste de la composition des vitesses. 
Il est aisé de vérifier qu’on peut la mettre sous la forme* 


1? S 


rasa Bll =O 


Nous obtenons ainsi a4 partir du principe de FERMAT une démonstration de la 
loi relativiste de la composition des vitesses. 
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Erzwungene Schwingungen in rotierenden Fhissigkeiten 


HANSJORG OSER 


Vorgelegt von H. GORTLER 


§ 1. Einleitung 

In einer frither erschienenen Arbeit zeigt GORTLER [1], daB beim Problem 
kleiner Schwingungen in einer gleichférmig rotierenden, inkompressiblen und 
reibungsfreien Fliissigkeit sich der Typus der beschreibenden Differentialglei- 
chungen andert beim Durchschreiten einer gewissen kritischen Frequenz der 
Schwingung. Bis heute fehlte aber der experimentelle Nachweis der im hyper- 
bolischen Fall auftretenden Charakteristiken, die, von den Unstetigkeitsstellen 
ausgehend, das ganze Strémungsbild aufteilen. Dem Verfasser ist es in seiner 
Dissertation* gelungen, diesen Nachweis zu fiihren und die Charakteristiken 
fotografisch festzuhalten. AuBerdem konnte ein der experimentellen Anordnung 
entsprechendes mathematisches Randwertproblem formuliert und exakt geldst 
werden. Wéahrend der physikalische Teil Inhalt einer anderen Veréffentlichung 
werden soll, wird die mathematische Behandlung des Problems hier in gekiirzter 
Form dargestellt werden. 


§ 2. Die Bewegungsgleichungen 
Die Eulerschen Gleichungen fiir die Bewegungen einer reibungsfreien, in- 
kompressiblen Fliissigkeit lauten auf ein mit konstanter Winkelgeschwindigkeit 
@ rotierendes Koordinatensystem bezogen: ‘ 
a +2:oxXv= ~ 7 Brad p. 
Hinzu tritt die Kontinuitatsgleichung: 
divp = 0. 


Hierbei ist » die Geschwindigkeit, gemessen im rotierenden System, 9 = const 
die Dichte der Fliissigkeit und 
p= p* — 20°70 + O82. 

p* ist der statische Druck. Es gilt nun unabhangig vom Koordinatensystem die 
Vektorbeziehung: de 

a 
Es seien u,v, w die Geschwindigkeiten in den Richtungen 7’, g’, bzw. z’ eines 
zylindrischen Polarkoordinatensystems. Setzt man die auftretenden Relativ- 


ov p2 
= — — — t-Vr 
Aa ered = vp X rot bv 


* Universitat Freiburg/Breisgau, 1957 
Arch, Rational Mech. Anal., Vol. 1 6 
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geschwindigkeiten als so klein voraus, daB die quadratischen Glieder vernach- 
lassigbar sind, so erhalten wir, wenn der Rotationsvektor mit der z’-Achse zu- 
sammenfallt, bei Rotationssymmetrie: 


du 1 op 
——— U = — - Tay 
at 20 0 67 
ov 
2mu + — = (0, 
ou + ay 
(1) Cr eileen Re 
at (the Cais 
1 ou ow 
—u + —-l == (OM 
y’ Oo” 1 Oz’ 


Wir legen nun das folgende physikalische Problem zugrunde: In der Rota- 
tionsachse der unendlich ausgedehnten, gleichformig um die z’-Achse rotierenden 
Fliissigkeit schwingt harmonisch eine Kreisplatte vom Radius 7, auf und ab. 
Die Ebene der Platte soll dabei stets senkrecht zur Rotationsachse stehen, die 
Ruhelage der Platte sei die Ebene 2’=0. Beschreiben wir die Plattenschwingung 
Rey z’= — bcos Pt 
und setzen a=fb, dann wird der folgende Ansatz fiir die Geschwindigkeits- 
komponente w/(r’, 2’, t) nahegelegt, um die Zeit zu eliminieren: 


wr’, 2,1) =4a4Wwy 2 psunpas 


Setzt man dies in die Gleichungen (1) ein, so ergeben sich folgende Phasen- 
beziehungen zu den anderen Funktionen (m ist ein MaBstabfaktor, der noch 
von fB und w abhangt): 


Ut Rote = U(r’, 2’) sin Bt, 


OM ase set) We sed Vr, zy eos Bt, 


m 2@ 
p(r’, 2',t) =aBomr, P(r’, 2') cosBt. 


Die bei den Funktionen U,V, W und P stehenden Faktoren sind so gewahlt, 
da8B diese Funktionen dimensionslos werden und unter ihnen besonders einfache 
Relationen entstehen. Wir bezeichnen U, V und W im folgenden kurz als Ge- 
schwindigkeiten und P als Druck. 

Fihrt man die Ansatze fiir u,v, w und # in die Differentialgleichungen (1) 
ein und eliminiert die verschiedenen Funktionen der Reihe nach, so entstehen 
zwei verschiedene Differentialgleichungen zweiter Ordnung. Davon gilt die eine 
fiir W und P: a2 42 eae, uae 

| ay a y’ or’ : a. ma (W, P) =0. 


Die andete gilt fiir U und V: 


(Cha 4 @ { B?—4@2 @ 
ler Ate Y or ye wee Be sah (U,V) =0. 


Beide Differentialgleichungen haben die Eigenschaft, elliptisch oder hyperbolisch 
zu sein, je nachdem ob 
B22. 
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Unter Beschrankung auf den hyperbolischen Fall setzen wir 


4w2 — f2 
oa = tee, 


m* = 


wobei w der Winkel ist, den die Charakteristikenkegel mit der Ebene z’=0 
bilden. Zur bequemeren Rechnung fiihren wir dimensionslose Koordinaten ein, 
indem wir fiir 0<B<2w setzen: 


B= 2-%y-™M 


mit der nicht wesentlichen Einschrankung, daB f von 0 und 2w verschieden ist *. 
In den neuen Koordinaten lauten die Differentialgleichungen: 


(ha i —@) oe? 
(ia) ee oh mone Reis 
oC? i) 1 o? i 


Ferner wird, wie man sich durch Einsetzen in die Differentialgleichungen (1) 
tiberzeugt: 


(2) Ulr,2) =—— [ (W).dr, 
(3) V(r,z) = Ulr,2), 
(4) P(r,z) =— f Waz. 


Randbedingungen sind zu stellen an der Oberflache der schwingenden Platte. 
Im Rahmen unserer linearisierten Theorie begehen wir keinen Fehler, wenn wir 
statt dessen die Randbedingungen in der Ebene z=0 stellen, da wir kleine 
Schwingungsamplituden voraussetzen. 

Die Randbedingungen lauten: 


dala 
a) we 0) =| 

OL en S41, 
b) U(r, 0) = V(r, 0) = 0. 


U,V, W und P sollen verschwinden fiir y—> co und z+ ~. 
Die Anfangsbedingung soll lauten: 
w(r,z,0) =0 
(ist bereits durch den Ansatz erfiillt). 
Die fehlende vierte Randbedingung (fiir den Druck P) wird ersetzt durch die 


Symmetriebedingung : Wien eae ne ney. 
Hierzu aquivalent ist die Bedingung, da8 der Druck in z=0 verschwindet. 
Die Lésungen sollen iiberall eindeutig und regular sein, ausgenommen auf 
Charakteristiken. 


* Fiir B=0 liegt iiberhaupt keine St6rung mehr vor, wahrend fir B=2w die 
Loésungen sich so sehr vereinfachen (parabolischer Grenzfall), daB keine physikalisch 
sinnvollen Randbedingungen mehr erfiillt werden kénnen. Dort miissen die nicht- 
linearen Gleichungen untersucht werden. 
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§ 3. Losung des Randwertproblems 


Da die Randbedingungen auf Koordinatenflachen vorliegen, wird ein Sepa- 
rationsansatz fiir die Funktionen U, V, W und P nahegelegt. SchlieBen wir die 
in der Rotationsachse singularen Lésungsbestandteile von vornherein aus, so 
erhalten wir fiir W die partikulare Lésung 


T,(k1) (A, cos kz + B,sin kz) 
und fiir U und V ganz entsprechend: 
J, (kr) (A,cos kz + Bysin kz). 


Hierbei sind J, und J, die Besselschen Funktionen erster Art, A; und B; sind 
Koeffizienten, die noch von der Separationskonstanten k abhangen kénnen. 
Die Symmetriebedingung liefert sofort B,=0, das partikulare Integral fiir W 


lautet also: T,(k1) Ay (R) cos kz. 


Wir integrieren nun iiber k von 0 bis oo und verfiigen tiber 4,() so, daB die 
Randbedingung fiir W(r, 0) erfiillt wird. Es soll also sein: 


a On tei al 
W(r,0) =f Fo(hr) Ay(®) a =| 
0 4 Teal 
Nun gilt aber die Beziehung*: 
is O 4 
I Tolk”) i) ak =| 
0 4 are We 


Die Funktion J,(k) ist die einzige Lésung fiir A,(k), denn wenn es eine zweite 
gabe, miiBte man mit der Differenz die Randbedingung W(r, 0) =0 fiir alle 7 
erfiillen kénnen. Aus der Fourier-Besselschen Darstellung einer beliebigen Funk- 
tion ** folgt aber dann, daB diese Differenz null sein muB, also wirklich nur eine 
Lésung existieren kann. Damit wird also: 


(5) Wir, 2) = in: (kr) J, (R) cos (kz) dk. 


Durch Anwendung der Beziehungen (2) bis (4) ergibt sich daraus fiir die iibrigen 
Funktionen: 


(6) OG aie in; (kr) J, (Rk) sin (kz) dk, 
(7) Vir.2) =f (br) E.) sin (BD) dB, 
(8) P(r,2) = — [ ehler) J, (A) sin (kz) dk. 


Durch die angegebenen Ausdriicke werden samtliche Randbedingungen im End- 


lichen erfiillt. Da8 die Funktionen im Unendlichen alle verschwinden, wird in 
§ 4 gezeigt werden. 


* Watson [2], S. 406 
** SOMMERFELD [3], S. 112 
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1. Die Funktionen U und V. Die Fourier-Integrale (6) und (7), welche 
U und V darstellen, kénnen geschlossen ausgefiihrt und durch bekannte Funk- 
tionen dargestellt werden. Durch Spezialisierung einer allgemeineren Trans- 
formation * ergibt sich folgendes: 


0 Oar aa i 
1 VA = 72 Z 
U=V= at ae Vai) ae a ey 
et A 
= Se oa A+7<z< ©. 


P, und Q, sind gewohnliche Kugelfunktionen erster und zweiter Art, was man 
auch bestatigen kann, wenn man mit x= (7?+1—2z?)/27 in die Differential- 
gleichung (1a) eingeht und dann die 


Differentialgleichung der Kugelfunktionen 
vom Index } erhalt. Die Darstellungen 
fir U und V gelten zunachst nur fiir 7-0 
und ry+1. Die Grenzwerte fiir r=1—0 
und y=1-+ 0 sind aber gleich und ebenso 
die Ableitungen an dieser Stelle, so daB 
die Funktionen dort sogar analytisch sind. 
Fiir y= 0 ergibt sich aus der Darstellung (6) 
sofort : 
C052) = V(0;2)==0 

mit Ausnahme von z=1, wo das Integral 
nicht konvergiert. 


2. Die Funktion W. Fiir J,(k7) J, (-) He ue ; 
P " 4 " 5 - Abb. 1. Charakteristiken bei einer rotationssymmetri- 
ist die Fouriertransformierte leider nicht schen ebenen Stérung im Schnitt 


in geschlossener Form bekannt. Deshalb 

mu auf eine Darstellung durch elementare Funktionen verzichtet werden, wie 
es fiir U und V méglich war. Fiir gewisse Bereiche der Variablen kénnen aber 
auch bei der Funktion W Aussagen gemacht werden. Es ist fiir 7 =0** 


of 4 One —e4 
W (0, 2) mul) cos (k 2) i — | ES = — 1) t| <P E SS). 


Geeignete Spezialisierung zweier allgemeiner Beziehungen *** ergibt ferner aus 


dem Integral (5): Wi a eel, at eae? A ey <4) 


Wis) = 0) peallte a0 aren ih (Xe A) 
Damit sind bereits zwei wesentliche Feststellungen mdglich: 

a) Der kegelférmige Bereich IV iiber der Platte (s. Abb. 1) besitzt nur eine 
Geschwindigkeitskomponente in vertikaler Richtung. Die Fliissigkeit schwingt 
hier genau (d.h. nach Amplitude und Phase) mit der Platte. 

b) Der Bereich I ist véllig stérungsfrei. Hier ist U=V=W=0. 


us OBERHETTINGER [4], S. 169 
** ERDELYI [5] Gl. 1.12.3 
xxx ERDELYI [5] Gln. 1.12.26 und 1.12.25 
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In den Bereichen II und III ist zuniichst nichts zu erfahren tiber die Funktion W. 
Der Weg zu ihrer numerischen Berechnung wird in §5 angegeben. 


3. Der Druck P. Auch hier sind allgemeine Aussagen nur in den Bereichen I 
und IV und auf der Rotationsachse méglich. Die endgiiltige Berechnung von 
P haben wir wie bei W daher nach §5 verschoben. Es ist*: 

a) fiir 7=0: 

(<2 C24; 

P(Og2\ i= ES 
|—-@—-j2—1)  1422< 00. 
b) im Bereich IV iiber der Platte: 
Play SP ORS ETS TF PALS 
c) im stérungsfreien Bereich I: 
P(r,z) =0 Very WN ia 


Bei diesen Fallunterscheidungen treten die Charakteristiken z=1—7, z=r—1 
und z=1-+,7 jeweils als Grenzen der Giiltigkeitsbereiche auf. Man erfahrt daher 
nichts iiber die Funktionen auf den Charakteristiken selbst. Wir werden aber 
in §6 die Umgebung der Charakteristiken genauer untersuchen. 


§ 4. Die Erfiillung der Randbedingungen im Unendlichen 


Es fehlt bisher noch der Nachweis, daB die angegebenen Lésungen im Un- 
endlichen verschwinden, was man zwar aus dem bisher Gewonnenen schon ver- 
muten kann, was aber noch keineswegs bewiesen ist. Fiir 7— oo ertibrigt sich 
ein Beweis: Alle Funktionen verschwinden fiir y>z-+1 identisch. Um das Ver- 
schwinden der Funktionen U und V fiir z— oo zu zeigen, kénnen wir von ihrer 
Darstellung als Kugelfunktionen ausgehen. Es ist fiir z>7+1: 


SUT Nc ee ge — 727 —4 
Ut Ce eee Ol eae! 
Nun gilt fiir Q,(x) in der Umgebung von x= 
cy 1 
1 (XG fiers meer 
21) 4/2 # 


Fiir z— co ergibt sich also: 
U(r, 2) = V(r,2)  — 7/23. 


Um das Verschwinden der Funktion W(y, 2) fiir z—> co zu zeigen, miissen wir 
auf die Integraldarstellung (5) zuriickgreifen. Wir werden gleichzeitig die Existenz 


des Integrales (5) mit aufzeigen. Zu diesem Zweck spalten wir (5) in folgender 
Weise auf: 


(*) Wie ae ey ae 


m 


Die Aufspaltung erfolgt so, daB beim Eintragen des ersten Gliedes der asymptoti- 
schen Entwicklung der Besselfunktionen der Fehler fiir Rm unter einer festen 
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Schranke bleibt. Das erste Integral wird durch partielle Integration umgeformt: 
1 7 ; |k=m 4 r d : 
(#x) [ dk =~ folky) (sinks % f ap Volk”) A(t) sin ke dk. 
0 0 
Das letzte Integral existiert und hat einen endlichen Wert. Es ist namlich 
ae 1 
ap olkr) L(®)) = — 7 Klkr) hlk) + Joh) lk) — = John) A) 
eine im Intervall O< k< m stetige (und damit beschrankte) Funktion*. Daher 
sind bei beliebigem endlichem m alle Ausdriicke auf der rechten Seite der 
Gleichung (**) beschrankt. Fiir 7 co geht also die rechte Seite gegen Null. Es 
bleibt weiter zu zeigen, da8 auch das zweite Integral von Gleichung (*) ver- 
schwindet fiir zoo. Man tragt das erste Glied der asymptotischen Ent- 
wicklungen fiir J,)(k7) und J,(k) in dieses Integral ein und erhalt nach einigen 
trigonometrischen Umformungen (fiir z>7 +1): 


fi pe 2 (es m/4) cos(k — 32/4) coskz de 
GP k 


™m. mm 


: -{Ci[m(r +2+1)] + Ci[m(y —2z—1)]+ 


~ 227 
+ stl[m(r +2 —1)] — st[m(z —r + 1) ]}. 
Bei beliebigem, aber festem m>0 gehen die Funktionen Ci(mz) und _ s7 (mz) 
gegen Null fiir z— oo **. 

Es ist damit nicht nur die Erfillung der Randbedingung im Unendlichen 
gezeigt, sondern zugleich die Existenz des Integrals bewiesen, weil fiir festes 
yund z (==7 +1) das Restglied beliebig klein gemacht werden kann. Der Beweis 
fiir die anderen Integrale (6) bis (8) lauft genau entsprechend. Im Zwischen- 
gebiet |y—1]<z<r-+1 tritt nur eine unwesentliche Vorzeichenanderung bei 
den Argumenten der Funktionen Cz(x) und si(«) auf. Man sieht weiter, daB 
fiir z=r-+1 die Funktion W unendlich werden muB, da Cz (0) = — ov ist. 

Um das Verschwinden der Funktion P(r, z) fiir zoo zu zeigen, kénnen wir 
einen Satz aus der Theorie der Fourierschen Integrale benutzen. Er lautet ***: 


Ist 


+00 
J \f(@| dx <0, 
dann konvergieren die Integrale 
+00 +00 
{t(x)coszxdx und ff(x)sinzxdx 


gegen Null fiir z— oo. 
Die Voraussetzung ist fiir unsere Funktion P(r, z) erfiillt. Wir zeigen dies, 
indem wir das betreffende Integral (8) wieder aufspalten: 


oo 4 E co 1 ; 
fl-ldes | jet FACIL +f Io(hr) J, (R)| ah. 


* Die Entwicklung von /,(k) beginnt mit k/2 
** Man vergleiche JAHNKE-EMDE [6], S. 3 
xxkk TITCHMARSH [7], S. 11 
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Das erste Integral existiert wegen der Stetigkeit des Integranden und hat einen 
endlichen Wert. Das zweite Integral konvergiert, weil der Integrand tberall 
beschrankt ist und fiir k-> co wie 1/k? abnimmt. Um in vdlligen Einklang mit 
der Voraussetzung zu kommen, miissen wir /(x) =0 setzen fir <0. Dann folgt 
die Behauptung, daB P(r, z) verschwindet fiir z— ov. 


§ 5. Numerische Berechnung der Funktionen 


Die Funktionen U und V sind durch die halbzahligen Kugelfunktionen 
P,(cos #) und Q;(*) (*>1) ausgedriickt. Nun stehen fiir diese leider nur 
Tafeln fiir x>1 zur Verfiigung, so daB die Darstellung durch die Funktion 
P,(cos 8) nur einen formalen Gewinn bedeutet. Wir kénnen diese Funktion 
aber durch elliptische Integrale darstellen: 


(0) P,(x) = P, (cos #) = * [2E(#) — K(#)]. 


Auch fiir die im folgenden noch auftretende Funktion P_, (cos #) ist eine solche 
Darstellung méglich: 


(10) P_,(x) = P_,(cosd) = * K(h). 


1 
2 


Dabei sind K(k?) und E(k?) die sog. vollstandigen elliptischen Integrale erster 


und zweiter Gattung und ats ore pps 


2 Zz 


Die angegebenen Beziehungen konnen sehr leicht bestatigt werden, wenn man 
die in Frage stehenden Funktionen durch die hypergeometrische Funktion dar- 
stellt und Potenzreihenvergleich anstellt. Es ist mit F=,F, 


P, (a) =F(—mjn+134;45*), 


ia 


(11) ed alirec ( 


2 


la on 


4 
A a is Tle ie oe, 


Damit sind U(7,z) und V(r, z) tatsachlich auf Funktionen zuriickgefiihrt, die 
vertafelt vorliegen. 


Die Berechnung von W(r, z) und P(r, z) geht nun zweckmiaBigerweise nicht 
von den Integralen (5) und (8) aus, da die numerische Auswertung dieser un- 
endlichen Integrale in Abhiangigkeit von zwei Parametern einen betrachtlichen 
Aufwand erfordern wiirde. Wir fiihren diese Funktionen vielmehr mittels der 
Beziehungen (2) und (4) auf U(y, z) zuriick. Die Differentiation in (2) kann auf 
die Kugelfunktion ausgeiibt und das Ergebnis daher wieder analytisch dargestellt 
werden. Es bleibt somit nur eine numerische Integration fiir W und eine weitere 
fir P zu leisten tibrig. Da numerische Integrationsverfahren im allgemeinen 
gegentiber Rundungsfehlern recht unempfindlich sind, kann mit einem maBigen 
Stellenaufwand gerechnet werden. Es ist nach (2) 


3 


Pree Wk! 
W= [ 0), dz. 


« 
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Die untere Grenze bleibt unbestimmt. Setzen wir sie gleich y—1, dann miissen 
wir eine Funktion /(7) als Integrationskonstante hinzufiigen. Wir haben sie 
durch Auswerten des Integrals (5) fiir verschiedene Werte von 7 bestimmt. In 
§ 6 werden wir durch eine OSTA RSUMMETAS, I As z)=f(r) exakt bestim- 


men. Setzen wir U in (2) ein und beachten die Differentiationsregel fiir die Kugel- 
funktionen, so erhalten wir fiir 
a) |r—1|<z<r4+1 


& 


is 1 eat 2) A ee) Aha 
(12) W = Aryy i = eee ea az . 
Fithren wir gemaB (11) die elliptischen Integrale ein und setzen auBerdem noch 
KE 
7 D 


(auch diese Funktion ist in den Tafeln der elliptischen Normalintegrale ent- 
halten), so ergibt sich schlieBlich: 


z 


x 1 2E E Lee, moles eae 
(13) W ae {is L(y 1)|D ! alee te ee 


P(r, z) erhalten wir durch nochmalige Integration von —W nach z. 


gale ye 


b) 1447<z%< co, Mit y= a —* erhalten wir ganz entsprechend: 


1 & +9) O49) — (1 +19) Qy(9) 
(14) Wie ea ae BT are — dé. 
Fiir die numerische Auswertung giinstiger ist die Darstellung: 
Lee Oo ee 
(15) We sr fe y—1 T (7 4) — ya — 1 ape 


P(r, z) wird wieder nach Gleichung (4) berechnet. 

Bei der numerischen Berechnung nach (13), (14) und (15) wurde die untere 
Grenze der Integrale jeweils y—1 gesetzt. Um die Integrationskonstante (die 
eine Funktion von 7 ist) zu bestimmen, mu8ten wir die Integrale 


(5) W= f Alkr) J (A) cos(kz) dk, 
0 
(8) P=— [ 5 Jolhr) h(t) sin (he) ak 
0 
fiir y=1, 2,..., 10 fiir je ein z<7v+1 und z>7+1 numerisch auswerten. Dabei 


wurde folgendermaBen vorgegangen: Die Integranden wurden vertafelt im 
a k: 0 (0.1) 10 bei W 
k: 0 (0.1) 7 bei P (bessere Konvergenz) 


und nach der Simpsonschen Regel integriert. Der Rest wurde durch Ein- 
tragen der asymptotischen Entwicklungen der Besselschen Funktionen in den 


90 HANSJORG OSER: 
betreffenden Integralen abgeschatzt zu: 


Ry (m) = ae {Colm (2 + P= 1) + Co (m(r + 1—2)) + 
+ si(m(r —1+2)) —si(m(z—7+1 yes 


(@+r+4) 
Rp(m) = Sa {(e +7+4)|Ci(m(z+7r-+1)) ae : ) 
j (y+1—2z))] 
— (ry +1—2) |Ci(m +1 —2)) “oe — ; ) | 
== il 
Sera) Wee: Mpa eee et 


— (=r +4) [si(m@—r +0) 4 omer “I 

Fiir die Funktion W war es méglich, auf Grund der in §6 angestellten Unter- 
suchungen die Funktion /(7), welche die Rolle der Integrationskonstanten ver- 
tritt, analytisch herzuleiten. Die Ubereinstimmung mit der numerisch gewonnenen 
Funktion war bis auf wenige Promille genau. 

Eine Schwierigkeit bei der numerischen Behandlung von (13) bis (15) ist 
noch zu erwahnen: Die Funktionen U, V, W werden alle unendlich bel z=r7-+1. 
Beim Integrieren iiber diese Stelle hinweg mu8 daher verlangt werden, daB ent- 
weder die Singularitat so schwach ist, daB sie integrabel ist, oder aber die beiden 
unendlichen Betrage miissen sich gegenseitig wegheben, wenn man dariiber 
hinwegintegriert. Es liegt hier der letztere Fall vor. Da wir, bedingt durch die 
Vertafelung, nicht von beiden Seiten gleich nahe an die singulare Stelle heran- 
kommen, miissen wir diesen Beitrag beriicksichtigen. Zur numerischen Erfassung 
dieser Unsymmetrie gehen wir zweckmabigerweise von (12) und (14) aus, da 
die Handhabung der Kugelfunktionen am leichtesten ist. 


1 (lr) P+ r—2)Pyl@) : 
Wel Fee i ts 2 ee Ke = Fe 
1 pf &—*#)Q_\(—*% —U =x) 0,(—2) 
ea ap 


Hier haben wir —x statt y gesetzt. Fiir z>7+1 geht x+—1. Wir werden 
demnach P, und P_; durch die asymptotischen Entwicklungen um —1, die 
Funktionen Q, und Q- , durch die Entwicklungen um +1 ersetzen und an- 
schlieBend den Cauchyschen Hauptwert des Integrals zwischen 1+7—e, und 
1-++r7-+ e, bestimmen. Wir erhalten so die KorrekturgréBen : 

Ay =- Z = log &2/&,, 


27 


| ee ee : 
Ap eal {(e aE £5) — &, log &, — & log €y}. 


Mit Hilfe dieser Korrekturen war es dann méglich, W und P bis etwa z=12 
(Tafelgrenzen) zu berechnen. 


§ 6. Die Umgebung der Charakteristiken 


Von besonderem Interesse ist das Verhalten der Lésungen in der Nachbar- 
schaft der Charakteristiken. Leider geben die numerischen Lésungen dort keine 
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feineren Einzelheiten wieder, so daB8 wir analytische Untersuchungsmethoden 
heranziehen miissen. 


Setzen wir in die Integrale (6) und (5) die folgenden Beziehungen ein*: 


J,(e7) J, (k) =a J hee cos pag, 


Soar) hk) = = fat (1 —rcos y) dp 


mit 
ae ere, Lae 


Wegen der Existenz der Integrale (6) und (5) kénnen wir die Integrationsfolge 
nach g und k vertauschen und erhalten: 


(16) Ul ~ {cosol fia (kw) Jsinkadk) dy, 


0 


(16a) W= f Se [ fi tboycoskzdn) ap. 
0 0 

Da der Integrand nur von cos q, nicht aber anderweitig von g abhangt, kann 

die Integration auf das Intervall von 0 bis x beschrankt werden. Die Integration 

nach k laBt sich nun in den Integralen (16) und (16a) geschlossen ausfiihren. 

Es gilt: 


1 

r SS Ore 
I Jy (ko) sin cette | — 
| 0 wo>2, 

2 ah 

a 1 a Ose 
Pk tes ats | ot a) 
: es O>Z. 


Bei der Integration tiber y lauft w von |y —1| bisv+1. Fiir z<|v—4| erhalten 
wir also wieder die fritheren Ergebnisse: 


1) Ci Omur 


% 


Le ireeuly, 
2) Ww! [= 1808 F dey = 
aS 2 Lo i >A. 
0 


Fiir |1—7|<z<1+~7 erhalten wir mit w (gy) =z 


(17) Se “eee ao, 


sf ret 


Hu) wp? V2 = Oy 
0 


Po 


& VE ete ey ae 
oe u 


. 
(18) haa 
| 


* SzEcé [8], S. 423ff. 


92 HANSJORG OSER: 
Fiir 1-+7<z< 00 ist pyp=7 zu setzen. 
4. Die Charakteristiken z=1—7 und z=7—1. . 
Aus den Darstellungen (18) erhalten wir durch Grenzbetrachtung, wenn wir 
e=1—r+e baw. z=r—1-+¢ einsetzen und im Integranden nur Glieder erster 
Ordnung in é beriicksichtigen : ; 
hi Wi(r,2) =1-——, 


Won 
lim AG a) = oF 


g—r—-l+ 


Fiir U ergibt sich nach S. 85: 


lin OU@32)3 im Um) r= 


z—~1—-r+0 g—>r—1+ 


Aya 
2yr ° 


z-(r-1) — 
2 Be ¥ 5 


SS 


Abb. 2, Die Funktionen U und V in Abhangigkeit von z—(r—1) 


Nahert man sich also von oben her den Charakteristiken, so erhalt man als 
Steigung der Stromlinien: 


zZ=1—7+0: WiU=2)7 4, 

z=r—1+0: Wit 45 
Beim Uberschreiten der beiden Charakteristiken bleibt der Druck P stetig, 
dagegen erleidet seine Ableitung einen Sprung. 


2. Die Charakteristik z=1-++7. 


Alle Geschwindigkeitsfunktionen U, V,W werden hier unendlich und ver- 
halten sich in dieser Umgebung asymptotisch wie: 


1 
anyr 8 lz— (+1). 
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Man leitet sich dies leicht aus den Beziehungen des § 5 her. Der Druck ist auch 
hier wieder stetig, hat aber eine unendliche Ableitung, was man der Darstellung 


P(r, 2) = ne (1+7—2) {logy +1—2z) -—1}4--- 


20 
entnehmen kann. Durch die Punkte sind weitere regulare Glieder angedeutet. 


z-(r-1) _ 


0 7 


¥ 3d 


SSS 


Z-(r-) — 
3 


Abb. 4. Die Funktion P in Abhangigkeit von z—(r—1) 


Wir stellen also fest, daB in der Umgebung der Charakteristik z=1-+-7 die 
Voraussetzung kleiner Geschwindigkeiten nicht mehr erfiillt ist. Die Lineari- 
sierung der Differentialgleichungen ist also zumindest in dieser Umgebung nicht 
statthaft. Jedoch werden hier nicht nur die Tragheitsglieder, sondern vielmehr 
auch die Zahigkeitsglieder der allgemeinen Navier-Stokesschen Gleichungen 
wirksam werden. Am starksten ins Gewicht fallen werden diese Einfliisse in der 
Nachbarschaft des Punktes (0,1) — dem Schnittpunkt der Charakteristiken — 
weil hier zur logarithmischen Singularitat noch die algebraische y~? hinzutritt. 
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Obenstehend sind die Funktionen U, V, W und P graphisch dargestellt. 
Als Abszisse wurde z—(r—1) gewahlt. Auf diese Weise riicken die Unstetig- 


keiten in die Abszissen 0 und 2 (fiir 7>1). 


§ 7. Die Stromfunktion 
Fiir rotationssymmetrische Probleme wird die Stromfunktion folgender- 


maBen definiert : pee : ep 3 x ve 


Aus den bekannten Funktionen U und W ergeben sich die Richtungselemente 
der Stromlinien ¥Y =const sehr einfach zu 

dzidr = —4AY = Wt. 
Mit Hilfe dieser Beziehung wurde Y%=const graphisch aus den Richtungs- 
elementen aufgebaut. Die Abb. 5 eingezeichneten Pfeile entsprechen der momen- 


y= const 


7. 


7 


2 


re 


Abb. 5. Die Stromfunktion yp (yr, z) = const 


tanen Strémung bei Aufwartsbewegung der Platte. Die Stromlinienabstande 
sind willkiirlich gewahlt und bedeuten nicht, daB iiberall dieselben DurchfluB- 
mengen vorliegen. 

Offen ist noch die Frage, ob im Punkt (0,1) die Kontinuitat erfiillt ist, 
weil dort unendlich viele Stromlinien ein- und auslaufen. Um den Kontinuitits- 
beweis zu erbringen, wollen wir eine geschlossene Flache um den Punkt (0, 1) 
legen und den gesamten Durchflu8 durch diese Flache berechnen. Wir wahlen 
hierfiir den Zylinder, der die schwingende Platte als Grundflaiche und die Fliche 
y=1 als Mantel hat. Die Deckflache soll bei z= 00 liegen. Rechnen wit Ein- 
stromung negativ und Ausstrémung positiv, so erhalten wir: 


~~ 
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1) Grundflache 


= ff (-—W)df=—x, 
2) Mantelflache von z=0 bis z=2. 

; Q2.= Sf U df. 

Mit 
aah COS Pp 

2 — rae 
und 

df=rdidz 


ergibt sich, wenn man die Winkelintegration mit einem Faktor 22 beritick- 
sichtigt, sofort fiir 71: 
2 


Qo 2 | fa et dgiaz: 
: ae V2 —2 +4 2cosp is 


Man setzt cos m = 2 cos?/2 —1 und kann das innere Integral sogleich ausfiihren * 


0. 2 J (2E(«, 2/2) —F(a,2/2))dz  (k=sina). 


Hierbei sind E und F unvollstandige elliptische Integrale, die wegen des Argu- 
mentes z/2 in die vollstandigen elliptischen Integrale E(k?) und K(k?) iiber- 
gehen**. Dann wird (mit k=2/2) 

1 


ee [2E(k?) — K(k?)] dk, 


=) fli-9) dk? = 4k E(R2)) = 4. 


0 


3) Mantelflache von z=2 bis z= ov. 


= ee ee anda. 
Os= of ee . 


Hier wurde die Winkelintegration schon mit einem Faktor 2a beriicksichtigt. 
Man findet ganz entsprechend wie im Fall 2): 


Q=2 f z[B(e,, /2) — (1 — 8) F(a, 2/2)] — 


x} 


F(a, 2/2 )p ae. 


Hier ist k=sin «=2/z. Beim Ubergang zu den vollstandigen elliptischen Inte- 
gralen ergibt sich 


Qs=2 f {FZ (EW) — (1 — #) K(P)] — RK) dz, 


2 


1 
= 2 f C(k) de = — 4B (kf = — 442. 
0 


* MEVER ZUR CAPELLEN OA Se RSS ic 
** JAHNKE-EMDE [6], S. 55 
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4) Die Deckflache z= co liefert keine Einstromung, da dort W vest . 
Es. ist: ae 
Ques = Qr+ Qo+ Qs = 0. 


Damit ist die Frage, ob im Punkt (0, 1) Fltissigkeit verschwindet oder neu ent- — 
steht, beantwortet: Es flieBt genau so viel Fliissigkeit ein, wie von dort auslauft. 
Die von der Platte her in den Zylinder einstrémende Flissigkeit flieBt vom f 
Kegelmantel z=1—y7 wieder ab. Zwischen den Charakteristiken z=r7—1 und 
z—=yr+41 kénnen wir daher drei Bereiche unterscheiden. Rechnen wir alles auf — 
dem Zylindermantel y=1, dann flieBt zwischen z=0 und z=z, die von dea 
Platte herkommende Fliissigkeit. Zwischen z, und z, strémt die von ‘oben in 
die Kegelspitze flieBende Fliissigkeit. Zwischen z =z, und z= 2 ist der DurchfluB 
insgesamt null: alle Stromlinien biegen um und verlassen den Zylinder wieder. 
Da Q, in Abhangigkeit von z gegeben ist, kann z, und z, aus den Gleichungen ; 


Q2(k,) = 42, E() =x, 
O,(R)'=4k, ERS) = 4. 
bestimmt werden. Es ergibt sich 
F107 Oi Mins Zora eke) 


In Abb. 5 sind diese beiden Trennstromlinien gezeichnet. Die untere hat die — 
Steigung —1 im Punkt (0,1), die obere beginnt mit einer Steigung von 0.12 | 
(Naherungsrechnung). : 
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